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ii 
ABSTRACT 


Development of an Electromagnetic and Mechanical Model for a Resonant Column 


and Torsional Shear Testing Device for Soils 
by 


Inthuorn Sasanakul, Doctor of Philosophy 


Utah State University, 2005 


Major Professor: Dr. James A. Bay 
Department: Civil and Environmental Engineering 

Two weaknesses in resonant column (RC) and torsional shear (TS) tests are the 
back electromotive force (back emf) created in the typical coil-magnet system of RC/TS 
equipment and the nonuniform stress-strain that occurs over the radius of specimen in 
RC/TS tests. This study focused on the developing of methods dealing with these 
problems more precisely. Methods proposed in this study are the use of electromagnetic 
model to describe the back emf effect, and the stress integration approach accounting for 
the nonuniform stress-strain for soil. | 

The electromagnetic model presents the back emf effect in terms of the equipment 
spring stiffness (keq) and equipment viscous damping coefficient (Ceq). The model was 
verified by the measured model parameters. Variation of each model parameter was 
evaluated over a range of frequencies. In the RC test, the back emf effect increases the 


resonant frequency and damping, whereas, in the TS test, the model verified that the 
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effect on the torque measurement could be neglected. The analytical model shows a 
large effect of back emf on damping measurements. The model predicts equipment 
generated damping, Deq, as a function of resonant frequency and the polar moment of 
inertia of the system. The effect of back emf is critical for the low-strain RC test of soil 
with low resonant frequency because the Deg is high and material damping of soil itself is 
small. 

A stress integration approach was developed to account for the nonuniform stress- 
strain in RC/TS testing. This involved integration of an assumed soil stress-strain model 
to develop a theoretical torque (T) and rotation (8) relationship. The theoretical T-@ 
relationship was developed using either the closed form or numerical integrations for 
different soil models, and curve fitting was performed to match the theoretical T-O 
relationship to the measured T-O relationship. The model providing the best fit to the 
experimental data was selected to represent the stress-strain behavior of the soil. The 
equivalent radius approach was modified using the stress integration method to improve 
the damping measurement. The equivalent radius ratio curve based on damping was 
developed and used to calculate strains for damping curve. 


(297 pages) 
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CHAPTER 1 


INTRODUCTION 


1.1 BACKGROUND 

An important task in geotechnical earthquake engineering is the measurement of 
dynamic soil properties. Although there are a wide variety of field and laboratory 
methods to measure soil properties, each has different advantages, disadvantages, and 
limitations. Most field methods, particularly the wave propagation techniques, are 
limited to measure low-strain soil properties. These tests include reflection, refraction, 
spectral analysis of surface wave (SASW), crosshole, downhole (and uphole), and 
seismic cone. In laboratory methods, soil properties can be determined over a wide range 
of strain levels. However, laboratory tests are performed on a small volume of soil and 
the in situ stress condition is almost impossible to reproduce. Dynamic soil properties are 
also influenced by many factors, such as density, soil structure, age, stress and strain 
history, and cementation, all of which can be easily destroyed by sample disturbance. 
Cyclic triaxial and cyclic simple shear tests measure high strain behavior. Resonant 
column, and cyclic torsional shear tests measure low to intermediate strain behavior. 
Bender element tests measure low-strain behavior. Variations of field and laboratory 
measurements with strain level are summarized in Table 1.1. Selection of a test method 
for measuring dynamic soil properties requires careful consideration and understanding 


of the specific tests. 
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Table 1.1 Typical Strain Levels Associated with Field and Laboratory Measurements 


Strain Level (%) 10° | 107 | 10° 
Seismic Reflection Test 
Seismic Refraction Test 


— 
Q 
= 


— 
S, 


— 
© 


i SASW Test 
Field UET 
Measurement | >@!SiUe Crosshole Test 


[Seance 
Cries | [| >} 


Laboratory 
Measurement 


The two sil properties used to describe soil behavior under cyclic loading are 
stiffness and damping. Stiffness and damping of soils are typically presented in terms of 
shear modulus, G, and material damping, D. Shear modulus and damping of soils are 
strain dependent properties. When the strain level is below 107 %, soil behaves as a 
linear elastic material and the shear modulus and material damping remain constant with 
strain amplitude. The shear modulus of soil in the linear range is known as Gmax. For the 
range of strain between 10*-10°%, the soil behavior becomes nonlinear elastic and the 
shear modulus decreases and material damping increases as the strain increases. At strain 
levels above 107%, soil behaves nonlinearly and the properties tend to change not only 
with strain level but also with the loading cycles. The variations of G and D over the 
entire range of strains are typically presented by the plots of G/Gmax versus log y and D 
versus log y, where y is a shearing strain. In geotechnical earthquake engineering 


problems, the evaluation of dynamic soil properties is required from low to high strains. 
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As shown in Table 1.1, the resonant column (RC) and torsional shear (TS) tests are the 
most effective methods for measuring the dynamic properties of soils over a wide range 
of strains. 

Utah State University recently acquired a Stokoe-type resonant column and 
torsional shear (RC/TS) apparatus. The apparatus was developed in the late 1970’s at 
The University of Texas at Austin and has been improved over the intervening years. 

The equipment is fix-free type, with the bottom of the specimen fixed and torsional 
excitation applied to the top. Previous studies reported a problem that occurs in the 
RC/TS drive system. As the magnet moves through the coil, an electromotive force is 
induced in the opposite direction to the magnet motion. This force is known as back emf. 
The back emf effect causes an error in measuring the material damping of soils. The 
error can be significant because the damping of soil is very small. In the past, empirical 
approaches have been used is quantify the equipment-generated damping due to the back 
emf effect. Using this approach, the material damping of soil is simply corrected by 
subtracting the equipment damping from the measured values. A better approach for 
predicting the effects of back emf is the use of the electromagnetic model as one 
developed in this study. 

Determination of strain in RC/TS tests is complicated due to the non-uniform 
stress-strain that varies linearly with the radius in a soil specimen in torsion. Previous 
work (Chen and Stokoe 1979) used the equivalent radius approach to account for the non- . 
uniform strain problem. This approach is adequate when calculating strain at low to 


intermediate levels, however, this approach is less accurate when performing the tests at 
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4 
higher strains. In this study, a stress integration approach was developed to account for 


the non-uniform strain problem more precisely. 

The purpose of this study is to develop approaches to account for the two 
weaknesses in RC/TS tests. These weaknesses are: (1) the back emf effect created in coil 
and magnet system of RC/TS equipment, (2) the nonuniform stress-strain that occurs 
over the radius of specie: The results from this study provide a simple and reliable 
method for predicting the equipment generated damping. The stress integration approach 
was proposed to account for the nonuniform strain effect. These improvements increase 
the accuracy of measuring the nonlinear stress-strain behavior using the RC/TS testing 


device. 


1.2 OBJECTIVES 

The first objective was to develop si deona model to describe the 
RC/TS drive system. Electrical, magnetic, and mechanical properties of the RC/TS drive 
system were used to describe the back emf effect. There are four parameters in the 
model. These are: (1) the coils resistance, Re, (2) the coil inductance, Le, (3) the torque- 
current factor, kj, and (4) the back emf-rotational velocity factor, kp. 

The second objective was to develop independent measurements to measure each 
model parameter and the polar moment of inertia of the drive plate. To perform these 
measurements, calibration specimens were used. The calibration specimens have a wide 
range in resonant frequencies. The effect of frequency on the model parameters and polar 
moment of inertia of the drive plate was investigated using these specimens. In order to 


minimize the adverse effect of ambient vibrations, this testing was performed on an 
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5 
isolation table, which was designed and built as part of this study. The model parameters 


and drive plate calibration were also performed on a drive plate from The University of 
Texas at Austin. The results of model parameters and drive plate calibration from the 
USU drive plate and the UT drive plate are compared and discussed. 

The third objective was to verify the model prediction of the equipment generated 
damping by comparing model predictions with actual equipment generated damping. 
Open and closed circuit tests were performed using the calibration specimens to 
determine the actual equipment generated damping. The predicted equipment generated 
was calculated using the model and compared with the measured values. Again, the 
equipment generated damping for the USU drive plate and the UT drive plates are 
compared. 

The fourth objective was to develop the stress integration approach to account for 
the nonuniform stress-strain in the RC/TS testing. This approach is used to develop the 
theoretical relationship between torque (T) and rotation (6) for an assumed soil stress- 
strain model. The assumed soil models used in this study are the hyperbolic model, the 
modified hyperbolic model, and the Ramberg-Osgood model Either closed form or 
numerical integration is used to develop these T-0 relationships. Curve fitting techniques 
are used to match the measured T-6 relationship to the theoretical T-6 relationship. The 
model providing the best ft to the measured T-O relationship is used to develop the stress- 
strain relationship for a soil. The equivalent radius approach (Chen and Stokoe 1979) 
was modified using the stress integration approach to account for the nonuiform strains in 


damping. 
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The fifth objective was to apply the electromagnetic model and the stress 
integration approach in RC/TS testing of Bonneville clay. The RC/TS testing was 
performed on Bonneville clay to develop the modulus reduction and damping curves. In 
this test, the electromagnetic model was used to correct for equipment generated damping 
and the stress integration approach was used to calculate stresses and strains. The 
modulus and damping curves for Bonneville clay are compared with generic modulus 


reduction and damping curves. 


1.3. ORGANIZATION OF DISSERTATION 

Chapter 2 contains a literature review of RC/TS testing. The review is divided 
into four sections; the behavior of soil during cyclic loading, the normalized modulus 
reduction and damping curves, the nonuniform stress-strain in torsional testing, and the 
effect of equipment generated damping in RC/TS equipment. 

Chapter 3 contains an overview description of RC/TS test equipment and method 
of analysis. The basic analysis methods used for abe results for both RC and TS 
tests are presented. 

Chapter 4 presents the two basic calibration of the RC/TS equipment. These 
calibrations are the drive plate calibration and the deformation sensors calibration. The 
drive plate is calibrated to determine the polar moment of inertia of the drive plate. 
Results of drive plate calibration for the USU and UT are presented. Finite element 
method was used for the analysis of nonrigid body motion of the USU drive plate. 
Results of this analysis are discussed and compared with the results from calibrations. 


The deformation sensors are also calibrated. The deformation sensors used in RC/TS 
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testing are proximitors, accelerometers, and the linear voltage displacement transducers 
(LVDT). The calibration factors used for each sensor are summarized and discussed. 

Chapter 5 presents the electromagnetic model development. The basic concepts 
of electromagnetism are used to describe toque and the back electromotive force (back 
emf) generated in the RC/TS drive system. The development of the electromagnetic 
model is then presented. The behavior and effects of back emf in RC/TS tests are 
discussed. 

Chapter 6 presents the measurements of electromagnetic properties of the RC/TS 
drive system. These properties are the model parameters described the electromagnetic 
model. An independent measurement was made of each parameter. Factors affecting the 
value of each model parameters were correlated to the strength of drive plate magnets. 
These measurements and the evaluation of the model parameters were performed for both 
the USU drive plate and the UT drive plates. Differences and similarities of both drive 
plate are discussed. Calculations of the drive plate efficiency are also included. 

Chapter 7 presents the verification of the electromechanical model. Comparisons 
of the predicted and measured equipment damping are included. Limitations of the 
proposed model are discussed. An explanation of discrepancies in the literature between 
predictions of equipment generated damping in RC and TS tests is presented. 

Chapter 8 presents the development of stress integration approach used to account 
for the nonuniform stress-strain in RC/TS tests. Closed form and numerical procedures | 
for integrating stress-strain models are presented. Curve fitting technique used for the 


stress integration approach is discussed. An evaluation of the equivalent radius approach 
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8 
is presented. An improvement of the equivalent radius for determining strain in damping 


measurements is also presented. 

Chapter 9 presents RC/TS testing performed on Bonneville clay to develop 
modulus reduction curves and damping curves. Index properties of the soil specimens 
tested in this study are summarized. The RC/TS testing procedures are presented. The 
analysis procedures developed in this study are applied to the RC/TS testing of 
Bonneville clay. These analysis procedures are the correction of equipment generated 
damping and the application of stress integration approach are presented. Finally, the 
factor affecting the modulus reduction and damping of Bonneville clay are compared 
with generic modulus and damping curves. Results are discussed. 

Chapter 10 contains a summary and the conclusions of this work. 


Recommendations for further research are also presented. 
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CHAPTER 2 


LITERATURE REVIEW 


2.1 INTRODUCTION 

The literature reviews presented in this chapter are separated into four sections, 
(1) behavior of saab during cyclic loading, (2) normalized modulus reduction and 
damping curves, (3) nonuniform stress-strain effects in torsional testing, and (4) effects of 
equipment generated damping in RC/TS equipment. 

The behavior of soil during cyclic loading is first described in detail. Three 
nonlinear stress-strain soil models are presented. Basic aspects of these models are 
discussed. 

The normalized reduction and damping curves for soils developed empirically by 
the past research studies are summarized. Factors affecting these two curves are 
presented and discussed. 

The nonuniform stress-strain in torsional testing is presented. Approaches 
accounting for this problem are presented and discussed in details. 

Finally, the effects of equipment generated damping in RC/TS equipment are 
presented. A summary of research works on the equipment generated damping is 
presented and discussed. 

2.2 BEHAVIOR OF SOIL DURING 
CYCLIC LOADING 
When soil experiences strain in the medium range, approximately between 10° % 


and 10°! %, it behaves as an elastoplastic material. Shear modulus tends to decrease as the 
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shear strain increases and hysteresis cause energy dissipation during loading cycles. 
Damping is used to represent the energy absorbing properties of soils and the cyclic 
behavior of soils is commonly described using the theory of viscoelasticity. However, 
unlike viscoelastic materials, the energy dissipation in soil is nearly independent of 
loading frequency. A typical nonlinear hysteresis loop for soil is shown in Figure 2.1. 
The hysteresis loop representing the stress-strain path during cyclic loading is 
related to the monotonic loading stress-strain path called the backbone curve. The rules 
for constructing a hysteresis loop using the backbone curve are called the Masing rules. 
These rules are illustrated in Figure 2.1. Suppose loading is increased from point O to a 
level indicated by point A, with a stress reversal occuring at point A. The unloading path 
is obtained by enlarging the loading path or backbone curve, OA, by a factor of two. The 
curve is then folded and the initial point, O, is translated to the point of stress reversal. 
The unloading curve extends to point B, which is axis-symmetrical to point A. The 
reloading path is obtained by enlarging the backbone curve in the same manner as for the 
unloading curve. The reloading curve starts at point B and intersects the backbone curve 
at point A, which was the point of initial stress reversal. These unloading and reloading 
curves complete the hysteresis loop which represents the nonlinear stress-strain behavior 
and energy dissipating characteristics of the soil. 
The nonlinear stress-strain relation is usually described by the secant shear 
modulus, G. This parameter can be defined as the slope of a line connecting the origin 
and a point on the backbone curve. For nonlinear soil, the shear stress is expressed as a 


function of shear strain thus the secant modulus can be obtained from: 
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Shear Strain 


Figure 2.1 Nonlinear hysteresis loop 


cst, (2.1) 


The energy dissipated per cycle is proportional to the area enclosed by the 
hysteresis loop, AW. Damping ratio is defined as the ratio of the energy dissipated per 


cycle to the maximum stored energy, W. Thus the damping can be expressed as: 
po oN (2.2) 
For nonlinear soil behavior, the maximum stored energy, W, is assumed to be 


equal to the area of the triangle bounded by a straight line defining the secant modulus as 


shown in Figure 2.1 and obtained from: 


W=srf(r). (2.3) 
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The area of the hysteresis loop (AW) representing the energy loss per cycle can be 


obtained from: 
Y 
AW = fdv- w]. (2.4) 
0 
By substitution and rearranging, the damping ratio can be written as: 


2 ffo)dy 
2| è 
pes ait. (2.5) 
x| fy) 


The damping obtained from the Masing behavior is considered as the frequency 
independent and hysteretic in nature. There is another type of energy loss due to 
viscosity that can be modeled using a linear spring-dashpot model. The simplest and 
most widely used model is called Kevin model consisting of a spring representing the 
elastic soil properties and a dashpot representing the damping connected in parallel. The 
loss of energy or damping generated in this model increases linearly with increasing 
frequency regardless of the strain amplitude. However, the damping properties of soil are 
practically frequency independent. Therefore, the application of the spring-dashpot 
model to soils requires judgment and should be restricted to special cases (Ishihara 1996). 

As the basic derivations for the secant modulus and damping ratio (G and D) are 
established, it is obvious that these two parameters are strongly related to the backbone 
curve and expressed as a function of shear strain amplitude. The backbone curve is the 
key characteristic of the nonlinear stress-strain relationship for soil under cyclic loading. 


Two models have been commonly used to describe the backbone curve are the hyperbolic 
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model and a four-parameter model known as the Ramberg-Osgood model (R-O model). 


Recently, Darendeli has proposed a modified hyperbolic model as an improvement over 
the conventional hyperbolic model (Darendeli and Stokoe 1997, and Darendeli 2001). 


The basic aspects of these models are discussed in following sections. 


2.2.1 Hyperbolic Model 

The hyperbolic model is the most widely used model to described nonlinear 
stress-strain relationship of soils. The stress-strain curve is bound by two straight lines; 
the first line is a tangent at small strains and second line is a tangent at large strains, as 
illustrated in Figure 2.1. The slope of the first line represents the small-strain elastic 
shear modulus known as Gmax. The second line is horizontal asymptote indicating the 
upper limit of the shear stress or shear strength, tz. The strain level where these two lines 


intersect is called reference strain, ,, which defined as: 


ae (2.6) 


The hyperbolic model can be expressed in the following equation: 


r= Se (2.7) 
pe 
Yi: 
where Gmax = small-strain shear modulus. 


The normalized modulus, G/Gmax calculated using: 


ee (2.8) 


Reproduced with permission of the copyright owner. Further reproduction prohibited without permission. 


14 
can be substituted into Equation 2.5 to obtain the following equation for damping 


(Ishihara 1996): 


D=- E Ea = (2.9) 
Yr Y: 

It should be noted that the damping ratio calculated from Equation 2.9 has 
limitations for predicting damping at very small and very high strain amplitude. When 
the shear strain becomes infinitely large, the damping ratio tends to converge to the value 
of 2/n = 0.637. The calculated damping ratio deviates from the actual soil behavior at 
very high strain hence overestimating the damping ratio. On the other hand, at very small 
strain the damping ratio gradually decreases as the strain decreases and converges to zero. 
Numerous EERE on soils have shown that the small stain damping is usually 
higher than the damping calculated from Masing behavior, thus it underestimates the 
actual damping of soil at small strains. 

Only two parameters, Gmax and Yr, are required to define the hyperbolic model. In 
many cases, these two parameters are insufficient to characterize shear modulus and 
damping ratio of soil. As a result, numerous researchers have attempted to adjust and 
refine the model to achieve a good fit to experimental data for different soil types (e.g., 
Hardin and Drenevich 1972a, 1972b, Seed et al. 1986, Vucetic and Dobry 1991, 


Darendeli and Stokoe 1997) 


Reproduced with permission of the copyright owner. Further reproduction prohibited without permission. 


15 
2.2.2 Modified Hyperbolic Model 


The modified hyperbolic model has been proposed by Darendeli and Stokoe 
(1997). The two-parameter hyperbolic model was modified by simple approach in order 
to obtain the normalized reduction curves for all soil types. In addition to the two 
parameters (Gmax and y,) the parameter, a so called a curvature coefficient, a, was added 


into the conventional hyperbolic model as shown below: 


r= Smt (2.10) 


(7) 
Y: 


As discussed by Darendeli (2001), the advantage of the curvature coefficient, a is 
that the stress-strain curve can be calculated even if the strain level is far below failure 
thus the value of t¢ is not identified. In this approach, the reference strain, y, is assumed 
to be the strain amplitude at shear modulus equal to one half of Gmax (Yr = Yo/Gmax-0.5): 
This approach simplifies and provides better approximation for the modulus reduction 
curves obtained from experimental data. 

For the damping curve, Darendeli (2001) also proposed the modified model using 
two additional parameters (b and Dinin) as expressed as: 

o1 
D= {=| D masing + D min > (2.11) 


where D masing = Masing behavior damping described the damping of soils at high 
strain. 


It is noted that Dmasing for coefficient of curvature, a, equal to 1.0 is the damping 
obtained from conventional hyperbolic model and derived in similar method as discussed 


earlier, thus it can be calculated from Equation 2.9. In Derendeli approach, Dmasing for 


Reproduced with permission of the copyright owner. Further reproduction prohibited without permission. 


16 
any value of coefficient of curvature, a, is assumed to be a polynomial function of 


Dyasinga=1.0- The numerical integration was used to determine the polynomial function 
terms (c1, c2 and c3) using sixty damping curves of soils and they were assumed to be 
functions of curvature coefficient. Dmasing is determined from the expression below: 
D masing = CiD masing a=1.0 + C2D masing a=1.0 + C3 D masing a=1.0 (2.12) 
where cı =-1.1 143a? + 1.8618a + 0.2523, 
c2 = 0.0805a’ —0.0710a —0.0095, and 
c3 = -0.0005a° + 0.0002a + 0.0003. 

The additional two parameters are the scaling coefficient, b and the small-strain 
material damping ratio, Dmin. The scaling coefficient, b is used to adjust the damping at 
high strains where the Masing behavior tends to overestimate the damping value. The 
small-strain damping ratio, Dmin represents the damping of soil at small strains. This 
modified hyperbolic model proposed by Darendeli (2001) was called the four-parameter 
model. The parameters y,, a, b, and Dmin are used to describe the model. The empirical 
relationships for these parameters were developed using statistical analysis based on the 
experimental data of four different soil types; clean sands, sands with high fine content, 
silts, and clays (Darendeli 2001). The factors affecting these model parameters are 


confining pressure, soil plasticity, overconsolidation ratio, loading frequency, and 


number of loading cycles. The relationships are shown below: 


y, =[p, +9,(PD(OCR)” b*, (2.13) 
a=9;, (2.14) 
D =|, +9, (PI(OCR™ ) b* [1+ 9, In(frq)], (2.15) 
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b=@,,+9,InW), (2.16) 


where © © = mean effective confining pressure (atm), 
PI = soil plasticity index (%), 
OCR = overconsolidation ratio, 
frq = loading frequency, 
N — = number of loading cycles, and 
6, through 6:2 = model parameters tabulated in Table 2.1. 

Although the damping ratio at very high strains and very small strains has been 
adjusted in this model using the parameters, b and Din, the model does not accurately 
predict damping at small strains at very low loading frequencies. Neither hysteretic 
damping nor the viscous damping can explain the soil deine behavior for this special 
case. A complete model describing the actual damping behavior for soils has not yet 


been developed. 


Table 2.1 Mean Values and Variances of the Model Parameters after Darendeli (2001) 
Parameters 
C a En | _9.99E-07_| 
C & | 1018-03 4166-09 
C o | 3.25E-01_|_2.855-03_| 
| ġa | 3.48E-01 | 2.20E-04 | 
| gs | 9.19E-01 | 6.78E-05 | 
C o | 801E01 | 1.73E-03_| 
0 


3 
1.29E-02 3.82E-06 
-1.07E-01 2.49E-03 
1 


2.92E-01 | 7.66E-04 
bii 6.33E-01 | 2.23E-04 
-5.66E-03 | 5.02E-06 


i 
02 
c) 3 
1} 
{} 5 
D 
A es 
C o | 28o | 4968-04 
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` 2.2.3 Ramberg-Osgood Model 


The four-parameter model was originally proposed by Ramberg and Osgood in 
1943 and modified for soil (Streeter et al. 1974) into the form of: 


11 
T 


Ty 


Tatia 


Yy ty 


: (2.17) 


where Ty = shear stress, 
Yy = shear strain, 
a = constant, and 
r = constant. 
Idriss et al. (1978) applied the model to soil by defining ty and yy as the shear 
strength, te and reference strain, yr, respectively. The reference strain is the same as the 
reference strain used in the hyperbolic model. Making this substitution, Equation 2.17 


can be written as: 


eee) ey (2.18) 


t-l 
T 
i+a 


G 


max Yr 
Unlike the hyperbolic model, the Ramberg-Osgood model is calculates shear 
strain as a function of shear stress, y = g(t). The damping can be determined as the 


function of shear stress (Ishihara 1996) using: 


2 fedr 


Dase, 
G 0) 


(2.19) 


Substitution of Equation 2.18 into Equation 2.19 gives: 
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r-l 


T 
-1 (|G 
patty Smt] (2.20) 
qr +1 T 
i+a 
G max Yr 
where a= -1, and (2.21) 


Yr 
yt = shear strain at failure. 
If the shear modulus at failure, Gg, and damping at failure, Dg are known, the 


parameter, r can be determined from: 


nD 1 
|e sical eee 
2 1-9 
Sar aaah (2.22) 
1- max 


The four parameters of the Ramberg-Osgood modelare Gmax, Yr, &, and r. Values 
for the a and r parameters are usually obtained from curve fitting laboratory data. The 
empirical relationships correlated between these parameters and soil properties are not in 
the literature. 

2.3 NORMALIZED MODULUS REDUCTION 

AND DAMPING CURVES 

The nonlinear cyclic behavior of soil can be characterized using normalized 
modulus reduction curves (G/Gmax versus y) and damping curves (D versus y). Over the 
past three decades, several researchers have investigated the soil properties that effect the 
shape of these curves. The factors affecting these two curves are characterized into two 


categories; nature of the soil and loading conditions. The first category includes type of 
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soil, plasticity index, overconsolidation ratio (OCR), void ratio, degree of saturation, 


grain characteristics, mineralogy, and soil structure. The second category includes strain | 
amplitude, mean effective confining pressure, number of loading cycles, and loading 
frequency. 

The earliest work on the factors affecting the shear modulus and damping of soils 
was conducted by Hardin and Drenevich (1972a, 1972b). They were the first to propose 
the hyperbolic function for stress-strain relationship of soils. They also proposed an 
empirical relationship for predicting Gmax was also proposed as a function of void ratio, e, 
plasticity index, PI, overconsolidation ratio, OCR, and mean effective confining pressure, 
©. They developed approximate shapes of material damping curves using two 
additional parameters; the maximum damping ratio of soil, Dmax, and the hyperbolic 
strain, yn. These vary with soil type, mean effective confining pressure, Go , number of 
loading cycle and loading frequency. 

Seed proposed a series of normalized modulus reduction and damping curves for 
sand (Seed and Idriss 1970, Iwasaki et al. 1978, Seed et al. 1986). The curves shown in 
Figure 2.2 represent the culmination of that series (Seed et al. 1986). The upper and 
lower ranges in the curves are attributed to variability in the characteristics of the 
granular particles (shape, size, gradation, and mineralogy), variability of nonlinear soil 
behavior, effect of confining pressure and accuracy in measurement. These curves 
represent the behavior of clean silty sand with confining pressure ranging from 0.25 atm 


to 4 atm. 
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--- Upper Bound 
— Average 


— - Lower Bound 


0.0001 0.001 0.01 0.1 1 
Shear Strain, y (%) 


Figure 2.2 Empirical relationships for normalized modulus reduction and material 
damping for sand after Seed et al. (1986) 

Sun et al. (1988) proposed the curves for soils with plasticity as shown in Figure 
2.3. Their normalized modulus reduction curves account for the effect of plasticity on 
nonlinear soil behavior while material damping curves are presented in terms of one 
mean curve and an upper and lower bound for all soils within the total range of soil 


plasticities. 
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Sun et al. (1988) 
—— PI=5-10 
---- PI = 10-20 


~--— PI = 20-40 


Sun et al. (1988) 
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0.0001 0.001 0.01 0.1 1 
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Figure 2.3 Empirical relationships for normalized modulus reduction and material 
damping for clay after Sun et al. (1988) 

Idriss (1990) proposed two “average” normalized modulus reduction curves for 
sands and clays and a single material damping curve for all soil types. 

Vucetic and Dobry (1991) proposed normalized modulus reduction and material 
damping curves accounting for the effect of plasticity on both modulus and damping as 


shown in Figure 2.4. However, the values of damping at low strain have been left 
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undefined due to the lack of data. The data used by Vucetic and Dobry (1991) were 


obtained from soils tested at confining pressures ranging from 0.25 to 4 atm. 


Vucetic and Dobry (1991) 
—— Sand 
D (%) === PI =15 


oer. - — 
o TEO Ta a A 


0.0001 0.001 0.01 0.1 1 
Shear Strain, y (%) 


Figure 2.4 Empirical relationships for normalized modulus reduction and material 
damping after Vucetic and Dobry (1991) 
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24 
Ishibashi and Zhang (1993) proposed a set of equations generate normalized 


modulus reduction and material damping curves accounting for confining pressure, ©, 


and soil plasticity, PI. The equations proposed by Ishibashi and Zhang (1993) are: 


= =K(y, P5, , and (2.23) 
0 333(1+6 00") G 2 G l 
D=———*_—_— 30,586 — | -1.547 — | +1}, (2.24) 
2 G max G max 
where m(y,PI)-m, = oar - nn 000886) tee : 
Y 


0.492 
K(y,PI) = zi + ann p 2000102 +n) Ht and 
| y 


0.0 for PI=0 
3.37x10°PI'™ for 0 <PI<15 
7.00x107 PI for 15<PI<70 
2.70x10° PI" for PI>70 


nPI) = 


The normalized modulus reduction and damping curves generated by the equations 
proposed by Ishibashi and Zhang (1993) are presented in Figures 2.5-2.7. Significant 
effect of confining pressure are observed on sand and clay as shown in Figure 2.5 and 
2.6, respectively. The data used to develop these curves are the results of tests performed 
at confining pressure less than 10 atm. Unrealistic relationships are observed at high 
confining pressure as shown in Figures 2.5 and 2.6, where the normalized shear modulus 
value begins to exceed 1.0. Figure 2.7 shows the effect of PI on normalized modulus 


reduction and damping curve at a confining pressure of 1 atm. Unrealistic values (greater 
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than 1.0) of normalized shear modulus are obtained for PI = 15 and 30 in the intermediate 


strain level. More realistic values are obtained for PI = 50, 100, and 200. The equations 
proposed by Ishabashi and Zhang (1993) do not accurately represent the normalized shear 


modulus and damping of soils at high confining pressures or the range of PI’s from 0-40. 


1.2 
1.0 
0.8 
G/G max 
0.6 
0.4 
0.2 
0.0 
20 
Ishibashi and Zhang (1993) 
15 Sand Curves 
D (%) --- 0.25 atm 
— latm 
10 — -4atm 
5 
0 
0.0001 0.001 0.01 0.1 1 


Shear Strain, y (%) 


Figure 2.5 Empirical relationships for normalized modulus reduction and material 
damping for sand after Ishibashi and Zhang (1993) 
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Ishibashi and Zhang (1993) 
Clay curves (PI = 30) 
--- 0.25 atm 


— I atm 
— -4atm 


0.0001 0.001 0.01 0.1 1 
Shear Strain, y (%) 


Figure 2.6 Empirical relationships for normalized modulus reduction and material 
damping for clay after Ishibashi and Zhang (1993) 

As discussed in the previous section, Darendeli (2001) proposed a modified 
hyperbolic model to characterize normalized modulus reduction and material damping. 
Curves generated using Darendeli’ s equations are shown in Figures 2.8-2.10. Soil 
specimens used in Darendeli (2001) database were taken from Northern California, 


Southern California, South Carolina and Lotung, Taiwan. The curves proposed by 
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Ishibashi and Zhang (1993) 
—— Sand 
~--- PI=15 


0.0001 0.001 0.01 0.1 1 
Shear Strain, y (%) 


Figure 2.7 Empirical relationships for normalized modulus reduction and material 
damping show the effect of plasticity index after Ishibashi and Zhang (1993) 


Darendeli (2001) show larger effect of confining pressure than the curves proposed by 
Ishibashi and Zhang (1993). The material damping curves at small strain account for the 
effect of confining pressure. The effect of plasticity index is shown in Figure 2.10. 
Darendeli shows that the damping at low strains for soil with high plasticity tends to be 


higher than soils with low plasticity. 
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Sand Curves 
PI = 0, N = 10 cycles 
frq = 1 Hz, OCR= 1 


20 

15 Darendeli (2001) 
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10 — -4atm 


0.0001 0.001 0.01 0.1 1 
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Figure 2.8 Empirical relationships for normalized modulus reduction and material 
damping for sand after Darendeli (2001) 


2.4 NONUNIFORM STRESS-STRAIN IN 

TORSIONAL TESTING 

One of the weaknesses of resonant column and torsional shear method is the fact 
that stress-strain distribution is not uniform in the radial direction of the horizontal plane 
of a soil specimen during torsional loading. Theory of elasticity was utilized originally 


because the torsional shear test was developed to measure the small strain soil behavior. 
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As described earlier, soils behave nonlinearly when subjected to strain levels as small as 


10° %. The nonuniform stress-strain effect becomes significant when testing on soil at 
medium to high strain levels. Using hollow rather than solid specimen can minimize this 
effect however the hollow specimens are not commonly used because of trimming and 


handling difficulties. 


1.2 


1.0 


0.8 
G/G 
ar 6 Clay Curves 
` PI = 30, N = 10 cycles 


0.4 frq = 1 Hz, OCR =1 


0.2 


0.0 


Darendeli (2001) 
0.25 atm 

— latm 

— -4atm 


— - 16atm 


0.0001 0.001 0.01 0.1 1 
Shear Strain, y (%) 


Figure 2.9 Empirical relationships for normalized modulus reduction and material 
damping for clay after Darendeli (2001) 


Reproduced with permission of the copyright owner. Further reproduction prohibited without permission. 


30 


Effect of PI 
6,=1 atm, N = 10 cycles 
frq = 1 Hz, OCR=1 


Darendeli (2001) 


0.0001 0.001 0.01 0.1 1 
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Figure 2.10 Empirical relationships for normalized modulus reduction and material 
damping show the effect of plasticity after Darendeli (2001) 


In torsional testing, measurements are made of torque (T) and rotation angle (0). 
Several approaches have been proposed to correlate the measured T-O relationship to the 


stress-strain relationship of the soils (Hardin and Drenevich 1972a, Taylor and Parton 


1973, Chen and Stokoe 1979). 
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Hardin and Drnevich (1972a) recommended using the average stress and strain in 


the soil specimen. For a cylindrical soil specimen with the height of L, the average 


strain, Yay can be expressed as: 


fyda EF 
Ta =F -26(% =) (2.25) 


~ fda 3L\(R-R? 
where 6 = rotation angle, 
R, = outside diameter of the specimen, and 
R; = inside diameter of the specimen. 
In the case of solid specimen, Rj = 0, and the average shear strain occurs at a 
radius of 0.667 times the sample radius. Therefore, the average shear stress and strain are 


obtained from following equation: 


es ~0.667R (2.26) 


tx 0.667R, (2.27) 


P 
where T = torque, and 
I, = area moment of inertia of soil specimen. 

This approach is based on the summation of linear variation of strain over the 
radius of soil specimen but does not account for the nonlinear stress-strain soil behavior. 

The equivalent radius approach dealing with the nonuniform distribution of strain 
in soil specimen was developed by Chen and Stokoe (1979). The hyperbolic model and 
Ramberg-Osgood model were used as known stress-strain relationship for soils and used 


to calculated T-O relationship using the following equations: 
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o=, and (2.28) 


T= JerdA =2n |tr'dr. (2.29) 
A 


0 

Chen and Stokoe developed T-0 curves by dividing the cross section of a 
theoretical soil cylinder equally into N concentric rings. The cross-sectional area (A) of 
each ring was equal, therefore the width of the rings varied with the distance from the 
center of the soil cylinder. For a given rotation angle, the average shear strain on each 
ring was calculated using Equation 2.28. Using the hyperbolic or Ramberg-Osgood 
stress-strain relationship, an average shearing stress was determined on each ring. The 
total torque was obtained from Equation 2.29. This procedure was repeated for different 
rotation angles using the results of this analysis, normalized plots of T versus O were 


developed as shown in Figure 2.11. The normalized values of T, and 0,, are defined as: 


T, =JG,,,,9,, and (2.30) 
y,L 

6=— , 2.31 

et 2.31) 


where yr is the reference strain from Equation 2.6. 

Figure 2.11 shows three different curves; one generated from the hyperbolic 
model, the other two from the Ramberg-Osgood model. For the hyperbolic model, the 
parameters Gmax and yr are equal to 50,000 kPa, and 0.2%, respectively. The Ramberg- 
Osgood model is used for the other two curves use the parameters of a = 0.7, r= 5.5, and 


a = 1.2, r = 4.1, respectively. The values T and 0 from these computed normalized T-0 
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Figure 2.11 Normalized torque-twist relationship from Chen and Stokoe (1979) 


curves were used to calculate the corresponded effective shear modulus (Get) using a 
basic torsion equation. Each value of Ge was then matched with the value of G 
generated by the theoretical stress-strain relationship. The corresponded 0 from Ge and 
y from G were compared and a correction is applied to radius of specimen to achieve the 
best match of y from experimental and theoretical relationships. This correction is 


expressed in term of the equivalent radius ratio (Req) or the ratio of the equivalent radius, 


Teq and the radius of specimen, R. 
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Figure 2.12 shows plots of the Reg with normalized 6 generated by Chen and 


Stokoe (1979). The value of Reg is as high as 1.0 for R-O model case and 0.8 for 
hyperbolic model at small strain. The value of Req decreases to about 0.75 for very large 
strain for both models. Based on this study, the Reg is recommended to be varied from 
0.82 for shearing strain less than 0.001% to 0.79 for shearing strain equals to 0.1%. 
Therefore, the nonuniformity of strain can be minimized by a simple correction of the 
radius of soil specimen. This approach is considered reasonable in practice however the 


calculation of strain for damping is still in question. 
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Figure 2.12 Variation of equivalent radius ratio with normalized twist in solid specimen 
from Chen and Stokoe (1979) 
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2.5 EFFECT OF EQUIPMENT GENERATED 
DAMPING IN RC/TS EQUIPMENT 


Numerous laboratory experiments have shown that the material damping of soils 
increases with the strain amplitude. The damping values range from below 1% at very 
low strains to above 20% at very high strains. Accurate measurements of material 
damping are challenging, especially at low strains where the soil damping is very small. 
The resonant column and torsional shear methods are typical used to measure the 
material damping of soils. There are several potential sources of error in resonant 
column and torsional shear damping measurements and many researchers have proposed 
some procedures to mitigate these problems. 

One of the major problems is related to the equipment generated additional 
damping due to the drive system. Typical drive systems in the resonant column and 
torsional shear device use a coil and magnet system. The additional damping is created 
from the movement of the magnet in the coil inducing an electromotive force (emf) 
opposing the driving motion called back emf. Kim (1991) and Hwang (1997) performed 
calibrations to determine the magnitude of this equipment generated damping and 
developed correction procedure for resonant column and torsional shear tests. 

Kim’s and Hwang’s ealibeasons were made by performing the tests on metal 
specimens that have resonant frequencies range from 14 Hz to 150 Hz. The method is 
based on the assumption that the material damping of the metal specimens is very small, 
and thus can be neglected. Results from the resonant column calibration tests indicate 
that the equipment generated damping can be as large as 4% at low resonant frequencies 


and decreases to as low as 0.3% at high resonant frequencies. Figure 2.13a shows the 


Reproduced with permission of the copyright owner. Further reproduction prohibited without permission. 


plot of the measured equipment damping versus different resonant frequency (Hwang 


1997). The material damping of soil specimen at the same resonant frequency is 


corrected by subtracting the equipment damping at the same frequency (Hwang 1997). 
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Figure 2.13 Equipment generated damping ratio using metal specimens from Hwang 
(1997) 
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37 
In the torsional shear test, the equipment generated damping was measured at the 


loading frequencies from 0.2 Hz to 20 Hz. The results show a trend opposite to the 
resonant column results, with the damping increasing with frequency as shown in Figure 
2.13b. 

Recently, other researchers have proposed alternatives to these empirical 
corrections for the resonant column method. Li et al. (1998) proposed an alternative 
testing method using the energy-injecting virtual mass resonant alanis system. The 
main idea of this approach is to maintain the resonant frequency and damping of testing 
material at a pre-specified values using a feedback system, thus minimal the back emf 
effect. The new system was verified by performing the test on the aluminum rod with 
known properties. This system works reasonably well however it is complicated and 
required a special control system. 

Cascante et al. (2003) presented a current-based measurement as an alternative to 
the conventional voltage-based measurement to minimized the back emf effect in 
resonant shea test. The derivation of new transfer functions for voltage and current 


measurement were introduced and expressed as: 


ee r | 
£ a (alo ) 
T @- A) = $e, (2.32) 
V @) oe - 
T Rade (nto +d) ai [eneg 
Bir r 
DA 
AW) _ i 
se ae (2.33) 
I Ko) (1-(a/o, )? + (ct +d a e)” 


o 


Reproduced with permission of the copyright owner. Further reproduction prohibited without permission. 


38 
where © = angular frequency, 


©o = resonant angular frequency, 

i = imaginary unit, 

A(q) = spectral density of the acceleration, 

V1(@) = spectral density of the total voltage, 

I(@) = spectral density of the current, 

Im = distance from the center of the specimen to the magnet, 

Ta = distance from the center of the specimen to the accelerometer, 
J = rotational mass moment of inertia of the specimen and drive plate, 
Bl = magnetic-force factor of the coils, 

cg = damping coefficient due to the eddy-current forces, 

d, = viscous damping of the specimen, 

R = resistance of the coil, and 

L = inductance of the coil. 

This method was evaluated by testing three aluminum specimens and comparing 
the results between the current and voltage measurement. Results showed that the 
damping measured by the voltage measurement was as much as 50 times higher than that 
measured using current measurements. The damping measurement in dry sand using the 
voltage measurement was also observed to be up to 300% of the damping measured using 
the current measurement. Cascante et al. (2003) recommended using the current-base 
measurement for the resonant column test and analyzing results using the current transfer 


function. This study used the curve fitting processes to obtain the model parameters. 
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Meng and Rix (2003) developed an electromechanical model to derive the 


equipment generate damping which is expressed as: 


c 2k k 1 
D =A= it (2.34) 
a c R +R 2Jo 


where Deq = equipment generated damping ratio, 

Ceq = equipment generated damping term, 

ce = critcal damping, 

ki = torque against current calibration factor, 

kp = induced voltage against angular velocity, 

R, = resistance of source, 

R, = resistance of coil, 

J = total mass polar moment of inertia of the drive system and specimen, 
and @, = resonant angular frequency. 

Meng and Rix (2003) also proposed a solution to reduce equipment generated 
damping in resonant column testing using an instrumentation operational amplifier circuit 
to convert the voltage-mode source into a current-mode source. From this study, the 
equipment generated damping is reduced by more than three orders of magnitude 
compared with the voltage-mode source. However, the current limitations of the 
operational amplifier will limit testing at high strain levels. 

Wang et al. (2003) proposed the electromechanical model used for presenting the 
back emf effect in terms of the “virtual inertia,” Igy and the additional electrical 
damping, Gini The electromechanical model derived by Wang et al. (2003) was 


presented in form of: 
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(I, — Tene) O+ (Coot + Cour) 0+ K,0 =T e", (2.35) 
apBL 
Ter = Reso’ and (2.36) 
aBR 
Cage =, 2.37 
EMF R? ies oL ( ) 
where I, = rotational moment of inertia of the system, 


& = angular displacement, 

Csoit = lumped viscosity for the soil, 

K, = torsional stiffness, 

T,e'' = applied torque, and 

a and B = proportionality constants. 
Development of the electromechanical model proposed by Wang et al. (2003) is very 
similar to model proposed in this study, however, the effect of back emf was presented 
using a different formulation. Frequency effects for a and B were not investigated by 


Wang et al. (2003). 
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CHAPTER 3 


TEST EQUIPMENT AND METHOD OF ANALYSIS 


3.1 INTRODUCTION 

A Stokoe-type resonant column and torsional shear (RC/TS) apparatus was 
employed in this study. This equipment is one of the most widely used to perform the 
dynamic testing on soils for research and commercial purposes. The equipment has been 
developed and modified to operate both resonant column (RC) and torsional shear (TS) 
testing (Isenhower 1979, Lodde 1982, Ni 1987, Hwang 1997, Darendeli 2001). A 
detailed description of the RC/TS testing method and equipment are presented in this 
chapter. The theoretical basis used for analyzing results for both RC and TS tests are 
included. 
3.2 OVERVIEW OF RESONANT COLUMN 

AND TORSIONAL SHEAR TESTS 

There are several versions of the resonant column test apparatus employing 
different conditions to constrain the top and bottom of specimen. The fixed-free type, 
where the soil specimen is fixed in place at the bottom and driving force is applied at the 
top, is the method most commonly used. The vibration response is monitored at the top 
of specimen. The ideal fixed-free configuration is shown Figure 3.1a. For actual testing, 
a torque must be applied to the top of soil specimen, thus a drive plate as well as 
deformation sensors must be attached to the specimen at the top. As a result, the top of 
soil specimen is not in the ideal free condition. The approach dealing with this 


complication is to lump the attachments into one mass and the boundary condition for the 
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a) Fixed Free Type-No Mass b) Fixed Free Type-Added Mass 


Figure 3.1 Fixed-free configurations for RC/TS test 


ideal fix-free configuration must be modified. Consequently, the vibration response of a 
soil specimen in RC/TS test is modeled as a cylindrical rod with one end fixed and an 
attached rigid mass on the other end as shown in Figure 3.1b. 

Two different types of testing can be performed on the same soil specimen in 
RC/TS apparatus. The first type of test is the resonant column (RC) test. The general 
principle employed in the RC test is to excite the soil specimen with a steady-state 
torsional motion over a range of frequencies to identify the first mode resonant 
frequency. The shear modulus can then be evaluated utilizing the well-defined boundary 
conditions and the specimen geometry and mass. Material damping is determined by one 
of two methods; the free vibration decay method or the half power bandwidth method. 

The second type of test is the torsional shear (TS) test. The difference between the 


TS test and the RC test is mainly in the excitation frequency. In the TS test, a slow cyclic 
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loading in the range of 0.01 to 10 Hz is applied to the specimen. Shear modulus and 


damping are determined based on the characteristic of the hysteresis loop. 
3.3. OVERVIEW OF RESONANT COLUMN 

AND TORSIONAL SHEAR EQUIPMENT 

Figure 3.2 shows a drawing of the RC/TS equipment. The RC/TS equipment is 
capable of testing a solid or hollow cylindrical specimen and a variety of specimen sizes 
by using interchangeable top platens and base pedestals. The equipment consists of three 
basic components. First, the drive system used to apply a torsional cyclic load at the top 
of soil specimen. Second, the deformation monitoring system for measure the rotational 
and vertical displacements of the top of the specimen. Third, the confinement system 
used to apply isotropic confining pressure to specimen during consolidation and testing. 
RC/TS testing is performed on a specially designed isolation table to minimize the effects 
of ambient vibration on the measurements. A schematic diagram of the RC/TS drive 
system and deformation monitoring system is presented in Figure 3.3. Details of each 


component are described as below. 


3.3.1 The RC/TS Drive System 

The drive system in RC/TS equipment consists of a drive plate and eight drive 
coils. The drive plate is a four-armed aluminum plate that has rectangular shaped 
magnets attached to the end of each arm as shown in Figure 3.2. The drive coil is a long 
copper wire wound in multiple elliptical loops. The drive plate is positioned such that 
each end of magnet is at the center of a coil loop. When current passes through the coils, 


the interaction between coil and magnet behaves similarly to an electrical motor creating 
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Figure 3.2 The RC/TS equipment configuration 
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Figure 3.3 Drive system and deformation monitoring diagram for RC/TS equipment 
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torsion on the drive plate. More details on coil and magnet interaction are presented in 
Chapter 4. 

A harmonic sinusoidal function is sent to the drive coils from a function 
generator. For the RC test, a constant amplitude signal is generated over a range of 
frequency by a dynamic signal analyzer (Hewlett Packard Model 35670A). For the TS 
test, the signal frequency, number of cycles, and amplitude are controlled by a function 
generator (Hewlett Packard Model 3314A). For both tests, the signal is transmitted 
through a power amplifier before being sent to the drive coils. The power amplifier used 
for this study is AE Techron Linear Amplifier (Model LVC 608). There are two 
advantages to using a power amplifier. First, the signal can be amplified up to about 25 
times when performing a high amplitude tests. Second, the power amplifier has special 
features for monitoring voltage and current it outputs. This feature is useful for torque 
calculations in TS test. The voltage and current signals monitored by the power amplifier 


are then sent to the dynamic signal analyzer to be recorded and analyzed. 


3.3.2 The RC/TS Deformation Monitoring System 

The second component of the RC/TS apparatus is the deformation monitoring 
system which includes the rotational and vertical displacement monitoring devices. The 
rotational monitoring device consists of both accelerometer and proximitors. 

In the RC test, the accelerometer is used for monitoring the rotation of the drive 
plate. Typically, the RC test is performed over frequency ranges of 10-400 Hz. An 
accelerometers perform well in this frequency range. The accelerometer is mounted on 


the drive plate so it senses the motion of the drive plate. The signal from the 
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accelerometer is transmitted through a charge amplifier for amplification and signal 


conditioning. The output signal from the charge amplifier is then sent to the dynamic 
signal analyzer. The real time record is converted into frequency domain by Fast Fourier 
Transform (FFT) and the frequency response is displayed on the dynamic signal analyzer. 

For the TS test, torque is applied at much lower frequencies than those used in RC 
testing (0.01-10Hz). Accelerometers are not suitable for this low frequency range so 
displacement transducers known as a proximitor is used to monitor rotation. The system 
consists of two proximity probes and a metal target as shown in Figure 3.2. The two 
proximity probes are fixed to a rigid post and the metal targets are mounted on top of the 
drive plate. Two proximity probes are used in the system to minimize the effect of 
bending on rotation measurements. The output time of the proximitor system is 
monitored using the dynamic signal analyzer. 

A linear vertical displacement transducer (LVDT) is used as a vertical 
displacement monitoring device for RC/TS equipment. The LVDT monitors the 
settlement when soil specimen as it consolidates. The LVDT is rigidly mounted to the 
same post used for the proximitor system and positioned at the center of the dive plate as 
shown in Figure 3.2. The output from LVDT is recorded manually using the 
multipurpose circuit box explained below. 

The multipurpose circuit box was designed and built for two different purposes. 
The first purpose of the box is to condition the proximitor signals. The box has a . 
subtraction circuit used to subtract a predetermined proximitor center voltage from the 
output voltage from each proximitor. The signal out of the subtraction circuit can be used 


for adjusting the gap between the proximitor and target. The gap distance is adjusted 
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until the circuit output is zero volts. The proximitor is then at the center of its linear 


range. During the test, a signal proportional to the difference between the two proximitor 
signals is generated using an instrumentation quality differential amplifiers. The signal 
from this amplifier is sent to the dynamic signal analyzer. The subtraction circuit for the 
proximitor signals developed in this study is presented in Figure 3.4. The second purpose 
of the box is to excite and monitor the LVDT. The LVDT output voltage is displayed 


using a 4 % digit voltage monitor. 


3.3.3 The RC/TS Confinement System 

The third component of RC/TS equipment is the confinement system. As shown 
in Figure 3.2b, the RC/TS equipment is enclosed inside a cylindrical cell where the 
isotropic pressure is applied to a soil specimen. Air pressure is applied inside the cell 
through the air line located at the bottom plate of the cell. The pressure is controlled and 
maintained by an air regulator. The cell is made from stainless steel and has the capacity 
to withstand pressure up to at least 150 psi. Inside the cell, the specimen is surrounded 
by a fluid bath which minimizes air defusing into the specimen. The fluid bath is a 
cylindrical plastic sleeve fixed to the base pedestal and filled with fluid such as silicone 
oil. During consolidation, soil specimen is drained using the two drainage lines located 


in the bottom pedestal. 


3.3.4 Isolation Table 
An isolation table was designed and constructed in this study to accommodate the 
RC/TS equipment as well as to isolate the measurement system from building vibrations. 


The table is isolated vertically and horizontally by a combination of cables suspension 
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and an air spring system. The isolation is triangular in shape and suspended from ceiling 


beams by three cables with the lengths of about 8 ft. The table is made from steel and 
filled with mill scale. Total weight of the isolation table is about 3000 lb. Three air 
springs are used to support the three corners of the tabletop. Air pressure in each air 
spring is controlled by valve system so that the table remains level. The isolation table 
was designed to have a natural frequency in vertical displacement of about 1 Hz and in 
horizontal displacement of about 0.3 Hz. Figures 3.5 and 3.6 show a photograph and 


drawing of the isolation table. 


Figure 3.5 Photograph of the isolation table 
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Figure 3.6 Drawing ofthe isolation table 


Reproduced with permission of the copyright owner. Further reproduction prohibited without permission. 


52 
3.4 METHOD OF ANALYSIS 


Different analysis methods are required for RC and TS testing. The RC test 
analysis procedure is based on the wave propagation in elastic media. The TS test 
analysis assumes static behavior to measure is the hysteretic behavior of elastroplastic 
material. Both tests are complicated by nonuniform stress- strain in the soil specimen. 
The nonuniform stress- strain effect is taken into account more precisely using a stress 
integration approach developed in this study. This approach is presented in Chapter 8. In 


this section, the basic analysis methods used in RC and TS testing is presented. 


3.4.1 The RC Test Analysis 

In RC test, a cylindrical soil specimen is excited torsionally in its first resonant 
mode. The wave velocity is determined from the resonant frequency and from 
dimensions and mass of the soil specimen. The test configuration is the fixed-free type 
with added mass applied on top of the specimen as shown in Figure 3.1b. Torque is 
exerted on top of soil specimen is equal to the inertia force of the added mass. This 


torque can be expressed as: 


00 0°60 
Le a (3.1) 
where G = shear modulus of specimen, 


I, = area moment of inertia of specimen, 


4 
€ z for a solid cylindrical specimen with radius, R), 


J, = mass polar moment of inertia of added mass, and 


6 = rotational angle. 
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By applying the fix-free end condition, the governing equation for torsional vibration of 


soil column can be written as: 


Ppa ae ad (3.2) 
J, V; Ng 
where ©n = undamped natural circular velocity (=27f,), 


f 
fa = undamped natural frequency | = ——=—— |, 
r ASY l V1—2D? 


fm = resonant frequency for constant force amplitude, 


J = mass polar moment of inertia of specimen, 
(= sme’ for a solid cylindrical specimen with mass, m),and 


D = damping ratio of specimen, 
vs = shear wave velocity of specimen, and 
L = length of specimen. 
The right hand side of Equation 3.2 can be written in the form of BtanB, where 
B = @nL/vs. By applying the predetermined values of the polar moment of inertia of 
specimen and added mass, the value of B can be solved graphically by plotting a curve of 
B versus J/J,. Once an appropriated value of B is selected, the shear wave velocity of 


specimen can be determined from: 


v, s^, (3.3) 


Shear modulus (G) is then calculated using the following equation: 


G=pv. 


s 3? 


(3.4) 
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where p = mass density of the specimen. 

The resonant column test produces the response curve that provides the 
information regarding the dynamic properties of the test specimen. The response curve is 
the transfer function presented in term of the ratio of the displacement (or rotation) and 
force (or torque) amplitude. For a single degree of freedom system, the transfer function 


is expressed as: 


H = — = -~ and (3.5) 
T 2 2 
O, O, 
2D — 
tang = —— : (3.6) 
Se ores 
where H = transfer function amplitude, and 


ọỌ = phase angle between torque and rotation. 
For constant force amplitude, the frequency at maximum amplitude, fin is given 


by: 


f, =f,v1-2D’. (3.7) 


and maximum transfer function amplitude, Hmax is given by: 


Haa = l x x for small D (<10%). (3.8) 


2DV1-—D? 


Utilizing the shape of the response curve, damping ratios can be determined using the 


half-power bandwidth method. The method is applied by selecting two points on the 
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response curve at H is equal to 33 Hmax. These points are called fı and fz as presented in 


Figure 3.7. The amplitude at fı and fy can be expressed as: 


1 1 1 1 
=H e 3.9 
a= ED F e re Pa 
1232 2 
1+( +| 4D 
fa fn 
Squaring and rearranging Equation 3.9 gives: 
fay eT 
e) -22 (-2D?’)+(1-8D°)=0. (3.10) 


Roots of Equation 3.10 become: 


2 
(=) =(1-2D?)+2Dv1-D? . (3.11) 


f, 


The second-order term of D is negligible since the D is small. Equation 3.11 reduces to: 


fY 
a #1+2D. (3.12) 


n 


Subtracting the roots from Equation 3.12 gives: 
2 2 
ETSEN KAD: (3.13) 
f, f, 


Dz (£, -f KE, +f.) i 
4f, 


Rearranging gives: 


(3.14) 


Assuming that the response curve is symmetric, (f2+f1) is equal to 2f,. Equation 3.14 


reduces to: 


(3.15) 
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Figure 3.7 Typical response curve for soil from the RC test 


Equation 3.15 is widely used to determine the damping ratio of soil in the 
resonant column test. It should be noted that the half power bandwidth method is based 
on the assumptions that D is small (less than 10%) and the response curve is symmetrical. 
However, when performing the resonant column test at high strain level, the response 
curve becomes distorted and is no longer has a symmetrical shape. Many researchers 
recommended using the free vibration decay method to evaluate damping of soil for high 
strain test (Kim 1991, Hwang 1997). Details regarding to the free vibration decay 
method are briefly described below. 

Free vibration decay is the method used to determine damping utilizing the decay 
of vibration. One of the properties of viscous damping is that the magnitude of the decay 
of vibration of any two successive peaks is a constant ratio. The natural logarithm of two 


successive amplitude of motion is called the logarithmic decrement, 5 expressed as: 


S=In—, (3.16) 
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where Zn and Z,+1 = two successive peak amplitudes. 


The ratio of peak amplitudes is given by: 


Za = exp(f, D25). (3.17) 


n+l m 
Substitution of Equation 3.7 gives: 


2nD 


Vi-D? 


Therefore, Equation 3.18 can be used for calculating the material damping, D in the free 


b= 


(3.18) 


vibration decay method. The weakness of the free vibration decay method is that there is 
an uncertainty in determining the associated strain level. Typically, the strain level is 
estimated from the first three cycles of the free vibration decay curve (Hwang 1997) and 
this method is used at higher strain level (> 0.01%) when the half power bandwidth is not 


applicable. A typical free vibration decay curve is shown in Figure 3.8. 


g 
D 
-2x107 
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Figure 3.8 Typical free vibration decay curve for soil 
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3.4.2 The TS Test Analysis 


In the TS test, torque is applied to soil specimen at a very low frequency where 
inertia is insignificant. Current flowing through the coils is measured and converted into 
torque using calibration factor. In this study, this calibration factor is called torque- 
current factor, k;. More details on the measurement of the kj value are included in 
Chapter 6. Rotation is measured by the proximitor system and calculated using a 
calibration factor presented in Chapter 4. The relationship between torque and rotation is 
then transformed into the stress-strain hysteresis loop as shown in Figure 3.9. As 
discussed in Chapter 2, the nonuniform stress-strain effect occurs in both the RC and TS 
tests. For TS test, the nonuniform stress-strain effect has an impact on the determination 
of both shear stress and shear strain. The equivalent radius approach (Kim 1991, Hwang 
1997) has been used to calculate stress and strain that correspond to the measured torque 
and rotation. A more precise method called the stress integration approach will be 
presented in Chapter 8. 

When the stress-strain (t—y) hysteresis loop as illustrated in Figure 3.9 is 
developed, shear modulus is simply calculated from the secant slope of any point on the 


loop. This secant modulus can be expressed as: 


G=2, (3.19) 
Y 
The damping in the torsional shear test is determined from the characteristics of 


the hysteresis loop. As described in Chapter 2, the damping is the ratio of the energy loss 


per cycle to the maximum stored energy. An example of hysteresis loop is illustrated in 
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Figure 3.9 Hysteresis loop measured in the TS test 


Figure 3.9. The area of triangle represents the maximum stored energy. The area of 


hysteresis loop represents the energy loss per cycle. The damping is then calculated from: 


D= Tar ; (3.20) 
4r A, 
where A, = area of hysteresis loop, and 


Ar = area of triangle. 


It should be noted that damping could also be determined from torque-rotation (T- 
0) hysteresis loop. The damping from T-6 relation is calculated using the same concept 
as the damping from t—y relation therefore Equation 3.20 is also applied in both cases. 


Differences between two methods are discussed in Chapter 8. 
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CHAPTER 4 


DRIVE PLATE AND DEFORMATION SENSORS CALIBRATION 


4.1 INTRODUCTION 

Like any other laboratory testing equipment, the RC/TS equipment needs to be 
calibration prior to performing tests. This chapter presents the basic calibrations of the 
RC/TS equipment. Two basic calibrations are presented in this chapter. The first is the 
drive plate calibration to determine the mass polar moment of inertia of the drive plate. 
The second is the deformation sensors calibration to determine the calibration factor for 
each deformation sensor. 

The mass polar moment of inertia of the drive plate is an important parameter in 
the analysis of RC tests. In this chapter, an experimental measurement was used to 
calibrate the polar moment of inertia of the drive plate. The calibration was performed 
for the USU drive plate and the drive plate No. 9 at The University of Texas at Austin 
(UT). Calibration of results for each drive plate are compared and discussed. 

The calibration specimens are used to measure the polar moment of inertia of a 
drive plate. Details on the calibration specimens are also included in this chapter. 

The finite element method was used for an analysis of nonrigid body motion of 
the USU drive plate. Results of this analysis are discussed and compared with the results 
from the calibration. 

Deformation sensors are used for monitoring the torsional rotation and vertical 


deformation of soil specimens in the RC/TS apparatus. These sensors are proximitors, an 
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accelerometer, and a LVDT. In this study, the proximitors and the LVDT were 


calibration. The calibration factor of the accelerometer provided by the manufacturer 
was used in this study. The rotation measured from the proximitors and the 


accelerometer is compared and discussed. 


4.2 CALIBRATION SPECIMENS 

Calibration specimens were designed and built for the drive plate calibration and 
for measurement of electromagnetic properties of the drive system presented in Chapter 
6. Requirements for the calibration specimens as well as their design and construction 
are presented. 

The calibration specimens have three requirements; (1) the specimens can be 
tested repeatedly without changing properties, (2) each specimen is designed to have a 
specific resonant frequency, and (3) the specimens must have linear frequency 
independent stress-strain behavior. In this study, the frequency range of interest is from 
0-200 Hz. Four metal specimens were designed to have resonant frequencies at about 14 
Hz, 62 Hz, 128 Hz, and 188 Hz. Each specimen was custom fabricated at Civil and 
Environmental Engineering Department (CEE) machine shop, Utah State University. 

The specimen with resonant frequency of 14 Hz was built from a solid steel rod. 
Figure 4.1a shows a drawing of this specimen. The top and bottom of the rod were 
rigidly welded to circular steel plates. There are four threaded holes at each top and 
bottom plates so the specimen could be screwed tightly to the drive plate and the bottom 


platen of RC/TS apparatus. Although the top plate is rigidly welded to the steel rod, the 
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Figure 4.1 Drawing and photo of 14-Hz calibration specimen 
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weld is relatively weak because the steel rod diameter is only 0.189 in. To protect the 


welded joints, and to prevent the rod from bending during set-up, a support system was 
built for this specimen. The supporting system was built from a hollow, steel, threaded 
rod with a top platen attached to it as shown in Figure 4. la. This top platen can be 
adjusted up and down. During setup, the top platen is adjusted to be against the top plate 
of the specimen. The supporting system during setup is shown in left photograph of 
Figure 4.1b. Once the installation is finished the supporting top platen is lowered so it 
does not interfered with the specimen during testing. The photograph on the right in 
Figure 4.1b shows the specimen of 14 Hz in ‘esting position. | 

The 62 Hz, 128 Hz, and 188 Hz calibration specimens were made from brass tube 
with different diameter as shown in Figures 4.2-4.4. The top and bottom plate of these 
specimens are built from circular brass plate connected to the tube by welding. There are 
four threaded screw holes at the top and bottom similar to the specimen of 14 Hz. Size 
and dimension of these specimens are presented in Figures 4.2b, 4.3b, and 4.4b. 

For accuracy in determining the mass polar moment of inertia of the calibration 
specimens, each part was precisely machined in geometric shape and weighed before 
assembly. Dimensions were carefully measured. The mass polar moment of inertia of 
each part was calculated accounting for all screw holes. These values are presented in 


Table 4.1. 


4.3 DRIVE PLATE CALIBRATION 
As explained in Chapter 3, the analysis method for the RC test is based on the — 


theory of wave propagation in a finite rod under the fixed-free conditions. In RC/TS 
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Figure 4.2 Drawing and photo of 62-Hz calibration specimen 
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Figure 4.3 Drawing and photo of 128-Hz calibration specimen 
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Figure 4.4 Drawing and photo of 188-Hz calibration specimen 
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Table 4.1 Properties of Calibration Specimens 


Mass Polar Moment 
Specimen of Density of Inertia 
Resonant Frequency | Material (b/in?) (ib-ft-sec”) 
od or Tube (J) | Top Plate (J;) 


equipment, a soil column is fixed at the bottom and torque is applied at the top of 
specimen. The torque is applied to the specimen through the drive plate and top platen 
that are fixed to the top of the specimen. As a result, the top of soil specimen is not 
completely free because of the drive plate and top platen attached to it. This attachment 
becomes an added mass that exerts an inertial torque on top of soil specimen. 
Accounting for this inertial torque requires the mass polar moment of inertia of the top 
platen and drive plate. Therefore, it is necessary to carefully measure the mass polar 
moment of inertia of the added mass for the RC test. 

The assembly at the top of the soil column in the RC/TS equipment consists of | 
two different added masses. First the drive plate, which consists of four-armed drive 
plate with magnets and all of the equipment attached to it such as the accelerometer, the 
counter weight, the proximitor target, and screws. And, second the top platen, which is 


attached to the drive plate contacts with the top of soil column. 
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The mass polar moment of inertia of the top platen can be determined directly by 


calculation because it has a simple geometric shape. However, the shape of the drive 
plate is much more complex making a direct calculation for the polar moment of inertia 
impractical. Instead, the mass polar moment of inertia of the drive plate is determined 


experimentally. Method for this experimental determination is described below. 


4.3.1 Drive Plate Calibration Method 

One of the complications in using the RC/TS equipment is the back-electromotive 
force generated in the coil and magnet in the drive system. This is typically called back 
emf effect. The back emf induces a current in opposite direction to the applied current. 
The typical method of testing is using forced vibration, hence the coil circuit is closed, 
allowing the flow of back-emf induced current. This current affects the resonant 
frequency and damping of the system. More details regarding the back emf effect is 
presented in Chapters 5-7. For this study, the drive plate calibration method was 
modified to eliminate the back emf effect. Rather than performing the test using a closed 
drive circuit, an open drive coil circuit was employed. This allows no current to flow 
through the circuit, hence no back emf effect. 

The electrical circuit for the open drive circuit test is shown in Figure 4.5. A fuse 
is connected between the voltage source and the drive system. Instead of driving the coil 
by applying a sinusoidal forcing function, the current is increased linearly resulting in 
increasing of torque on the drive plate. The current is increased until it blows the fuse, 
opening the coil circuit. The torque on the drive plate instantly decreases to zero, and the 


system freely oscillates at its resonant frequency without generating any back emf. 
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Figure 4.5 Electrical circuit for the open circuit drive plate calibration 


The drive plate calibration system can be modeled as a single-degree-of- freedom, 
spring-mass system. The calibration specimen represents the massless spring with 
torsional spring stiffness, Kg, and the top plate of calibration specimen and the drive plate 
represent the mass in the system. Additional analyses were performed accounting for the 
mass of the calibration specimen. These yielded nearly identical to the simpler single- 
degree-of-freedom assumption. Therefore, only the single-degree-of-freedom results are 
presented. 


In torsional mode, the natural frequency, fa can be expressed as: 


-1 | Ko 
" 2m4], +J, 


(4.1) 


where Ko = torsional spring stiffness of the calibration specimen, 


J, = mass polar moment of inertia of the drive plate and the assembly, and 
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J, = mass polar moment of inertial of the top plate of the calibration 


specimen. 


In free vibration, the resonant frequency, fm is: 


f, =f,v1-D’, (4.2) 
where D = damping ratio of calibration specimen. 


Substitution of Equation 4.2 into Equation 4.1 gives: 


K 
pee —t-y1-D’. (4.3) 
2n\VJ,+J, 


The Kg and J, are two unknowns presented in Equation 4.3. Solving for these 
unknowns requires at least one more equation. To obtain additional equations, the test 
was repeated on the same calibration specimen with additional masses of known mass 
polar moment of inertia is fixed on top of the drive plate as shown in Figure 4.6. The 


resonant frequency, and damping, D is then measured and Equation 4.3 is modified to 


become: 
1 K, 7 
=— _ /—__—"_y1-D’, 4.4 
"  In\5,+5, +A vo) 
where AJ = mass polar moment of inertia of an added mass. 


With two equations and two unknowns, the value Ky and J, can be calculated from 
Equations 4.3 and 4.4. 

In this study, three added masses were built for the drive plate calibration. They 
are steel circular disks with different thickness as shown in Figure 4.7. Dimension and 
density of the added masses were carefully determined. The polar moment of inertia of 


the added masses were calculated and tabulated in Table 4.2. The calibration was 
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Added Mass 


Calibration Specimen 


Figure 4.6 Test setup for the drive plate calibration with added mass 


performed using the four calibration specimens described in the previous section. Tests 
without added mass and tests with added masses were performed on each specimen. 


Therefore, there were four equations to solve for Kg and Jo. These equations are shown 


below: 

—— 41- 4.5 
e J, + T Ga 
JED; 4.6 
fnt Sa T +J, T ! Ga 
fa =— — K A ; 4.7 
i Bia ? 67) 
f, s — É _ fp? (4.8) 

= 2J, +J, +A, a l 
where fmo,1,2,3 = measured resonant frequency of the calibration specimen for the 


test with no added mass, mass no.1, 2, 3, respectively, 
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Do,1,2,3 = measured damping ratio of the calibration specimen for the test 
with no added mass, mass no. 1, 2, 3, respectively, and 


AJı 23 = polar moment of inertia of the added mass no. 1, 2, 3, 


respectively. 


02.9750 
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zZ 
9 
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dimensions are in inches 


Figure 4.7 Drawing for added masses for the drive plate calibration 
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Table 4.2 Properties of Added Masses for Drive Plate Calibration 


; Polar Moment of 
ee Material ; Inertia, AJ 
7 (Ib-ft-sec”) 


0.283 2.384E-4 


4.3.2 Drive Plate Calibration Results 

Solutions to Equations 4.5-4.8 give six solutions for values of Kg and J, for each 
calibration specimen. Table 4.3 presents the measured resonant frequency and damping 
of the system from each test. The damping ratio is very small because there is no 
additional damping from the back emf effect. 

Solutions for Ks and J, can be solved graphically by plotting all possible values of 
Jo versus Kg for each test. These plots are shown in Figure 4.8. The intersection point 
represents a solution of two equations and each solution has an average resonant 
frequency for the two equations associated with it. All of the solutions for values of J, 
versus Kg are presented in Table 4.4. It can be seen that solutions based upon tests with 
large differences in resonant frequency result in outlying values of Jo.. This deviation is 
more significant in the stiffest specimen because it has the most difference in resonant 
frequency between two tests. An average value of Kg was calculated from the average of 
the three of the closest resonant frequencies and this value was used to determine 


electromagnetic properties presented in Chapter 6. 
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Table 4.3 Drive Plate Calibration Results 


128-Hz Specimen 188-Hz Specimen 


Added 14-Hz Specimen 61-Hz Specimen 


Mass 
No.! ) 


No 
Added | 14.314 0.020 61.861 128.418 0.010 188.079 0.023 
Mass 


1 13.603 122.718 | 0.008 117.564 | 0.019 


12.988 180.775 0.007 173.917 | 0.017 
3 12.463 113.029 


Note: 1. Top plate is considered as part of specimen 


2.25x10° 2.3x10 
2.20 A 
F 32.2 
£2.15 & 
i 2 
2 = 
52.10 82.1 
2.05 62-Hz specimen 
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1.9 
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a) b) 


2.6x10" 
2.8 
é 5 
324 2 
a C=) 
£ 2 2.6 
£ ES 
[e] m 
ka 
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Figure 4.8 Solutions for Kg and J, for the calibration specimens 
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Table 4.4 All Solutions for Values of Kẹ and J, Obtained from Drive Plate Calibration 


Average Polar Moment ; 
| Resonant of Inertia of Torsiodal Average 
: Stiffness, Ke Kg 
| Frequency, fa | Drive Plate, J, (Ib-fi/rad) (Ib-fi/rad) 
(Hz) (Ib-ft-sec”) 


| 12.725 21393 | 18102 | | 
| fo j , 
| 13.295 


57.972 2.123E-3 345.78 
59.097 2.128E-3 346.55 
60.382 2.149E-3 349.74 


224 15588 


| 128 Hz e | 1560.6 
| 120.730 2.297E-3 1584.2 
122.993 2.317E-3 1597.1 


125.568 2.372E-3 1632.9 


Tane oe AATE 
188 Hz rae ee : : Ase oe 

| 177.938 2.606E-3 3855.2 

: 180.998 2.641E-3 3904.4 

| .73TE-3 


184.427 2.737E 4038.4 


Note: 1. An Average Ke is calculated from three of the closest values represented by 
the values in shading area 
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Figure 4.9 Variation of J, with frequency 


From Table 4.4, the solutions of J, for pair of tests are plotted versus the average 
resonant frequencies in Figure 4.9. Results show a trend of increasing J, with frequency. 
The equation presented in Figure 4.9 represents the average trend of the data. The 
increase of J, is about 28% from 14 Hz to 188 Hz. The increasing of J, with frequency 
found in this study will cause an error of the calculated shear wave velocity, Vs, if it is not 
accounted for. If too low a value of J, is used for calculating V, in a high frequency test, 


the error can be as high as 12%. Therefore, the increasing J, should be taken into 


account. 


Theoretically, Jo of the drive system should be independent of frequency if the 


drive plate moves as a rigid body. This assumption may not be realistic. 
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4.3.3 Finite Element Analysis of the Drive Plate 


Increasing values of J, of the drive plate with increasing frequency may be due to 
non-rigid body motion of the drive plate. As described in Chapter 3, the drive plate 
consists of four aluminum arms plate with a magnet attached to each arm. The non-rigid 
body motion of the drive plate possibly occurs at the connection between the arm and 
each magnet. To verify this assumption, finite element analysis, using a program called 
RISA —2D, was performed to prove this assumption. The drive plate and boundary 
condition were modeled to replicate the experiment setup of the drive plate calibration. 
The stiffness of the connection between the drive plate and magnet was investigated. 

The shape and dimensions of drive plate were modeled with a finite element 
meshes were generated using rectangular plate elements as shown in Figure 4.10. 
Material properties of drive plate and magnet are those of aluminum and steel, 
respectively. The calibration specimen was modeled as four springs. Each spring rotated 
freely in x and y direction and each of them has stiffness of one fourth of the stiffness of 
the calibration specimen. The springs were attached to the drive plate at position of 
screw hole where the drive plate is fixed on top of the calibration specimen. The 
accelerometer and the counter weight were also modeled as a round bar with steel 
material properties. 

The connection between the aluminum plate and magnet was assumed to cause 
the non-rigid body motion. For the RC/TS equipment used at USU, the magnet is 
connected to the drive plate as shown in Figure 4.11. The steel rod is epoxied into the 
magnet and connected to the drive plate with a setscrew. In the finite element analysis, 


the connection was modeled using very thin plate element. The stiffness for this 
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Figure 4.10 Finite element analysis results of the drive plate (exaggerated deflection) 
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connection was adjusted literately in the analysis until the finite element results match 
with the experimental results. The finite element analysis was performed with different 
spring stiffness representing each calibration specimen. The best match value of stiffness 
for this connection was found to be equal to about 220 ksi. This value is much lower 
than the stiffness of the aluminum plate which is about 10‘ ksi. Results of the resonant 


frequency and calculated value of J, are presented in Table 4.5. 


Set Screw 


Magnet 


Steel Rod 


Note: Not To Scale 


Figure 4.11 Connection of the drive plate and magnet 


Table 4.5 Results of the Finite Element Analysis of the Drive Plate 


Polar Moment 
of Inertia of 
Drive Plate, J, 
(Ib-ft-sec”) 


14 Hz 14.104 2.173E-3 
61 Hz 61.374 2.146E-3 
28 Hz 128.971 2.240E-3 


188 Hz 187.891 2.549E-3 


fin for FEM 
(Hz)) 


Specimen 
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In the finite element analysis, a deflection was observed at the connection 
between the magnet and the drive plate. As shown in Figure 4.10b, the magnet tends to 
rotate independently with the drive plate and the rotation angle was observed to be larger 
at higher frequencies. At higher frequencies, the movement of the drive plate and the 
magnets system no longer behave as a rigid body. 

The calculated values of J, from the finite element analysis are plotted along with 
with the experimental results shown in Figure 4.12. The agreement between results from 


both methods verifies that the non-rigid body motion of the drive plate causes increasing 


Jo with frequency. 
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USU Drive Plate 
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Analysis (RISA) 
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Figure 4.12 Comparison of results from the experiment and the finite element analysis 
for the drive plate calibration 
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It should be noted that the results of this analysis apply only to the USU drive 


plate. Other RC/TS equipment may behave differently depending on the type of 


connection between the magnet and drive plate. 


4.3.4 Comparison of Drive Plate Calibration 
of the USU and UT Drive Plate 


To investigate the effect of frequency on the polar moment of inertia for other 
drive plates, the same calibration method was performed using the same calibration 
specimens using drive plate No. 9 of The University of Texas at Austin (UT). The 


comparison between the USU and UT drive plate calibration is shown in Figure 4.13 


3.0x10” 
— USU Drive Plate: J, = 0.0021013+4.673e-10* freq ” 


~~ - UT Drvive Plate No. 9: J, = 0.0021429+1.4982e-8* freq 7 
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Figure 4.13 Comparison of the drive plate calibration for the USU and UT drive plate 
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Only a slight increase in J, with increasing frequency was observed in the UT 


drive plate. The error in calculated shear wave velocity due to the increasing of J, with 
frequency for the UT drive plate would be less than 5%. 

A physical difference between the two drive plates was observed at the 
connection between the drive plate and the magnet. The USU drive plate uses a steel rod 
and a set screw to connect the magnet to the drive plate. Figure 4.14 shows the 
connection of the USU drive plate. The UT drive plate is connected to the magnet using 
an adhesive as shown in Figure 4.15. The difference in connection method explains the 


difference in behavior of these two drive plates. 


Figure 4.14 Photo of drive plate and magnet connection for the USU drive plate 
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Figure 4.15 Photo of drive plate and magnet connection for the UT drive plate No. 9 
4.4 DEFORMATION SENSORS CALIBRATION 


4.4.1 Proximitor Calibration 

The proximitors are mainly used to measure the rotation at low frequency for the 
TS test. The proximitor system for RC/TS equipment consists of two proximitors and 
targets. The proximitor system is a non-contacting measurement. The proximitor probe 
generates a voltage proportional to the distance between the probe and the target. 
Although there is a calibration factor provided by the manufacturer, the calibration is 
required because the calibration factor and linear range vary with the target material. 

In this study, the proximitor calibration was performed using a milling machine to 
obtain an accurate measurement of the gap between the proximitor and target. The 


proximitor was mounted to the machine table and the target was fixed to the mill head. 
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The proximitor output voltage was measured at 5 mils (0.005 in) intervals as the mill 


head was moved away from the table. 

The results of the proximitor calibrations are presented in Figure 4.16. The 
average linear range for the two proximitors is about 0.085 in. Hence, the proximitors 
can measure strain levels up to about 1.5% for soil specimens with a 1:2 (diameter: 


length) aspect ratio. The slopes of the plots for two proximitors and the average value are 


also shown in Figure 4.16. 


Average Linear Range -1.4 to -19.1 V (~0.085 in) 
-20 Center Voltage ~ -10.25 V 

— Proximitor A (Slope = 206.38 V/in) 

- - - Proximitor B (Slope = 207.47 V/in) 

Average Slope = 206.925 V/in 
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an 


= 
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-5 


0 
0.00 0.02 0.04. 0.06 0.08 0.10 0.12 
Gap (in) 


Figure 4.16 Results for the proximitors calibration 
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To confirm the measured proximitor calibration factors, swept sine test was 


performed on the 14 Hz calibration specimen. A constant amplitude torque was applied 
to the specimen over a frequency range of 0 to 155 Hz. The resulting rotation was then 
measured with both the proximitors and the accelerometer. Using the accelerometer 
calibration factor supplied by the manufacturer, a proximitor calibration factor was 
calculated. This calibration factor is plotted versus frequency in Figure 4.17. 

As seen in Figure 4.17, the results were noisy at low frequencies because the 
accelerometer output is very low at low frequencies. At higher frequencies, the plot is 
smoother and the calculated calibration factor is constant. An average value for the 
calculated proximitor calibration factor of 215.02 V/in was obtained from the best quality 
measurement frequency band of 50-120 Hz as shown in Figure 4.17. This value is about 
4% higher than the measured value of 206.925 V/in from the calibration presented 


previously. 
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Figure 4.17 Proximitor calibration factor calculated from rotation using accelerometer 
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The difference of the proximitor calibration factors from both methods is as 


expected, The proximitor calibration on the milling ae is based on the assumption 
that the proximitor probe and target are perfectly parallel. In the actual measurement, the 
specimen is deformed in torsion, therefore, the target is rotated and no longer paralleled 
to the proximitor probe. An improvement for this method is to measure the output 
voltage with the rotation angle. However, this requires an accurate measurement of very 
small rotational angle. In this study, the proximitor calibration factor determined from 
the accelerometer measurement is considered more accurate. Therefore, the proximitor 
calibration factor of 215.02 V/in (0.003322 rad/V) is used. 

4.4.2 Linear Voltage Displacement Transducer 

(LVDT) Calibration 

The LVDT used in this study is Schaevitz Model HR500, which can measure 
vertical displacements up to one inch. The LVDT consists of two parts: the transformer, 
and the core. The core is inserted inside the transformer and the LVDT circuitry outputs 
a voltage proportional to the core position in the transformer. In RC/TS equipment, the 
steel rod is threaded to the core and positioned at the center of the drive plate. 

As in the proximitor calibration, the milling machine was used for the LVDT 
calibration. The core was mounted to the mill head and the housing was mounted to the 
machine table. The output voltage at 5 mils (0.005 in) intervals was measured as the table 
was move away from the mill head. The resulting calibration is plotted in Figure 4.18. 
The line represents the best least square fit to the data. The calibration factor for LVDT 


used for this study is 0.151256 in/V. 
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Figure 4.18 Results for the LVDT calibration 
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CHAPTER 5 


ELECTROMAGNETIC MODEL DEVELOPMENT 


5.1. INTRODUCTION 
This chapter presents the basic concepts of electromagnetism that describe how 
torque is generated in the RC/TS drive system. These same basic concepts also explain 
how back emf effects RC/TS measurements, These back emf effects are quantified. 
Once the back emf problem has been identified, an electromagnetic model is 
developed to mathematically describe the RC/TS drive system. This model can be used 
to describe, predict, and correct for the back emf effect is the RC/TS tests. An approach 


for using this model to correct RC/TS measurements is presented. 


5.2 DRIVE SYSTEM IN THE RC/TS EQUIPMENT 

Torque in the RC/TS drive system is generated from reaction force between coils 
and magnets. As shown in Figure 5.1, each of the four arms of the RC/TS drive plate is 
attached to a magnet. A coil encircles each end of the magnet. The eight coils are in 
series with the coils at either end of magnet wound in opposite directions. When 
energized, the coils generate magnetic fields. These magnetic fields interact with the 
fields from the magnets to generate forces of attraction or repulsion. The torque is 
transmitted from the drive plate through the top platen to the top of soil specimen. The 


magnitude and direction of force is quantified in next section. 
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Section A-A 


Figure 5.1 Drive system in the RC/TS equipment 


5.3 BASIC MAGNETISM 

Wadelsanding the torque created in the RC/TS drive system requires knowledge 
of the concept of magnetism (Halliday and Resnick 1981) including magnetic fields, 
magnetic forces and electromotive forces. These concepts are presented in following 


sections, with specific application to the RC/TS drive system. 
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5.3.1 Magnetic Field 


When a current flows through a wire it creates the magnetic field in the direction 
given by the right hand rule. Ampere’s law is used to determine the strength of this 
magnetic field: 

fBdl = pi, l (5.1) 
where B = magnetic field (tesla = 10° Gauss), 
1= length of the wire (m), 
i = current in the wire(A), and 
Mo = permeability constant = 4rx107 (tesla-m/A). 

Equation 5.1 can be applied to a coil, which is a long wire wound into loops. 
Figure 5.2 shows the lines of the magnetic field in a loosely wound coil. The magnetic 
fields tend to cancel each other between the wires. At points inside of the coil and far 
from the coil, the magnetic field is parallel to the axis of the coil. For the case of coil 
with tightly packed wires, the coil becomes a cylindrical current sheet. The right end, 
from which the magnetic lines emerge, behaves like the north pole of a magnet and left 
end behaves like the south pole of a magnet, as shown in Figure 5.3. Ampere’s law 
applies to a coil is: 


B= Np, (5.2) 
where N is the number of loops. l 


An additional term that needs to be defined is the magnetic flux, Pg. The 
magnetic flux is defined as the number of the magnetic lines that pass through a given 


surface, S and it is expressed as: 
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Coil 


© Current out of Plane 
® Current into Plane 


Figure 5.2 Magnetic field of a loosely wound coil 


. PREPREKE RKRN 


Figure 5.3 Magnetic field of a coil of finite length 
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©, = [B-ds. (5.3) 
The magnetic flux has units of tesla-m* or weber (Wb). If the magnetic field is constant, 
and at right angles to all point on a planar surface of area, A, then Dg becomes: 
®, =BA=Bk, (5.4) 
where 1 = length of the coil loop, and 
x = width of the coil loop. 
A coil carrying currents behaves like a permanent magnet with a north pole at one end 
and a south pole at the other end. A permanent magnet can be modeled as an identically 
shaped coil that has an equivalent surface current flowing around it. Figure 5.4 shows the 
current loops on the cross section of the permanent magnet. In the interior of a magnet, 
the currents are in opposite directions and cancel each other, leaving the net current loop 
at the surface as shown in Figure 5.4. Although these currents are auinaey they satisfy 


Ampere’s law. Therefore, Equations 5.1 — 5.4 are applicable to a permanent magnet. 


VII, 


a) b) 


Figure 5.4 Equivalent surface current in a permanent magnet 
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3.3.2 Magnetic Force 


The basic equation for magnetic force describes the force between two parallel 
wires carrying current as shown in Figure 5.5. The force exerted to the wire a can be 
expressed as: 

F, =i,1xB,. (5.5) 

Equation 5.5 shows that the force vector is the cross product between the current 
vector in wire a, i and magnetic field vector from wire b, By. The magnetic field, By is 
calculated using Equation 5.1. The magnetic field created due to a current on a long 


straight wire at distance, x can be expressed as: 


B, = fet. (5.6) 


Therefore, the magnitude of force in two parallel wires is: 


i „li,i 
Pr a 


(5.7) 


Figure 5.5 Force on two parallel wires carrying currents 
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As explained earlier, a permanent magnet behaves the same as a coil carrying 


current. Therefore, it is possible to describe the force on the magnet due to a nearby coil 
as the force between two coils interacting with each other: one is carrying apparent 
current and another is carrying real current. As shown in Figure 5.6, the magnet has 
equivalent surface current flowing around it and the coil sets up a magnetic field due to 
the real current. Force is applied to the side of the magnet due to the vertical line of the 
magnetic field from the coil with the vector B is at a right angle to the vector i. The cross 
product of the current and magnetic field vectors is the force vector having direction 
shown in Figure 5.6. In addition, vertical forces are also created due to nonvertical lines 


of the magnetic field. However, these forces cancel each other out. 


Figure 5.6 Reaction force between a magnet and coils 
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In RC/TS drive system, the direction of the current is opposite for a pair of coils 


because they are wound in the opposite directions. According to the right hand rule, 
magnetic fields set up by the two coils are in opposing directions. One coil produces a 
repulsive force with the magnet and the other produces an attractive force with the 
magnet as shown in Figure 5.6. The net force moves the magnet in one direction 
depending on the direction of the applied current. Similar results occur for each magnet 
and coils because they are set up in the same north-south pole arrangement all around the 
drive plate. All four magnets are subjected to the same direction of magnetic force thus 
they move in phase with each other creating torque on the drive plate as shown in Figure 
5.7. 

The equivalent surface current of the magnet was only used to verify the direction 
of the torque, it is not a real current that can be measured in the laboratory. Rather than 
calculating the force acting on the magnet due to the magnetic field of the coil, in this the 
force was found by calculating the force acting on the coil due to the magnetic field of 
the magnet, which is measurable. Considering only the lines of magnetic field at right 
angle to the direction of current flowing through the coil, the magnitude of the force is: 

F=iNI1,B, (5.8) 
where i= current in the coil, 
l = circumference length of coil loop, 
N = total number of loops in the coil in x and y direction (N = NxNy) 
illustrated in Figure 5.9, and 


B = magnetic field of the permanent magnet. 
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Figure 5.7 Torque generated on the drive plate 


From Equation 5.8, it can be seen that the force created in the RC/TS drive system 
is proportional to the current flowing through coils. Torque created in the system can be 


determined by introducing the moment arm of the drive plate, d, and modifying Equation 


5.8 to: 
T=(N,N,)I,Bdi=kji, (5.9) 
where k; = torque-current factor. And, 
k; =N,N,1,Bd. (5.10) 


Reproduced with permission of the copyright owner. Further reproduction prohibited without permission. 


97 


The factor k; is an important factor used for calculating torque generated in the 
torsional shear test. As shown in Equation 5.10, if the circumference length of the coil 
loop and the moment arm of the drive plate are constant, the value ki is proportional to 
the eae of loops in the coil and strength of the magnetic field of the magnet. This 
implies that torque in RC/TS system can be easily increases if either the numbers of turn 
increase or a stronger magnet is used. The second option is practical. However, there is 
also a disadvantage to using a stronger magnet to create more torque. The disadvantage 
is discussed in next section. Calculations of k; value using Equation 5.10 is also 


presented in Chapter 6. 


5.3.3 Electromotive Force 

When the magnet moves axially inside the coil, the movement of the magnet also 
induces an electromotive force (emf) according to Faraday’s law. The induced 
electromotive force in a circuit is equal to the rate at which the flux through the circuit is 
changing: 


dd 
Sat =a (5.11) 


The emf is a potential difference between two points in a circuit therefore it has a unit of 
joule/coulomb or volts. 

The direction of the induced emf is specified by Lenz’s law using the 
conservation of energy principle. Lenz’s law states that: the direction of the induced 
current will appear in such a direction that it opposes the change that produces it. Notice 
that Faraday’s law also suggests this opposition by using minus sign as presented in 


Equation 5.11. Lenz’s law refers to induced current (iems), which means that it applies to 
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closed circuits such as that shown in Figure 5.8. Ifthe north pole of the magnet is pushed 


toward the coil, an induced current sets up a magnetic field that opposes to the motion of 
the magnet. As shown in Figure 5.8, the side of the coil which the north pole of the 
magnet is moved toward becomes a north pole, and the two north poles will repel each 
other. The right hand-rule shows that the for the magnetic field setup by the coil to 
emerge to the right side of the coil, the induced current must be in the direction in Figure 


5.8. In RC/TS equipment, this phenomenon is called “back emf.” 


Enf R 


Figure 5.8 Electromotive force induced by the movement of magnet 
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The induced current due to the back emf always counteracts with the applied 


current in the coil. The back emf is a source of energy loss in the drive circuit. This 
energy loss acts as damping generated by the RC/TS equipment, therefore is it is called 
equipment generated damping. The equipment generated damping is combined with the 
material damping of the soil during RC/TS measurements. An error in soil material | 
damping due to the back emf effect has = one of the weaknesses in the RC/TS test. 

To quantify the back emf effect in the RC/TS drive system, the cross sections of a 
coil in the x, y, and z planes is illustrated in Figure 5.9. Assuming that the magnetic lines 
created by the magnet passes through the area of the coil loop, Ac, the flux Mg enclosed 
by the loop for every turns of the coil is equal to: 


D, =BA,. (5.12) 


x 


a) x-z section b) y-z section 


Figure 5.9 Cross sections of coil loop 
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Applying Faraday’s law to the coil, an emf appears for every loop of the coil. 


With the assumptions that the flux through every loop of the coil is added and that the 
coil is an ideal coil, the induced emf is given by: 


_d(N®,) _ d(N,N,BA,) 


€ 
a dt dt 


(5.13) 


For the case of coil and magnet in RC/TS drive system, the magnetic flux changes 
as the magnet moves. This change is the increase or decrease in total Dg through the coil 
loop caused by moving the magnet toward or away from the coil. The total Dg is 
represented by the terms N® , or NxNy®g in Equation 5.13. Since the magnet is only 
moved in y direction, the change in Dg, therefore, is only caused by a change of Ny. 
Equation 5.13 can be rearranged as: 


dN 
Emt = BAN, >. (5.14) 


The change of Ny or dN, is directly related to a moving distance dy as shown in 


the relationship below: 


y 


N 
dN, =— dy, (5.15) 
y 


where y = total width of the coil loop in y-direction as shown in Figure 5.8. 
Substitution of Equation 5.13 into Equation 5.12 gives: 


_BA,N,N, dy _ BA,N,N, 


€ = 
emf y d t y 


v, (5.16) 


where v = velocity of the magnet. 
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Equation 5.16 implies that emf is proportional to the velocity of the magnet. In 


the RC/TS system, it is more useful to express the emf in relation to the rotational 
velocity of the magnet. Equation 5.16 then becomes: 


BA, N,N 
Ee _2 c> x ae A (5.17) 
d y dt dt 


where A = rotational velocity of the magnet, and 


kg = back emf-rotational velocity factor which is equal to, 


2 BA.N.N, 
ks 
d sy 


(5.18) 

The factor kg derived in Equation 5.18 is an electromagnetic property of the 
RC/TS drive system. This factor represents the amount of back emf created in the drive 
system, hence higher value of kg means higher effect of back emf. From Equation 5.18, 
the factor kp is proportional to the number of loops in the coil and the magnetic field of 
the magnet which are the same as factors affecting the value of ki. More back emf is 
expected if numbers of loops in the coil are increased, or a stronger magnet is used. As 
presented in the previous section, increasing the numbers of loops in the coil and using a 
stronger magnet increases capacity of the RC/TS drive system in generating torque. 
Consequently, more back emf effect is also generated and resulting in higher equipment 
generated damping. However, if the back emf effect can be predicted and the equipment 
generated damping is taken into account precisely, the problem with the back emf is then 
solved. In this study an electromagnetic model of the RC/TS drive system is developed 
to describe the back emf effect that can be used to predict the equipment generated 


damping. 
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5.4 | ELECTROMAGNETIC MODEL 


An electromagnetic model was developed using the combined electrical, magnetic 
and mechanical properties of the RC/TS apparatus. The purpose of the model is to 
describe the behavior of the system including the effect of the back emf. The model is 
used to predict the equipment generated damping due to the back emf allowing a precise 
measurement for the material damping of soil. This section presents the step-by-step 
development of the electromagnetic model. The effect of the back emf in the model is 
described. In addition, the effects of back emf effect in the RC and the TS tests are 
considered separately. Methods are presented to account for the back emf effect in each 
test. 

Figure 5.10 shows the RC/TS drive circuit. A voltage source with internal 
resistance of R, generates voltage potential, Vs. The resistor Rs and inductor L, represent 
the total resistance and total inductance of the eight drive coils, respectively. The back 
emf (Eeme) appears in the circuit as an additional voltage source. The circuit is connected 
to ground at point b. Current, i, is assumed to flow through the circuit in the clockwise 
direction, as shown in Figure 5.10. At point a, the voltage potential V, is given by: 

V, = V, —iR,. (5.19) 

Applying Kirchhoff s voltage law to the drive circuit, starting at the voltage 


source and going around the circuit in the clockwise direction results in: 


Vie ot ey =0, (5.20) 
dt 
where L.(di/dt) = drop in voltage potential due to the coil inductance. 
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Voltage | ———s hn YE ee 
Source Drive Coils 


Figure 5.10 Drive circuit of the RC/TS equipment 


Substitution of Equation 5.18 and Equation 5.19 into from Equation 5.20 and gives: 
V, -(R, +sLc) =s0k,, (5.21) 
where s = Laplace variable (s = jæ), 
œ = circular velocity. 
Solving for the current in the coil yields: 


_ V, —s0k, 


= ; 5.22 
R, +sL, DRR 
As derived in previous section, the torque is proportional to the current. An 


expression for torque can be obtained by substituting Equation 5.22 into Equation 5.9: 


T=ki=k V, ~S0ky 


: ; 5.23 
! ' R,+sL, G2) 
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Soil in the RC/TS system can be described as a viscoelastic material represented 


as a single degree of freedom system. In this type of system, the elastic property is 
represented by a spring with stiffness, kọ and the damping is represented by a dashpot 
with damping coefficient, cp. The simplest and most widely used model to predict the 
viscoelastic behavior is called the Kevin model. The drive system of the RC/TS 
equipment is represented by a mass with mass polar moment of inertia, J attached to the 
model. Torque is applied to the system as shown in Figure 5.11. 

From the assumed single degree of freedom system, torque is evaluated from the 
following differential equation: 

T =k,0+c¢,s0+Js0 i (5.24) 
where ke = torsional spring stiffness, 
Ce = torisonal viscous damping coefficient, and 


J = polar moment of inertia of the system. 


Figure 5.11 Kevin model 
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Substitution of the torque obtained from the drive circuit (Equation 5.23) gives: 


x,| V2 82 Ke | 4 940,90 + 3870. (5.25) 
R,+sL, 


Division by 0 and rearrangement gives: 


Vs ee ee oe ae (5.26) 
0 (| R,+sL, R,+sL 


c c 


The complex terms ann in Equation 5.26 can be expanded and rearranged into a 
+S 


c c 


form a+sb as follows: 


kk, (i =sL, e k:k,R, —sk;k,L, ; (5.27) 


R, +sL, | R,-sL, R? +o} 


By substituting Equation 5.27 into Equation 5.26, and j@ for s, Equation 5.26 becomes: 


V. R,+joL,[(,  okk,L,) . kk,R 
a = ene || k + Be j jag) cg + 2 |- 0) J. 5.28 
(s mpe mariz) | m 


Therefore, the voltage transfer function, Hy for the rotation and applied voltage is 


expressed as: 


k; 
H, =S= R tjok. = (5.29) 
V, ik, +k,, )+ jac, +65 -a 7| 
where Keg = equipment spring stiffness, 


Ceq = equipment viscous damping coefficient, and 


œ = angular velocity (@ = 2zf). 
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The terms Keg and Ceq are: 


o’k,k,L, 
ag PETE ; (5.30) 
KKR, (5.31) 


C= 
“TR toL 
- The current transfer function for rotation and current, Hj can be expressed as: 


0 k, 


H. =—= ; 5.32 
i Ik, + jac, -o| enn) 


Equation 5.29 is the transfer function developed to predict the response curve for 
soil including the back emf effect in RC testing using a voltage measurement. The 
voltage measurement for RC test is commonly used because the amplitude of sinusoidal 
signal generated by the source generator is quantified in terms of voltage. Equation 5.29 
shows that an increase of the stiffness and viscosity of the system due to the back emf are 
expected as presented by the terms keg and Ceg, respectively. The term keq and Ceq are 
function of frequency and four parameters: coil resistance, Re, coil inductance, Le, torque- 
current factor, k;, and back emf-rotational velocity, kg. The value of keq tends to 
increases as the frequency increases while the value of Ceg tends to decrease as the 
frequency increases. The frequency effect on value of keg and Ceq is investigated 
experimentally in Chapter 7. 

Equation 5.32 shows that there is no back emf effect in RC/TS testing wien using 
a current rather than a voltage source. The back emf effect can also be eliminated by 


measuring current in RC/TS tests using voltage sources. 
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5.4.1 Behavior of Back Emf in RC/TS Drive System 


The single degree of freedom system composed of a spring-dashpot system called 
the Kevin model can be used to represent the RC/TS system. The same model can be 
used to describe the effect of back emf. 

In Equation 5.29, the additional terms keg, and Ceq were introduced to the transfer 
function for voltage measurement due to the back emf effect. The back emf is therefore 
San be described as an additional spring and an additional dashpot connected in parallel 
with the original spring and dashpot representing a soil specimen in the system as shown 
in Figure 5.12. The spring-dashpot system with back emf effect can be described with 


the following equation: 


T =(k, +k,,)0+(c, + Coq )SO +Js"0. (5.33) 


Figure 5.12 Modified Kevin model including the back emf effect 
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In Figure 5.12, the kọ and cg in parallel with keq and Ceq implies that an applied 


torque generated in the RC/TS drive system is not all transferred to kp and cg because it is 
partially transferred to keq and Ceg. Consequently, the measured kg and cg of the system 
does not represent the stiffness and damping of the soil behavior but a combination of the 
soil and the back emf. As a result, the soil stiffness and damping are overestimated if the 
back emf is not taken into account. Methods used to account for the back emf effect in 


RC and TS tests are presented in next section. 


5.4.2 Back Emf Effect in RC/TS Tests 

As described in Chapter 3, the RC and TS tests are performed at different loading 
frequencies. The RC test is performed at first mode resonance, which is typically higher 
than 10 Hz while the TS test is performed at very low frequencies usually in the range of 
0.01-1Hz. Therefore, the effects of back emf in the RC and TS tests are taken into 
account differently. 

In the RC test, Equation 5.29 can be applied directly to the response curve 
measured in the RC test using voltage source. To account for the back emf effect, the 
terms keg and Ceq can be calculated if the model parameters, Re, Le, ki, kg are known. The 
effect of back emf in resonant frequency and damping can be predicted using the 


following equations: 
Ape AL Ky + Keg h-2D0 +D.) SERET and (5.34) 
aSa N g Be Ne ae 


paa (5.35) 
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where D = damping of system, 


ce = critical damping coefficient (c= 4af, LJ), 


fa = undamped natural frequency of specimen (f, = ——=—— ), and 


XJ = total polar moment of inertia of the system including such as drive 
plate (Jo), top plate (J;). 
The Afm and Deg are the effects of back emf on resonant frequency and damping 
due to the equipment spring stiffness, keg, and the equipment viscous damping 
coefficient, Ceg described by the model. The Afm is an increasing of resonant frequency 
and the Deg is an equipment generated damping. 
In TS test, soil behavior is analyzed based on characteristics of the hysteresis 
loop. The hysteresis loop is constructed from a measured torque and rotation. Torque 
applied to the soil specimen is partially transmitted to the keg and Ceq as shown in Figure 
5.12. The torque can be corrected by rearranging of Equation 5.31 as shown below: 
T, =T,, - (k.,9 + ¢,,80) = k,0+¢,80 + Js’6 , (5.36) 
where Te = corrected torque, and 
Tm = measured torque. 
Equations 5.34-5.35 and 5.36 can be used to account for the back emf effect for 
RC and TS test, respectively, if the terms keq and Ceg are known. In order to determine the 
Keq and Ceq, the model parameters: Re, Le, Ki and kg must be evaluated. These parameters 
are the electromagnetic properties of the RC/TS drive system that can be measured 
experimentally using the measurement method developed in this study. Details on the 


measurement method used for each parameter and results are presented in Chapter 6. 
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CHAPTER 6 


CALIBRATIONS OF DRIVE PLATE ELECTROMAGNETIC PROPERTIES 


6.1 INTRODUCTION 

This chapter presents measurements of the electromagnetic properties of the 
RC/TS drive system. These properties are the parameters of the electromagnetic model 
developed in this study. The testing methods and analysis procedures used to evaluate 
each parameter are described. All measurements were performed on both the USU and 
the UT drive plates. A comparison of the USU and UT drive plates are included at the 
end of each section. The experimental procedures developed in this study are simple and 
use readily available equipment. Where possible, measurements were performed over a 
range of frequencies and amplitudes to evaluate frequency effects and nonlinearity. 

Additionally, the magnetic field strength of the drive plate magnet was measured. 
These magnetic field were then used to determine the theoretical values of the 
electromagnetic properties measured in this chapter. The theoretical value of the 
magnetic properties and the measured value were compared and the ratio was used to 


calculated efficiency of the drive plate. 


6.2 MODEL PARATMETERS 

As derived in Chapter 5, the proposed electromagnetic model consists of four 
parameters, coil resistance, Re, coil inductance, Le, torque-current factor, kj, and back 
emf-rotational velocity factor, kp. In order to verify and apply the model, these 
parameters need to be determined. An independent measurement is used to measure the 


model parameters. The measured model parameters are then used to predict the 
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equipment generated stiffness, Keq and equipment generated viscosity, Cog. These 


independent measurements are presented in the following sections. 


6.2.1 Coil Resistance and Inductance 

The coil resistance, Re, and inductance, Le, can be measured directly using an 
ohmmeter and an inductance meter. The ohmmeter measures the resistance under D.C. 
conditions, while the inductance meter measures the inductance at a single high 
frequency (typically 5-10 kHz). This study focused on measuring R, and L, over the 
range of frequencies typically encountered in RC/TS testing. The measurements were 
performed using a stiff aluminum specimen so that the drive plate rotation was negligible. 
This eliminated back emf effects. Figure 6.1 shows a drawing and a photo of the stiff 
specimen developed for the measurement of R, and Le. 

The drive system of the RC/TS drive system was modified in this study so that R. 
and L, could be measured in the same test. The modified circuit is shown in Figure 6.2. 
A reference resistor, Rien Was connected in series with the source and the coils. The 
measurement was performed using swept sine test or applying a sinusoidal voltage signal 
over a frequency range of 0 to 400 Hz to the drive coil circuit. Measurements were made 


of voltage, Vı and V2 , as shown in Figure 6.2. The current flowing through the circuit is: 


pee, (6.1) 
Rut 
where V1-V2 = voltage potential drop across the Rye; 


V1 = voltage output from the source, and 


V2 = input voltage to the coils. 
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60.1875 


3.7500 


5.5000 


3.7500 


dimensions are in inches 


b) 


Figure 6.1 Drawing and photo of stiff aluminum specimen 
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Vi Reef V2 
| Re | 
r Rs 
CM lg, i 
Voltage} | & fom = 0 | 
Source 


Drive Coils 


Figure 6.2 Electrical drive circuit for R. and L., measurements 


All the coils are connected to ground, therefore V2 can be obtained from: 
V, =i(R, +joL,). (6.2) 
Substitution of the current, i, and expansion gives: 


V-V: 


V, = R, + joL,). (6.3) 


ref 
In this equation, the coil resistance, Re is in the real part of the equation and the coil 
inductance, Le is in the imaginary part of the equation. As the value of V; and V2 are 


measured a any given Res the coil resistance, R, and inductance, L, can be obtained 


from: 
R,=re V: Re |» and (6.4) 
V, g V, 
1 V. 
L =—i 2 Rol. 6.5 
c O imag V, w V, a ( ) 


Reproduced with permission of the copyright owner. Further reproduction prohibited without permission. 


114 
Units of R, and L, are ohm (Q) and henry (H), respectively. The tests were 


performed with reference resistor, Rre of 25, 50, 100 ohms. Ref was varied to evaluate 
any nonlinearities with current. The linearity of the current was also evaluated by using 
the source voltage of 0.01, 0.1, 0.5, 2 volts. 

Measurements of coil resistance and coil inductance for the USU drive plate are 
shown in Figures 6.3-6.5. Results for the coil resistance show a slight increase with 
frequency, about 3% over the measured frequency band. Measurements using different 
reference resistors and source voltages show no significant nonlinearity in the coil 
resistance. 

Inductance measurements for the USU drive plate are also shown in Figures 6.3- 
6.5. There is no variation of coil inductance with frequency. There is a slight source 
voltage effect with the maximum effect observed in the tests at a voltage of 2 volts with 
Ref = 25 and 50 ohms. The outlying results have slightly higher inductance values, 
however, similar results were not observed at the test of 2 volts with Riet = 100 ohms. 
There is no adequate explanation of this effect, but the variations are very small and can 
be neglected. 

The resistance and inductance of the coils was estimated by averaging of all the 
results. The average value of coil resistance, Rc, is 45.6 ohms. The average value of coil 
inductance, Le, is 0.0256 henries. These values are used for the predicting the equipment 
generated damping using the electromagnetic model. 

Rotation of the drive plate was also measured during the coil resistance and 
inductance measurements. The siesta and proximitors were used to measure this 


rotation. The results are shown in Figure 6.6. Generally, accelerometers do not perform 


Reproduced with permission of the copyright owner. Further reproduction prohibited without permission. 


115 


9 
2 
Source Level (Rref = 26 ohms) 
—— 10 mV (Avg = 45.599 ohms) 
—— 100 mV (Avg = 46.637 ohms) 
—— 500 mV (Avg = 45.643 ohms) 
—— 2 V (Avg = 45.786 ohms) 
0 100 200 300 400 
30x10" 
S 
Q 
=] 


Source Level (Rref = 26 ohms) 

10 mV (Avg = 0.02562 H) 
—— 100 mV (Avg = 0.02563 H) 
—— 500 mV (Avg = 0.02565 H) 
—— 2 V (Avg = 0.02591 H) 


0 100 200 300 400 


Frequency (Hz) 


Figure 6.3 Coil resistance and inductance for Ref of 26 ohms at different source voltage 
levels for the USU drive plate 
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Figure 6.4 Coil resistance and inductance for Ryer of 50 ohms at different source voltage 
levels for the USU drive plate 
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Figure 6.5 Coil resistance and inductance Rrer of 100 ohms at different source voltage 
levels for the USU drive plate 
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well at very low frequencies, thus the noisy accelerometer signals at low frequency 
inthese tests are expected. Unlike the accelerometer, the proximitors only perform well 
at frequencies below 50 Hz. From these results, it is confirmed that the rotation is 
negligible during the measurements, with all roation being less than 4 x 10° radians. 

The R. and L, measurements were also performed on the UT drive plate No. 9 at 
The University of Texas at Austin. Results are shown in Figures 6.7 and 6.8 for Ref of 
24 and 100 ohms. The variation of R, and L, with frequency is very similar to the results 
observed for the USU drive plate. The overall average value of R, is 44.7 ohms and L, is 


0.02375 H. 


Accelerometer Measurement 
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Figure 6.6 Rotations measured during R, and L. measurement on the USU drive plate 
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Figure 6.7 Coil resistance and inductance for Ref of 24 ohms at different source voltage 
levels for the UT drive plate No. 9 
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Figure 6.8 Coil resistance and inductance for Ref of 100 ohms at different source 
voltage levels for the UT drive plate No. 9 
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6.2.2 Torque-current Factor 

The torque-current factor (kj) is used for cainllsane torque in RC/TS tests. The 
factor is also one of the parameters in the electromagnetic model used to predict the 
equipment generated damping from the back emf effect. As explained in the previous 
chapter, torque generated in the RC/TS drive system is proportional to the applied 


current, thus the factor k; can be expressed as: 
> (6.6) 


where k; = torque-current factor (Ib-ft/A), 
T = torque generated to the drive plate, and 
i= current in the drive coil circuit. 

Determination of the factor, kj requires the measurement of torque and current. In 
this study current was measured directly using a current proportional signal from the 
power amplifier. Torque was determined using two different methods. First, by 
performing the torsional shear tests using calibration specimens with known torsional 
stiffnesses. Second, by measuring torque directly using a torque sensor. Details of these 
two methods are presented below. 

In the first method, torsional shear tests were performed on each calibration 
specimen, each having different stiffness. The rotational stiffness, Ke, of each specimen 
was calculated along with the polar moment of inertia of the drive plate presented in 
Chapter 4. Average values of Kg for each specimen are used for this measurement and 


are presented in Table 4.4. 
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Torsional shear tests were performed on the calibration specimens at frequencies 


of 0.01, 0.05, 0.1, 0.5, and 1 Hz. Current and rotation were measured during the tests. In 
addition, the applied voltage was also measured to determine the torque-voltage factor, 
ky. Torque was calculated using: 
T=K,0, (6.7) 
where Kg = rotational stiffness of calibrated specimen (Ib-ft/rad), and 
8 = rotation. 

Plots of torque versus current and torque versus voltage are shown in Figure 6.9. 
Because material damping of calibration specimens is very small, the hysteresis behavior 
is negligible. A linear relationship between torque and current, and torque and voltage 
was observed. Factors kj and k, are determined from slope of the plots in Figure 6.9. 

The k; and k, values obtained from the torsional shear test on each calibrated specimen at 
different frequencies are tabulated in Table 6.1. 

Results for each specimen show consistency in the values of kj and k, obtained 
from the tests at different frequency. Slight differences exist in the results between the 
specimens. As explained earlier, the torque calculated in this test is based on the 
torsional stiffness of the specimen, which was estimated during the drive plate 
calibration. Errors in the value of Kg will result in errors in calculations of torque. 


Despite this, the variation of k; and ky is very small. The overall average results of this 


test are kj = 0.3172 Ib-ft/A, and ky = 0.007075 Ib-ft/V. 
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Figure 6.9 Torsional shear test results on 62-Hz calibration specimen 


Table 6.1 Factors ki and ky Determined from Torsional Shear Tests on Calibration 
Specimens 


14-Hz Specimen 62-Hz Specimen 128-Hz Specimen 188-Hz Specimen 
(Hz) 


0.30480 | 0.0068502 | 0.32037 | 0.0071377 | 0.32474 | 0.0072577 | 0.31993 | 0.0070971 
0.30476 | 0.0068519 | 0.31989 | 0.0071278 | 0.32210 | 0.0071593 | 0.32121 | 0.0071249 
0.30478 | 0.0068517 | 0.32006 | 0.0071404 | 0.32168 | 0.0071565 | 0.32043 | 0.0071105 
0.30507 | 0.0068595 | 0.31997 | 0.0071467 | 0.32112 | 0.0071576 | 0.3213 | 0.0071344 
1,00 0.30623 | 0.0068676 | 0.31993 | 0.0071460 | 0.32318 | 0.0072039 | 0.32106 | 0.0071241 


0.30513 | 0.0068546 | 0.32037 | 0.0071397 | 0.32256 | 0.0071870 | 0.32079 | 0.0071182 
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The second method to determine k; value is to measure torque directly using a 


torque sensor. The torque sensor used in this study was a Futek model T5100. It has 
capacity to measure torque up to 160 in-oz (0.83 ft-lb). The calibration factor of the 
torque sensor is 462.8 Ib-ft/V per excitation voltage. When performing this 
measurement, the torque sensor was mounted rigidly to the top of the aluminum column 
designed and custom built to fit the torque sensor and drive plate of RC/TS aatapinant 
The calibration factor of the torque sensor provided by manufacturer was used to 
determine torque. The factor k; was then obtained from the ratio of torque and applied 
current. 

The variation of kj versus frequency was not measurable using the torque sensor 
because the resonant frequency of the combined torque sensor-aluminum column-drive 
plate system was found at a frequency at about 30 Hz. Thus, the torque measurement 
with frequency encounters significant dynamic effects. As a result, the test was only 
performed under static conditions. Figure 6.10 shows the k; value at different source 
levels plotted versus time. The average value of kj from these tests is 0.31795 (Ib-ft/A). 

The values of factor ki from both methods are very close. However, the 
frequency effect on the factor ki is uncertain. As explained in Chapter 5, the factor k; is 
related to the magnetic field of a magnet. The possibility of magnetic field varying with 
frequency is discussed in next section. Despite this, the variation of k; with frequency 
that could theoretically occur does not cause a significant effect on the predicting 
equipment generated damping. This is discussed in Chapter 7. For this study, a value ki 
equal to 0.318 is used for calculating torque and as a model parameter, which is assumed 


to be constant at all frequencies. 
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Figure 6.10 Static test results for k; measurements of the USU drive plate 


The torque sensor was used to measure the factor k; for the UT drive plate. 
Results are plotted in Figure 6.11. An average k; value of 0.657 Ib-ft/A was measured. 
The k; of UT drive plate is about twice the value measured for the USU drive plate. The 
difference in kj values results from the different types of magnets in the drive plate. 
Referring to Equation 5.10, the value of kj is proportional to the numbers of loops in the 
coil and the strength of magnetic field of the magnet. The magnets used in the USU drive 
plate are made from Alnico 5 and the magnets used in the UT drive plate No. 9 are 
maded from Alnico 8. The Alnico 8 magnet generates higher magnetic field than the 


Alnico 5 hence, the higher value of k; for the UT drive plate. 
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Figure 6.11 Static test results for k; measurements of the UT drive plate No. 9 


6.2.3 Back Emf-— rotational Velocity Factor 
As described in Chapter 5, the back electromotive force (back emf), Eem is 


created from the movement of the magnet through coil. The back emf-rotational velocity 


factor, kg is obtained from rearranging Equation 5.13 to gives: 


€ 
k, ==, (6.8) 
o8 
Me 


where = rotational velocity. 


To measure the value of kp, the drive coils were rewired and separated into two 


groups. Each group has two pairs of coils located opposite each other. One group was 
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used as the drive coils and the other was used as the response coils. Configuration of the 


drive coils and response coils and magnet for kg measurement are shown in Figure 6.12. 


Drive Coils 


4 a 
T oO 
ő 3 
2 7 Function 
2 A Generator 
© 
Q = 
7 a 


Drive Coils 


Figure 6.12 Configuration of coils and magnets for kg measurement 
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The tests were performed using the calibration specimens. Electrical circuits of 


the drive and response coils are presented in Figure 6.13. Sinusoidal voltage with 
constant amplitude was applied to the drive coils over a frequency range of 0-400 Hz, 
creating the movement of the drive magnets and the magnets in the response coils. The 
back emf voltage occurring from movement of the magnets in the response coils and the 
drive plate rotation was measured. The value of kg that represents all four coils is simply 


calculated from: 


V V 
k, =2 26 / =2 Ta i (6.9) 
(0) 
Ya 4 
where Vout = voltage potential across the response coils. 
Vout 
05R, } 
i R ' 
Vs 0.5L, | 
Voltage Eemf ' 
Source Pae [ecscinsd 
Drive Coils 

a) Drive coils b) Response coils 


Figure 6.13 Electrical circuit setup for ks measurement 
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In this study, an effect of crosstalk voltage was discovered during kg 


measurement. The crosstalk effect occurs due to the voltage coupling between the 
response coils and the drive coils. The crosstalk tends to interfere with the measured Vout 
in the response coils. When the crosstalk is combined with the measured Vout, the 
calculated value of kg becomes unrealistic. As shown in Figure 6.14, the crosstalk effect 
was investigated by an evaluation of the phase shift between the measured rotation and 
the Vou. According to Faraday’s law, the back emf is proportional to the velocity or in 
this case, the rotation. This means that there should be no phase different between the 
back emf and rotation velocity. The phase shift of 180 degree observed in Figure 6.14a is 
clearly the effect of the crosstalk. If the crosstalk is not taken into account in this 
measurement, the value of kg would be as shown in Figure 6.14b, with negative value of 
kg at frequencies below the phase shift occurs. The physical meaning of the negative kg 
is that the back emf turns the magnet and coil into a self-perpetuated motor driving itself, 
a physically impossible condition. 

The crosstalk effect has to be taken into account for the kg measurement. The 
simplest approach is to measured the crosstalk voltage and subtract it out from the 
measured Vou. The crosstalk voltage needs to be measured independently without the 
back emf effect. To accomplish this, the crosstalk measurement was performed using the 
rigid aluminum specimen. As explained earlier, the rigid specimen was used so that the 
back emf effect can be neglected. The cross talk test was performed using the same 
circuit as the kg measurement, and Vout was measured and considered as the crosstalk 


voltage, Ve, as shown in Figure 6.15. 
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Figure 6.14 Crosstalk effect on kg measurement 
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Figure 6.15 Variation of crosstalk voltage with frequency 


The results shown in Figure 6.15 show that the crosstalk voltage developed during 
the measurement. The crosstalk increases up to about 4 mV as the frequency increases to 
400 Hz. To account for the effect of crosstalk voltage, the behavior of crosstalk was 
investigated. Results from crosstalk measurement showed that the crosstalk was 
observed to be a function of the rate of change of current, ay which can be expressed 
as: 

di 
V,=k,—=k,ijo. 6.10 

¢ ot oJ ( ) 


c 


The constant kÀ is called the crosstalk factor. The value of k, is a complex number 


that consists of two parts; real and imaginary. The plot of k, value versus frequency is 
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shown in Figure 6.16. There is a slightly variation in the real part of k, with frequency 


and the imaginary part ks is very small, thus it was neglected. Therefore, a ke value of 


—3.6x10° (V/A/s) is used for this study. 
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Figure 6.16 Variation of the crosstalk factor with frequency 
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The effect of crosstalk is taken into account by applying the crosstalk factor into the 


measured current when performing the kg measurement. As a result, Equation 6.9 


becomes: 


v 7 
k, =2 Vu tk jO (6.11) 
jo 


It is noted that the value of kg calculated from Equation 6.11 is also a complex 
number. As explained earlier, the back emf is proportional to the velocity of the magnet. 
After correcting the crosstalk voltage, the back emf voltage should be in phase with the 
rotational velocity hence the real part kg is a constant and the imaginary part is equal to 
Zero. 

Results of kg measurement of the USU drive plate showing plots of the kg value 
varying with frequency for four calibrated specimen are presented in Figure 6.17. It is 
observed the overall value of kg is much smaller after being corrected for the crosstalk 
voltage and no negative values of ks were measured. There is an increase of the real part 
kp with frequency especially at frequency above 200 Hz. The same variation but less 
significant is also observed for the imaginary part kg. The value of imaginary part kp is 
very small and is neglected. Test results from different specimens are identical except for 
the outlying results from the specimen A with the resonant frequency of 14 Hz. This 
specimen is very soft and rocking motion may affect the test result. Curve fitting was 
performed for the average, upper and lower bound for the value of kg varying with 


frequency, the curve fitting equations are presented in Figure 6.17. 
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Figure 6.17 Variation of kg with frequency for the USU drive plate 


Reproduced with permission of the copyright owner. Further reproduction prohibited without permission. 


135 
The same kg measurement method was performed to measure the factor, kg for 


the UT drive plate. Results for real part and imaginary part of kg are plotted with 
frequency in Figure 6.18. The increasing trend was also observed for the UT drive plate. 
There is a less significant effect of the calibrated specimen on the kg measurement, 
therefore, the curve fitting was performed on the average real part kg. The curve fitting 
equation is presented in Figure 6.18. To compare this result with the result for the USU 
drive plate the curve fitting functions for both drive plates are plotted as shown in Figure 
6.19. The factor kg of the UT drive plate are found to be almost twice times of the factor 
kg for the USU drive plate at low frequencies. The increasing rate of kp values for the 
UT drive plate was lower than for the USU drive plate. 

Referring to Equation 5.18, the value of kg is a function of numbers of loops in 
the coils and magnetic field of the magnet. As discussed in previous section, different 
types of magnet were used for the UT drive plate and the USU drive plates. The stronger 
magnet is used for the UT drive plate. Therefore, a higher value of kg for the UT drive 
plate is expected. 

There is no clear explanation for the variation of kg versus frequency. One 
possibility is that the magnetic field may vary with frequency. The change of magnetic 
field with frequency is not yet fully understood and requires further study. For the USU 
drive plate, the curve fitting equation for the average kg is used and the difference. 
between values of the upper and lower bound kg is considered as an uncertainty. The 
uncertainty in kg value of the USU drive plate is evaluated and discussed in Chapter 7. 
For the UT drive plate, the average curve fitting equation is used to represent the value of 


kg with frequency. 
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Figure 6.18 Variation of kg with frequency for UT drive plate No. 9 
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Figure 6.19 Comparison of kg values for the USU drive plate and the UT drive plate 


6.3 DRIVE PLATE EFFICIENCY 

In this section, the strength of magnetic field, B, of the magnet for the USU and 
UT drive plate was measured. The value of B can be used to determine theoretical k; and 
kp for ideal coils and magnets. The calculated value of kj and kg were then compared 
with the measured value and efficiency of the drive plate can be evaluated. 

Strength of magnetic field of each magnet was measured using a Gauss meter. 
The meter used was a DC Gauss meter Model DCM Serial 1487 manufactured by 
AlphaLab, Inc. Strength of magnetic field generated by all four magnets for the USU 
drive plate and UT drive plate No. 9 was measured. The measurement was performed at 
multiple positions on the cross sections of the north and south poles of the magnet. The 
values of B at different position varied slightly, however, the average B value for north 
and south pole of each magnet is differed by less than 10%. Therefore, an average B of 


all four magnets was used to represent the magnetic field of the drive plate. For the USU 
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drive plate, the average B was equal to 733 Gauss. For the UT drive plate No. 9, the 


average B was equal to 1463 Gauss. The B value for the UT drive plate No. 9 was about 
twice the value for the USU drive plate. This as expected due to the different magnet 
material used for each drive plate (Alnico 5 for the USU drive plate and Alnico 8 for the 
UT drive plate). These values were used for calculating the theoretical kj and kg for an 
ideal coil and magnet for each drive plate. 

To calculate the theoretical value of k; and kg, Equation 5.13 and 5.18 are also 


presented in this chapter as shown below: 


k; = N,N,1,Bd, (6.12) 
and 
2 BA,N,N 
k, =—_———"—_,, (6.13) 
d sy 
where B = magnetic field of the magnet, 


Nx = number of loops in the coil in x-z plane, 

Ny = number of loops in the coil in y-z plane, 

L. = circumference length of the coil loop, 

A, = cross section area of the coil loop in x-z plane, 

y = width of coil loop in y direction, and 

d = moment arm of the drive plate. 

The cross section area and coil in x-z plane and y-z plane and the dimensions of 

coil are presented in Figures 6.20 and 6.21 for the USU and UT drive plate. Average 
circumference length of the coil loop and cross section area in x-z plane were determined. 


The number of loops in the coil in x and y direction were counted. The calculation of 
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theoretical of k; and kg was then performed based on this information and results are 


tabulated in Tables 6.2 and 6.3. It is important to note that unit of B is in Gauss and 
voltage is also in Joule/Coulomb, which are the SI unit, therefore in calculation, the unit 
conversion must be considered. 

The theoretical value of k; and kg were calculated and compared with the 
measured value obtained in previous sections as shown in Table 6.2 and 6.3. The 
theoretical k; values are equals to 0.426 and 0.828 for the USU and the UT drive plate. 
And, the theoretical kg values are equals to 1.178 and 2.290 for the USU and the UT 
drive plate. The different between the theoretical and experimental value can be 
considered as an efficiency of the drive plate. The efficiency represents how well the 
drive plate performs in comparison with the theoretical prediction. The kj, efficiency 
represents how much torque can be generated in an ideal coil and magnet compared to 
with the actual case, therefore, the k;, efficiency represents torque efficiency of the drive 
plate. The torque efficiencies of 75% and 80% were obtained for the USU drive plate 
and the UT drive plates, respectively, as shown in Table 6.2. 

For the case of kg, the efficiency represents back emf efficiency of the drive plate. 
In this case, the measured value of kg increases as the frequency increases. For this 
study, the kp values in frequency range of 0-200 Hz is considered because the resonant 
frequencies of most soil type fall in this range. Therefore, a back emf efficiency was 
calculated for frequency range of 0-200 Hz as shown in Table 6.3. The back emf 
efficiency of 32-50% was obtained for the USU drive plate and 38-42% was obtained for 


the UT drive plate. 
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Figure 6.20 Cross sections of coil for the USU drive plate 
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Figure 6.21 Cross sections of coil for the UT drive plate No. 9 
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Table 6.2 Calculation of Torque Efficiency 


Drive Plate ki o kbA o O 
Ideal | Measured | Efficiency 
t ae G 
3 (or eels ms) = PET ses 


Teena CORDS eT an 


Table 6.3 Calculation of Back Emf Efficiency 


kp (V/rad/s) 
Drive : 
Plate Gn’ (in an (eas) (Gauss ari a o 
B ‘0 


Coo [ass oss| w | so | ns [ii| osmose [moose 
arros ioa nsis a | a [es aao] onsas [se 


Note: 1. The kg values for frequency range of 0-200 Hz. 


The efficiency for kg is much lower than for the case of k; because the calculation 
for kg is based on the assumption that the magnetic field is uniform over the cross 
sectional area of the coil loop. If the torque efficiency is used as the ground truth, the 
actual cross section area with a uniform magnetic field is only about 50% of total area as 


the back emf efficiency is about 50% less than the torque efficiency. 


6.4 SUMMARY 

In this chapter, the measurements of parameters of the proposed electromagnetic 
model are presented. The experiment methods were developed to determine the variation 
of the parameters with frequency. 

The coil resistance, R, shows a slight variation with frequency. This is very small 
variation and can be neglected. The average coil resistance for all 8 coils is 45.6 ohms 


for the USU drive plate. 
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The coil inductance, L, is constant for the frequency range of 0 to 400 Hz. A 


value of 0.0256 H was measured for the USU drive plate. 

The torque-current factor, k; measured from both low frequency and static 
condition are in agreement. Both methods of measurement are unable to evaluate the 
variation of k; with frequency. In this study, the factor k; is assumed to be a constant and 
equal to 0.318 Ib-ft/A for the USU drive plate. 

The crosstalk effect was found in the back emf-rotational velocity, kg 
measurement. The crosstalk needs to be taken into account to determine the kg values. 
The factor kg increases from about 0.4 up to 1.4 V/rad/s for frequencies from 0 to 400 
Hz. The uncertainty of kg values measured from different calibrated specimen is shown 
in Figure 6.22. The curvefitting equation for the average kg is used and the uncertainty is 
evaluated in Chapter 7. 

The theoretical value of kg and k; can be predicted using the equations derived in 
Chapter 5. The theoretical k; values are equals to 0.426 and 0.828 for the USU and the 
UT drive plate. And, the theoretical kg values are equals to 1.178 and 2.290 for the USU 
and the UT drive plate. The torque and back emf efficiency represents the different 
between the measured and theoretical value of kj and kg. The torque effieciency for both 
drive plate was in the range of 75-80%. And, the back emf efficiency for both drive plate 
was in the range of 32-50%. The lower back emf efficiency is due to the lack of a section 


uniform magnetic field in the coil cross section. 
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Figure 6.22 Uncertainty for the kg measurement of the USU drive plate 
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CHAPTER 7 


VERIFICATION OF ELECTROMAGNETIC MODEL 


7.1. INTRODUCTION 

In this chapter, measured equipment generated damping is compared to the 
predicted equipment generated damping using the electromagnetic model developed in 
Chapter 5. The model parameters used for these predictions were described in Chapter 6. 
The uncertainty in predicting equipment generated damping from the model parameters is 
evaluated. The model predictions for the USU drive plate are compared with the UT 
drive plate. Procedures to apply the electromagnetic model to RC and TS testing are 
presented. An explanation of discrepancies in the literature between predictions of the 


equipment generated damping in RC and TS test is presented. 


7.2 | MEASURMENT OF BACK EMF EFFECT 

The electromagnetic model derived in Chapter 5, predicts that the back emf 
effects the stiffness and the damping of the system. The back emf effect tends to increase 
the measured resonant frequency and damping in the RC test and causes an error in 
torque measurement in the TS test. The effect of back emf on the TS test was found to be 
very small and negligible. This is discussed later in this chapter. This study primarily 
focuses on the back emf effect in the RC test. To verify the electromagnetic model 
prediction developed in this study, the back emf effects on the resonant frequency and 
damping need to be quantified. 

The first step in quantifying the effects of back emf on resonant frequency and 


damping is to measure those properties without the back emf effect. One method that can 
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be used to eliminate the back emf effect is to perform the RC test with an open drive 


circuit. This is the same approach that was employed in Chapter 4. The drive plate 
calibration was performed on the calibration specimens with free vibration and an open 
drive circuit. The resonant frequency and material damping of each calibration specimen 
is reported in Table 4.4. These results are repeated in Tables 7.1-7.4 to show the resonant 
frequency and damping can be effect. A slight effect of added mass on material damping 
can be observed in Tables 7.1-7.4. Adding mass tends to slightly decrease the material 
damping of the calibration specimens. However, the damping values of the calibration 
specimens are very small, about 0.02%. 

Once the stiffness and material damping of the calibration specimen was 
determined, the effects of back emf were measured. Typical RC tests using a closed drive 
circuit were performed by applying a constant force vibration over a band of frequency 
spanning the resonant frequency. The vibration response of the specimen at first mode 
resonance was evaluated. Because this experiment was performed with a closed drive 
circuit, the back emf effect was included. These tests were performed on each calibration 


specimen with added masses in the same manner as the drive plate calibration. 


Table 7.1 Summary of Back Emf Effect for 14-Hz Specimen 


Ko = 17.942 Ib-ft/rad 


Ts Without Back Emf With Back Emf Af De 
es m ei 
fm (Hz) D (%) f m (Hz) D/ (%) (Hz) (%) 


No Added 14.314 | 0.020 14.321 0.739 | 0.0075 | 0.719 
Mass 
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Table 7.2 Summary of Back Emf Effect for 62-Hz Specimen 


Kg = 343.80 Ib-fi/rad 


fa (Hz) | Dæ% sis a 
Meee, 61.861 0.020 | 61.887 0.0262 | 0.158 


Test 


Mass No.3 54.064 0.011 


Table 7.3 Summary of Back Emf Effect for 128-Hz Specimen 


Ke = 1560.60 Ib-ft/rad 
Without Back Emf With Back Emf Af, 


Dex 
fa (Hz) D (%) f m (Hz) D' (%) (Hz) (%) 
es 128.418 | ooo | 128.474 | 0.083 | 0.0558 | 0.073 

122.718 0.008 122.758 0.074 0.0400 0.066 


Table 7.4 Summary of Back Emf Effect for 188-Hz Specimen 


Ko = 3768.53 Ib-ft/rad 


Without Back Emf With Back Emf Af D 
Test 7 m %) 
fn (Hz) D (%) f (Hz) | DA (Hz) (%) 


a 188.079 | 0.023 | 188.133 | 0.060 
ass 

Mass No.1 | 180.775 | 0.019 | 180.816 | 0.052 
Mass No.2 | 173.917 | 0.017. | 173.959 | 0.050 
167.797 | 0.015 167.838 0.051 0.0406 
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Figure 7.1 shows the comparison between the response curves from tests with and 


without back emf effects on the 62 Hz calibration specimen. The response curve with the 
back emf effect clearly shows a much broaden shape hence higher damping. The 
resonant frequency for the test including the back emf is also slightly higher than the test 
without the back emf. Similar effects were observed for all of the other calibration 
specimens with and without added mass. Results of all the these tests are tabulated in 
Table 7.1-7.4. Note that symbol f'm and D’ represent the measured resonant frequency 


and damping that included the back emf effect. 


E = 61.887 Hz [622 Specimen 
£ = 61.887 Hz, D = 0.176% 62-Hz Specimen 


Normalized Amplitude 


Without Back emf 
fa = 61.861%, D = 0.018% 


61.6 61.7 61.8 61.9 62.0 62.1 


Frequency (Hz) 


Figure 7.1 Typical back emf effect on the response curve in RC test 
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The increase of resonant frequency, Afm and the equipment generated damping, 


Deq, are obtained from the difference between the resonant frequency and damping with 
and without back emf effect. The Afm represents the additional stiffness created from the 
effect of the back emf. This increase in stiffness is very small, less than 0.05%. The Deg 
decreases as the resonant frequency increases. The increase of stiffness and damping are 
observed for all specimens. The variation of the Afm and the Deg with resonant frequency 
are plotted in Figure 7.2. 

The back emf effect on resonant frequency and damping of the UT drive plate 
were also evaluated. Similar tests were performed using the same calibration specimens 
and the values of Afm and the Deg were determined. Figure 7.3 presents the plots of Afm 
and the Deg with resonant frequency for both the USU drive plate and UT drive plate No. 
9. In addition, the Deg of UT drive plate No.4 measured by Hwang (1977) is also 
included in Figure 7.3. As shown in Figure 7.3, values of Deq for UT drive plate No.4 
measured Hwang are very similar to those measured on drive plate No. 9 it this study. 
The Afm of the UT drive plate shows an increasing trend except for the one data point at 
about 14 Hz. The values of Afm and Deg for the UT drive plates are much higher than the 
values observed for the USU drive plate. The higher values are explained by the different 
magnet material thus the different value of k; and kg, as presented in Chapter 6. 

In the E TA model, the effect of back emf is presented in terms of keq 
and Ce}. The term keq and Ceq represents the equipment spring constant and equipment 


viscous damping coefficient due to the back emf. The measured keg and Ceq can be 
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Figure 7.2 Effect of back emf on the USU drive plate 
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0.20 Afan = 0.25001-0.303566 0009 7167fm 
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m 
-0.011868f 
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Figure 7.3 Comparison of back emf effect for the USU drive plate and the UT drive 
plate 
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evaluated from the increases of resonant frequency and equipment generated damping, 


Deq presented in Tables 7.1-7.4. The measured Keg and Ceq are calculated from equations 


shown below: 


where 


2 
/ 
ka [e ZJ-K,, and (7.1) 


vi—2D” 
f’ 


Cog = Daq ss -= ? (7.2) 


Ko = torsional spring stiffness of calibration specimen 


2J = total mass polar moment of inertia of the system. 


The measured values of keg and Ceq calculated using Equations 7.1 and 7.2 for the 


USU and the UT drive plates are compared with the predicted values of keg and Ceq in 


next section. 


7.3 


ELECTROMAGNETIC MODEL 


7.3.1 Model Prediction 


The predictions were made of the equipment spring stiffness, keg, and equipment 


viscous damping coefficient, Ceg, due to the back emf using the electromagnetic model 


developed in this study. The terms keg and Ceq are parameters in the transfer function for 


a voltage measurement derived in Chapter 5. This transfer function is: 


k, 


1 


R,+joL, 


He (7.3) 
i [(k+k,,)+ja(c+c,)—07J] 
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okk,L, 
“Rae (7.4) 
Coq = Ekki. (7.5) 


where @ = 2nf. 

Four model parameters; Re, Le, ki, and kg can be used to predict the values of keq 
and Ce}. These parameters were measured using procedures described in Chapter 6. 
Results of the measurements of model parameters are summarized in Table 7.5 for the 
USU drive plate. Using these parameters, and Equations 7.4 and 7.5, values of keg and Ceq 
can be calculated for the back emf effect at any loading frequency. 

Model parameters were measured at frequencies up to 400 Hz, however, 
verification of keq and Ceq is limited for frequencies below 200 Hz. This limitation is due 
to the lack of the measured values of keq and Ceq at frequency higher than 188 Hz. 


Predicting of the back emf effect for frequency above 200 Hz requires further study. 


Table 7.5 Measured Model Parameters for the USU Drive Plate 


requency 


Coil inductance, 0.0256 


Upper Bound 
kp = 0.39908+2.1769e-6f 2”! 
Back emf-rotataional . : Average 
velocity Factor, ke Virad/s | Increasing | 4. = 0.3776+3.12276-7£ 5" 
Lower Bound 
kg = 0.35126+1.7291e-8f °” 
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Figure 7.4 shows the comparison of measured and predicted values of Keg and Ceq. 


The upper and lower bounds in the prediction represent the uncertainty in the 
measurement of kg. The measured and predicted values of keg increase with frequency in 
similar fashion. However, the model prediction underestimates the value of keq at 
frequency between 60 and 130 Hz. This discrepancy is of little concern because the 
effect of back emf on resonant frequency is very small, less than 0.05%. The measured 
values of Ceq scattered and reasonably constant over the frequency range of 0-130 Hz, 
with a slight decrease at higher frequency. The prediction shows slightly underestimates 
Ceq for frequencies below 60 Hz. 

There are two possible explanations for the underprediction of Ceq values using the 
electromagnetic model. First, the model was developed based on the assumption that the 
equipment generated damping is solely the product of the back emf effect. In actual | 
system, additional energy losses can be created from other sources. Typically, these 
additional energy losses only have a very small effect on damping, thus it can be 
neglected. Second, the variation of k; with frequency is uncertain. As discussed in 
Chapter 6, the measurement method used for measuring the factor k; was unable to obtain 
the variation of k; with frequency and the parameter was assumed to be a constant with 
frequency. To evaluate a possible variation of k; with frequency, Equation 7.5 was used 
to back calculate possible values of k;. These values are plotted in Figure 7.5. This 
figure shows that a slight variation in kj values will improve the agreement between the 


measured and predicted Ceq. 
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Figure 7.4 Comparison of the measured and predicted values of keq and Ceg for the USU 
drive plate 
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© k=0.318 measured in Static Condition 


Back Calculation 
® Maximum k; 
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@ Average k; 
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0 50 100 150 200 


Frequency (Hz) 


Figure 7.5 Back calculated ki from measured Ceq using the electromagnetic model 


A summary of the model parameter for the UT drive plate is presented in Table 
7.6. The comparison of the measured and predicted keq and Ceq of the UT drive plate is 
also plotted in Figure 7.6. For this drive plate, ceq was again slightly underpredicted at 
low frequencies. The model performs very well at higher frequencies in predicting both 
Keg and Ceg. The error in predicting Ceq values at low frequency for the UT drive plate 


may come from the sources as discussed earlier. 


7.3.2 Model Verification 


The electromagnetic model is verified by comparing the model prediction of Afm 
and Deg with the measured values. The primary focus of this study is the prediction of 


Deg because the Afm is very small and can be neglected. The measured values of Afm and 
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Table 7.6 Measured Model Parameters for the UT Drive Plate 

Parameter Unit Variation with Measured Value 
Frequency 

ee je 


Coil Resistance, R, 


Torque-current Factor, k; | [b-ft/A Uncertain 


Back emf-rotataional 
velocity Factor, kg 


Virad/s | Increasing | kg = 0.87919+6.7644e-8f ©” 


Deq were presented earlier in this chapter. The predicted Afm and Deg are determined 


from the predicted values of keq and Ceq presented earlier. The predicted values of Afm 


and De, can be determined from Equations 5.34 and 5.35 which are also presented below: 


1} /K,+k K 
Af. =— 2e 2 D+D.) —,/—21-2D? |, and 7.6 
m N XJ D+Da) -Y3 5 (T0) 


Dy =—*- (1.7) 


In this case, the torsional stiffness and damping of calibration specimen and polar 
moment of inertia of the system are known therefore, the Afm and Deg can be evaluated 
with frequency from Equations 7.6 and 7.7. 

Figure 7.7 shows the measured and predicted effects of back emf on Afm and Deg 
with resonant frequency for each calibration specimen, with and without the added mass 
for the USU drive plate. The values presented in Figure 7.7 are the average predicted 


value and the error bars represent the uncertainty in the model prediction. 
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Figure 7.6 Comparison of the measured and predicted values of keg and ceq for the UT 
drive plate No. 9 
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An exponential function was fitted to the measured and predicted values. Predictions of 


Afm are not especially good, however, this is of little concern because Afm is very small. 
The model accurately predicts Deg. It slightly underestimates the values of Deg for 
resonant frequencies between 14 and 62 Hz, and overestimates Deq for resonant 
frequencies between 128 and 188 Hz. However, these errors are inconsequential. 

The model performs very well in predicting the Afm and Deg for the UT drive plate 
as shown in Figure 7.8. In these tests, large amounts of ambient noise corrupted low- 
frequency measurements. Ifthe Afm measurement at 14 Hz is ignored, the model 
predicting Afm performs better for the UT drive plate than the USU drive plate. The 
model underpredicts the Deg at 14 Hz by about 0.7%, however, the matches between 
measured and predicted values is very good for the higher frequencies. One possible 
explanation for the better performance of the model for the UT drive plate is that back 
emf effects are larger for the UT drive plate, so small errors and inconsistencies in the 


model have less effect. 


7.3.3 Model Application for RC and TS Test 

As explained in Chapter 5, different methods should be used to account for the 
back emf effect in RC and TS tests. The electromagnetic model prediction presented 
earlier was verified for the RC test. The model can also be used for the TS test. In this 
section, an example of the application of electromagnetic model to account for the back 


emf effect in RC and TS is presented. 
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Figure 7.7 Electromagnetic model prediction of Afm and Deg for the USU drive plate 
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Figure 7.8 Electromagnetic model prediction of Afm and Deg for the UT drive plate No. 9 
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For the RC test, the increasing of resonant frequency due to the back emf effect is 


very small and can be neglected. The electromagnetic model is mainly employed for 
predicting the back emf effect on damping measurements. As explained earlier, the 
model presents the back emf effects on damping in term of Ceg. The value of Ceq can be 
predicted by using measured model parameters and Equation 7.5. The equipment 
generated damping due to the back emf effect, Deg, is evaluated from a ratio of Ceq and 
critical damping coefficient, ce, presented earlier in Equation 7.7. ce is a function of 
resonant frequency, fm, and total mass polar moment of inertia of the system, XJ, as 
presented below: 

c, =4af, XJ, (7.8) 


where f = fn 


" V1-2D7 
Equation 7.8 implies that the Deg is not only a function of resonant frequency but also a 
function of XJ. Therefore, the different damping corrections should be applied for 
different sized specimens, and different drive plates and top platens. Equipment 
generated damping can be calculated by determining Ceq from the electromagnetic model, 
ce from the resonant frequency of the specimen and moment of inertia of the drive plate 
and top platen. Measured damping, D’ can be corrected using: 
Corrected D = D’ — Dea (7.9) 

In addition, the model can be used to predict the response curve that included the 

back emf effect. Equations 7.3-7.5 are used to generate a transfer function for a voltage 


measurement in RC tests. The response curve is a plot of this transfer function versus 
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Figure 7.9 Predicting response curve for RC test using the electromagnetic model 


frequency. Figure 7.9 shows an example of model prediction of response curve for the 
RC test on the 62 Hz calibration specimen. 

For the TS test, loading frequency is much lower than in the RC test, hence, the 
rotational velocity of the magnet is much lower. As described in Chapter 5, torque 
credie in the RC/TS drive system is partially transmitted to the dashpot, thus the actual 
torque applied to the drive plate needs to be corrected. The measured torque can be 
corrected by using Equation 5.36, which is repeated here: 

T, =T,, -(K,,9 + ¢,,88) = k,0 +080 + JS°0 , (7.10) 
where Ts = corrected torque due to the back emf, 


Tm = measured torque, 
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8 = rotational angle, 
kg = torsional spring stiffness of soil specimen, 
Ce = torsional coefficient of viscosity of soil specimen, 
J = mass polar moment of inertia of the system, and 
s = j@ (@ = 2xf). 

Because the loading frequency is very low in the TS test, the ceq has a very small 
effect to the applied torque. Figure 7.10 shows the measured and corrected torque 
calculated from Equation 7.10 for a TS test on soil specimen at 0.333 Hz. Both curves 
are virtually identical, therefore the stiffness and damping calculated from these torques 


are exactly the same. As a result, the effect of back emf in the TS test is very small and 


can be neglected. 


1.0 
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Figure 7.10 Effect of back emf on TS test 
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7.3.4 Model Limitations 


The model proposed for this study is limited to predict Deg at frequencies less than 
200 Hz. The uncertainties in the model parameters seem to have a very small effect in 
predicting Deg, The model performs very well in predicting Deg for this range of 
frequency. It should be noted that although the model is limited for frequencies less than 
200 Hz, the values of Deg above 200 Hz are very small and most likely can be neglected. 
7.4 DISCUSSION ON DISCREPANCIES 

IN EQUIMENT GENERATED DAMPING 

IN RC/TS TESTS 

The measured Deg reported by Hwang (1997) decreases with increasing resonant 
frequency in RC tests. This same trend was observed in this study. However, in TS test, 
Hwang (1997) reported an increasing of Deg with loading frequency. This discrepancy is 
investigated in this section. 

In the TS test, the hysteresis loop is created from the plot of torque versus 
rotational angle. The damping is determined from area of hysteresis loop which is 
directly related to the phase difference between the torque and rotation. As explained in 
Chapter 5, torque created in the RC/TS drive system is proportional to the current and 
torque should be calculated from the current signal using the torque-current factor, ki. 
However, voltage measurements are more commonly used for determine torque in TS 
testing. To investigate the effect of using voltage rather than current in torque 
calculations, the correlation between current (i) and voltage (V) is shown below: 


V 


pea a, (7.11) 
R, + joL,) 
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Equation 7.11 implies that the voltage and current signals are not in phase with 


each other due to the coil inductance. This phase shift tends to increases as the frequency 
increases as shown in Figure 7.11. As a result, an addition of apparent damping is 
created when the voltage is used for calculating torque. Figure 7.12 illustrates the area of 
hysteresis loops due to the phase shift between the voltage and current. The increasing of 
loop area represents the error in measuring damping as the frequency increases when 
voltage is measured instead of current. Therefore, the phase shift of voltage and current 
is a source of error in measuring damping. If voltage is measured, a correction for phase 


shift is required when calculating torque. 


Phase Shift = tan” (oL /R,) 


Phase Shift (rad) 


— R, = 45 ohm, L, = 0.025 H 


0 20 40 60 80 100 
Frequency (Hz) 


Figure 7.11 Phase shift between voltage and current in RC/TS equipment 
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Figure 7.12 Variation of additional damping due to voltage measurement for torque 


This discussion explains the discrepancy of the increasing of equipment generated 
damping with frequency in the TS test as oppose to the resonant column test reported by 
Hwang (1997). The damping from the voltage and current phase shift for each frequency 
is calculated from the loop presented in Figure 7.12 and compared with the assumed 
equipment generated damping measured by torsional shear test using voltage for 


calculating torque as shown in Figure 7.13. The phase shift damping matches very well 
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Figure 7.13 Comparison of damping from a phase shift and measured equipment 
generated damping for TS test by Hwang (1997) 


with the Hwang’s experimental data. As a result, the damping reported by Hwang is the 
damping created from using torque calculated by voltage without a phase shift correction. 
The actual equipment generated damping created from the back emf is negligible in TS 


test. 
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CHAPTER 8 


| ACCOUNTING FOR NONUNIFORM SRESS-STRAIN IN RC/TS TESTS 


8.1 INTRODUCTION 

Nonuniform stress-strain effects occurring over the radius of soil specimen are 
one of the weaknesses in the RC/TS testing. An equivalent radius approach proposed by 
Chen and Stokoe (1979) has been used to account for this effect. As discussed in Chapter 
2, the equivalent radius approach uses an equivalent radius to calculate stress and strain 
when analyzing RC/TS test results. In this study, a method called the stress integration 
approach was developed to account for the effect of nonuniform stress-strain more 
precisely. The stress integration approach utilizes three assumed soil stress-strain models 
to develop representative stress-strain relationships for soils in RC/TS testing. 

This chapter presents the development of the stress integration approach. 
Procedures used for each assumed stress-strain soil model are described. Limitations of 
the equivalent radius approach are discussed. Finally, the equivalent radius approach was 
modified using the stress integration method to improve damping measurements. 

8.2 STRESS-STRAIN DISTRIBUTION 

IN SOIL FOR RC/TS TESTING 

A soil specimen subjected to torsional loading in a RC/TS test is analyzed as a 
cylindrical bar twisted by torque, T, acting at the end of the bar as shown in Figure 8.1. 
Circular cross sections of the bar are identical at every cross section and subjected to the 


same internal torque, T. Torque on a right-circular cylinder causes as condition of pure 
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a) Transverse section 
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b) Longitudinal section 


Figure 8.1 Distribution of shearing stress in soil column for RC/TS test 
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shear. This means that every cross section remains planar and circular and all radii 


remain straight. Shear strain at the outer surface or maximum shear strain, Ymax can be 


expressed as: 


Yax = —> (8.1) 


where 6 = rotational angle, 

R = radius of the bar, and 

L = length of the bar. 
Because the radii in the cross sections of the bar remain straight and undistorted during 
twisting, the shear strain within the interior of the bar is given by: 


mee 
TD ORI? 


(8.2) 
where r = interior radius of the bar. 

Equations 8.1 and 8.2 are the equations for the torsional deformations of circular 
bars that are based upon the geometry of the deformed bar, therefore they are valid for 
any material whether it behaves linearly or nonlinearly. The rotation, 0, along a 
longitudinal section of the soil column in a RC/TS tests varies linearly because of its 
uniform cylindrical shape. However, the distribution of shear stress over the radius of 
soil specimen is non-uniform as shown in Figure 8.1a, due to nonlinear behavior of soil at 


medium to high strain levels. This nonuniform stress strain distribution must be taken 


into account when analyzing RC/TS tests. 
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8.3 STRESS INTEGRATION APPROACH 


The stress integration approach is a new, more general approach to account for 
nonuniform stress-strain than the conventional equivalent radius approach (Chen and 
Stokoe 1979). The equivalent radius approach applies an equivalent radius ratio to 
calculate stress and strain corresponding to the measured torque (T) and rotation (6) in 
RC/TS tests. The nouniform stress-strain can be accounted for more precisely by 
integrating the stress over the radius of soil specimen to obtain a T-@ relationship. 

The stress integration approached developed in this study can be derived using 
basic equations of mechanics. The approach starts by relating the shear stress acting on a 


circular cross section to the applied torque, T, using the equation: 


R 
T= Í dM = ferda = [2nr?cdr ; (8.3) 
A 0 
where M = resultant moment over the entire cross section area, 


A = cross section area, and 
qt = shear stress. 
If the stress-strain relationship is known, then for any given value of 0, the shear 

stress at any point in the specimen can be determined. This requires using Equation 8.2 
to calculate the shear strain and calculating the corresponding shear stress from the 
known stress-strain relationship. Thus, the distribution of shear stress and strain over the 
entire cross section of the soil can be evaluated. Since the shear strain varies linearly 
with the radius, the distribution of the shear stresses has the same shape as the stress- 


strain relation. Torque then can be obtained from the integral is Equation 8.3. 
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There are three procedures that must be performed in order to apply the stress 


integration approach to RC/TS testing. First, an experimental T-0 relationship must be 
measured. Second, an assumed soil stress-strain model is applied to Equation 8.3 and a 
theoretical T- relationship is developed using either closed form integration or 
numerical integration. And third, curve fitting techniques are then used to match the 
theoretical T-6 relationship with the measured T-O relationship to obtain the best fit 
model parameters. These model parameters are then used to develop the stress-strain 
relationship for soil. 
8.3.1 Measurement of T-@ Relationship 

in RC/TS Testing 

In RC testing, a resonant frequency of the system and a corresponding peak 
rotation at resonance are measured. If the mass polar moment of inertia of the specimen 
and drive plate are known, the resonant frequency of the system can be related to the 
shear wave velocity, vs of the soil specimen obtained by using the BtanB method, as 
presented in Chapter 3. An equivalent linear elastic shear modulus or an effective shear 
modulus, Ger, is calculated from: 

Ga =PV, » (8.4) 

where p = mass density of soil specimen. 

By linearizing the problem, Ger and peak rotation at resonance, 8, can be related 
to an applied torque, T, using the equation: 


G,,1,8 
arta (8.5) 
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where I, = area polar moment of inertia, and 


L = length of specimen. 
Equation 8.6 is used to generate measured T-6 relationship in RC testing. 
In the TS testing, the T-O relationship is measured directly and a hysteresis loop is 
developed 2 the T-9 plane. The measured T-9 relationship in TS testing is generated 


from the maximum T and 9 in each loading level. 


8.3.2 Development of Theoretical T-0 Relationship 

The three soil stress-strain models used in this study were the hyperbolic model, 
the modified hyperbolic model, and the Ramberg-Osgood model. Each model is 
evaluated corresponding to theoretical T-6 relationship using either closed form or 
numerical integration. Closed form integration was used for the hyperbolic model and 
the numerical integration was used for the modified hyperbolic model and the Ramberg- 
Osgood model. Procedures for the closed form integration and numerical integration are 
presented in the following sections. 

8.3.2.1 Closed Formed Integration. For the hyperbolic model, a closed form 
solution is used to develop the theoretical T-0 relationships. The hyperbolic stress-strain 


relationship can be expressed as: 


T= =, (8.6) 
1+% 
Yr 
where Gmax = small strain shear modulus, and 


Yr = reference strain. 
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For a cylindrical specimen in torsion, the shear strain at any point within the interior of 


soil specimen is defined in Equation 8.2, therefore, shear stress can be related to the 


interior radius as follow: 


Guð E 
t=—— E. (8.7) 
1+—80— 
y, L 
Substituting Equation 8.7 into Equation 8.3 gives: 
R G nax9 2 
T= [2ar? —— ar. (8.8) 
lr 
2 1+—0— 
Y L 
By integrating Equation 8.8, the following equation for a T-9 relationship for the 
hyperbolic model is obtained: 
-1 E lee a2 41.3 oR 
T= : May —3R( - ) + 6( F yo [+ 2G 6 =] Iny Ing, + : ) |. (8.9) 


A flow chart showing the stress integration approach using the hyperbolic model 
is presented in Figure 8.2. A measured T-9 relationship from a RC/TS test is shown in 
Figure 8.2a. First, the model parameters Gmax and yr are assumed, and a theoretical T-0 
relationship is calculated using Equation 8.9. The theoretical T-O curve is plotted with 
the measured T-6 curve as Figure 8.2b. Curve fitting techniques are then used to match 
the theoretical and measured T-O curves by adjusting values of Gmax and yr iteratively. 


The best fit parameters are used to develop the stress-strain relationship presented in 


Figure 8.2c. 
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Figure 8.2 Flow chart showing the stress integration approach using the hyperbolic 
model 
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8.3.2.2 Numerical Integration. Numerical integrations are used to develop the 


theoretical T-9 relationships for both the modified hyperbolic model and Ramberg- 
Osgood model. The model parameters for the modified hyperbolic model are: Gmax, Yr 


and a. The modified hyperbolic stress-strain relationship can be expressed as: 


R (8.10) 


And, the model parameter for the Ramberg-Osgood model are: Gmax, Yr a, and b. The 


Ramberg-Osgood stress-strain relationship can be expressed as: 


b-1 
(8.11) 


Equation 8.10 and 8.1 1 cannot be integrated like the two parameter hyperbolic model. 


T T 
y=——|1+a 
G max G max Yr 
Therefore, a numerical integration is used for these models. Two different methods for 
numerical integration can be used. 
The first method involves integration relative to strain. From Equation 8.2, the 


interior radius, r, relative to strain, y, is: 


=—y, 8.12 
ea (8.12) 


The differential of strain, dy, can be related to the differential of radius, dr, using 


Equation 8.12 as follows: 
dt = <dy. (8.13) 


Equations 8.12 and 8.13 are substituted into Equation 8.3, and changing the upper limit of 


integration tO Ymax to correspond with the new independent variable, y, gives: 
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L 3 Yeon 
T= 2E) Foar, (8.14) 


where t(y) = t-y relationship. 

Figure 8.3 presents a flow chart illustrating the stress integration approach using 
integration relative to strain, y, for the modified hyperbolic model. A measured T-0 
relationship from a RC/TS test is shown in Figure 8.3a. The stress-strain relationship 


described by the modified hyperbolic model is substituted in Equation 8.14 to obtain: 
3 Yeux 
T= 2a{ = f Gmt ly? dy (8.15) 


Torque is determined by numerically integrating the stress-strain relationship as shown in 
Figure 8.3b. The resulting theoretical T-0 relationship is shown in Figure 8.3c. Curve 
fitting is then performed to determine the best fit model parameters. 

The second method involves integrating with respect to radius, r. Figure 8.4 
presents a flow chart illustrating the stress integration approach using integration relative 
to radius, r, for the Ramberg-Osgood model. A measured T-9 relationship from a RC/TS 
test is shown in Figure 8.4a. The Ramberg-Osgood model presents y as a function of 
shear stress, t. For a given rotation, 0, there is a one-to-one correspondence between y 
and r, shown in Equation 8.12. In the case of Ramberg-Osgood model, the radius, r, can 


then be expressed as: 


b-1 
| (8.16) 
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Figure 8.3 Flow chart showing the stress integration approach for the modified 
hyperbolic model using integration relative to strain 
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Figure 8.4 Flow chart showing the stress integration approach for the Ramberg-Osgood 
model using integration relative to radius 
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A plot of shear stress, t, and radius, r, can be generated using Equation 8.16. Torque is 


determined by numerically integrating the t-r relationship using: 
R 
T = 2n [u(ryrdr , (8.17) 
0 


where t(r) = t-r relationship. 
Figure 8.4b presents the integration of the t-r relationship. The resulting theoretical T-O 
relationship is shown in Figure 8.4c. Curve fitting is then performed to determine the 


best fit model parameters. 


8.3.3 Curve Fitting Technique 

Curve fitting techniques were performed using an assumed stress-strain soil 
model following the procedure illustrated in Figures 8.2-8.4. The program Igor Pro 4.01 
was used to perform the curve fitting. The details on curve fitting techniques used in Igor 
are presented in this section. 

In curve fitting, the theoretical T-0 relationship calculated from the assumed soil 
stress-strain model is fit to the measured T-@ relationship by adjusting values of the 
model parameters. For the hyperbolic model, the curve fitting function is a closed form 
solution. For the modified hyperbolic model and the Ramberg-Osgood model, the curve 
fitting function requires numerical integration. The Igor curve fitting program searches 
for the values of model parameter that best match the experimental data. The 
“best” values of the model parameters are the ones that minimize the value of chi-square. 


chi-square is defined as: 
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chi -square = J22) , (8.18) 


where y = fitted value for a given point, 
yi = original data value for the point, and 
oi = standard deviation of each data value 

Igor uses the Lavenberg-Marquardt algorithm to search for the minimum value of 
chi-square using nonlinear, least square curve fitting. The search starts from a supplied 
initial guesses of the model parameters. The initial guesses must be provided manually 
and the curve fitting algorithm searches for the best fit starting from these initial guesses. 
As the fit proceeds and better values are found the chi-square value decreases. The fit is 
finished where the difference between experimental data and fitted data are minimum. If 
the initial guess is poor, the fit may conclude with an unsatisfactory result. This can be 
observed visually. The initial guesses must be adjusted until the fit achieves the lowest 
chi-square value. 

The chi-square factor is an indicator of the quality of the fit. However, judgment 
must be exercised in evaluating the quality of the fit. In some cases, two models may 
return similar values of chi-square, however, one model may match the experimental data 
very well at low to medium strain levels and the other may match better at high strains. 
In this case, neither of the models are good for predicting the soil behavior over the entire 


range of strains and special considerations must be taken. 
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8.4 EVALUATION OF THE EQUVALENT 


RADIUS APPROACH FOR PREDICTING 

STRESS-STRAIN RELATIONSHIP FOR SOIL 

The equivalent radius approach is based on the assumption that the representative 
stress and strain in soil specimen in torsion occurs at a radius called the equivalent radius, 
re To apply the equivalent radius approach, an equivalent radius ratio, Reg, which 
defined as the ratio of the equivalent radius and the outside radius of soil specimen is 
used to calculate the stress and strain. According to Chen and Stokoe (1979), Req value 


varies from 0.82 for strains below 10° % to 0.79 for strains at 10! % for a solid 


specimen. Shear strain is calculated from: 


Y= Roq X Yma - (8.19) 
Shear stress, t, is calculated using: 
T=R a X Trex (8.20) 
where Tmax = Maximum shear stress for a linear elastic specimen: 
t= E, (8.21) 


P 
To evaluate the equivalent radius approach, theoretical stress-strain curves were 
developed using three soil models. The thress models are the hyperbolic model, the 
modified hyperbolic model and the Ramberg-Osgood model. The model parameters 
assumed for each soil model and the stress-strain curves are presented in Figure 8.5. 
From the theoretical soil models presented in Figure 8.5, the T-O relationships 
were determined. The closed form solution was used calculate the theoretical T-O 


relationship for the hyperbolic model. Numerical integration was use to calculate the 
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theoretical T-6 relationships for the modified hyberbolic model and Ramberg-Osgood 


model. The theoretical T-O relationships for each soil model are presented in Figure 8.6. 

Once the theoretical T-6 relationships were developed, the equivalent radius 
approach was then applied to determine the stress-strain relationship using Equations 
8.19 and 8.20. The stress-strain relationship from equivalent radius approach was then 
compared with the actual stress-strain relationship to evaluate how well the equivalent 
radius approach predicts the stress-strain relationships. Those curves are presented in 
Figure 8.7. Results are plotted in terms of dimensionless or normalized stress and strain 
so that the results are independent of model parameter: Gmax and yr. Shear stress was 
normalized by tf, where tr is defined as: 

Te =G nax Yr- (8.22) 
And, shear strain was normalized by yr. 

Figure 8.7 shows that the equivalent radius approach slightly overestimates shear 
stress, T, at high strain levels for all three soil models. This overestimation indicates the 
equivalent radius approach has limitations when performing RC/TS tests at high strains. 
At strains below 0.1% the equivalent radius approach is satisfactory. 

The equivalent radius approach was developed based upon moduli calculations. This 
approach does not necessarily apply equally well for damping measurements. The 


following section presents an approach used for damping measurements. 
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Figure. 8.5 Three soil models used for evaluation of equivalent radius approach 
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Figure 8.6 Theoretical T-6 relationship of the soil models presented in Figure 8.5 
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Figure 8.7 Comparison of actual stress-strain curves and the stress-strain curves using 
the equivalent radius approach 


8.5 ACCOUNTING FOR NON-UNIFORM 

STRES-STRAIN IN DAMPING 

The nonuniform stress-strain effect on soil specimens in torsion must be taken 
into account when developing damping curve in RC/TS testing. In this study, the 
equivalent radius approach was modified to determine the strain corresponding to a 
measured damping value more precisely. Values of Reg based on both shear modulus and 
damping are determined separately. These results are compared and normalized curves 


are presented. 
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8.5.1 Procedure to Determine (Reg)c Based 
upon Shear Modulus 


Figure 8.8 illustrates a procedure to evaluate the (Req)g based on shear modulus. 
The hyperbolic model was used to develop a stress-strain relationship for a theoretical 
soil. This stress-strain relationship was then used to generated modulus reduction (G 
versus y) curve shown in Figure 8.8a. The stress-strain for the theoretical soil was used 
to calculate a T-9 relationship using Equation 8.9. With this T-6 relationship, the 
ien shear modulus, Ger can be calculated using Equation 8.5. The Geff versus 
O curve is shown in Figure 8.8b. For a given, yı, the corresponding value, G1, can be 
obtained from Figure 8.8a. A value 0; that is associated with the Gerr value equal to G; is 
identified in Figure 8.8b. Then using the following equation, a value of (Req)ci is found 


for 6; and yı as shown in Figure 8.8c. 


Ree (8.23) 


8.5.2 Procedure to Determine (Req)p Based 
upon Damping 


Figure 8.9 illustrates a procedure to evaluate (Req)g based on damping. Again, the 
hyperbolic model was used to develop a stress-strain relationship for a theoretical soil. 
As described in Chapter 3, hysteretic damping of soil can be obtained from the 
characteristics of the hysteresis loop. Assuming the theoretical soil behaves according to 


Masing rules, the damping can be determined from: 
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Figure 8.8 Determination of Req based on shear modulus 
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Figure 8.9 Determination of Req based on damping 
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2f t(y)dy 
p- 1l. (8.24) 
m| yy) 


Equation 8.24 was used to develop the D-y curve can be generated as shown in 
Figure 8.9a. The stress-strain relationship for the theoretical soil was then used to 
calculate T- 0 relationship using Equation 8.9. This T-@ relationship was then integrated 
using: 


2{T@)d0 


2 
D. = 4| “———- 1]. 8.25 
eff 1 6T(@) ( ) 


With this T-9 relationship, the effective shear modulus, Deff can be calculated 
using Equation 8.25. The Der versus O curve is shown in Figure 8.9b. For a given, yı, 
the corresponding value, D;, can be obtained from Figure 8.9a. A value 0; that is 
associated with the Deg value equal to D; is identified in Figure 8.9b. Then using 


Equation 8.23, a value of (Reg)p1 is found for 0; and yı as shown in Figure 8.9c. 


8.5.3 Instability in Closed Form Solution 

Instability of the closed form solution for the T-9 relationship was observed at 
low values of 8. This occurred when shear modulus was in the linear range where Geff 
was equal to Gmax as shown in Figure. 8.10. The instability of the closed form solution 
results in unrealistic values of the calculated Ge where Ges becomes larger than Gmax 


which is theoretically impossible. Consequently, an error can be created in the evaluation 
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Figure 8.10 Instability of the closed form solution 


of Reg ratio due to this instability. As shown in Figure 8.10, the instability occurs when 
value of Geg/Gmax is approximately one. For these low strains, the value of Reg is 
independent of strain level. In this study, Reg was analyzed only in the range of usable 


curve data. Reg was assumed to be constant for strains in the linear range. 


8.5.4 Evaluation of Effects of Soil Nonlinearity on Req 

Effect of soil nonlinearity on Reg based on both G and D was investigated. Three 
hyperbolic soil models were used. These soil models have reference strains, y,, of 0.01%, 
0.1%, and 1%, and have same value of Gmax, 1.0E6 lb/ft”. The modulus E T E 
and damping curves for all soil models are presented in Figure 8.1 la and Figure 8.12a, 
respectively. The Reg ratio curves based on G and D for these soil models were 


developed following the procedure shown in Figures 8.8 and 8.9. 
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Figure 8.11 Modulus curve and their associated Req curves based on G with different 
soil models 


Reproduced with permission of the copyright owner. Further reproduction prohibited without permission. 


192 


Hyperbolic Model (G nax = 1.0E6 psf) 


y, = 0.01% 
— y, =0.1% 
en Y =1% 
9 
A 
A y (%) 
a) Damping curves 
0.85 
0.80 
7 
ne 
0.75 
0.70 
0.65 
10° 10° 10° 10° 10° 107 
b) R,, curves oe) 


Figure 8.12 Damping curve and their associated Reg curves based on D with different 
soil models 
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Figure. 8.13 Req curves based on both G and D 
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(Req) curves based on G are presented in Figure 8.11b. All (Req)g curves 


converge to a value of about 0.8 in their linear range. As the strain and rotation increase, 
values of (Req)g decrease. The points where these curves decrease associated with the 
value of yr. 

(Req)p curves based on G are presented in Figure 8.12b. Again, all (Req)p curves 
converge to a value of about 0.8 in their linear range. The curves also fall off similarly to 
the (Reg) curves but decrease to lower values. 

It is more useful to present the results of Req in terms of normalized rotation, 0/0,, 
so that results are independent of the soil nonlinearity. The reference rotation, 0,, is 


defined as: 


L 
o, = 6.26) 
All Req curves based upon G and D, are plotted versus the 0/8,, and the results are 
presented in Figure 8.13b. Now all of the (Req)g curves based on G plot on top of each 
other, and all of the (R.q)p curves based upon D plot on top of each other. All curves 
converge to a value of about 0.8 at low strains. At high strains, the (Reg)p curves are 
significantly lower than the (Reg)c curves. This suggests that using the same Reg values 
for both G and D is not appropriate. | 
Figure 8.14 presents a comparison of Reg values generated in this study with the 


Req values suggested by Chen and Stokoe (1979). Three weaknesses with the approach 


suggested by Chen and Stokoe (1979) are observed in Figure 8.14. First, the soil 
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0.85 
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Req 


0.70 


0.65 


—-—- Equivalent radius approach (Chen and Stokoe, 1979); 


0.60 
0.0001 0.001 0.01 0.1 1 


y) 


Figure 8.14 Comparison of Req values from this study and Chen and Stokoe (1979) 


nonlinearity is not taken into account. Second, a single value of Reg does not accurately 
represented the actual Req for a wide range of strains. Third, values of (Req)p based upon 


damping are significantly lower than the values of (Reg)g based upon G at high strains. 


8.5.5 Using Modified Equivalent Radius Approach 
to Generate Damping Curves 


Figure 8.13b can be used to develop damping curves from RC/TS tests. To do 
this, reference strain, Yr, must first be determined. This value would typically be 
identified during curve fitting to generate modulus reduction curves. Using the yr, 
reference rotation, 0,, can be calculated using Equation 8.26. For every rotation 


measured in RC/TS testing, a value of (Reg)p can be obtained from Figure 8.13b. Strains 
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for each rotation can be determined using Equation 8.18. In the TS testing, the damping 


associated with that strain can be obtained from the properties of the T-O hysteresis loop. 
In the RC testing, the damping associated with the calculated strain can be determined 
from either the half power bandwidth method or free vibration decay method. 
8.5.6 Limitations of Modified Equivalent 

Radius Approach 

Three limitations in using hysteretic damping to describe the material damping for 
soil are recognized. First, the hysteretic damping is zero at very small strains because the 
area of hysteresis loop becomes zero in the linear range. Actual material damping of soil 
has never becomes zero due to damping from other mechanisms. Second, the hysteretic 
damping overestimates the damping for soil at high strains. The hysteretic damping 
increases up to a maximum value of 63.7 % at strain levels above 50%, which does not 
happen in the real soil behavior. Third, the hysteresis loop is governed by the Masing 
behavior which does not accurately describe the behavior of soil. These limitations 


require further study. 
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CHAPTER 9 


RC/TS TESTING OF BONNEVILLE CLAY 


9.1 INTRODUCTION 

RC/TS testing was performed on Bonneville clay to develop modulus reduction 
curves and damping curves. Previous to this study, no measurements have been made of 
modulus reduction and damping in Bonneville clay. In this testing, the electromagnetic 
model was used to correct for equipment generated damping and the stress integration 
approach was used to calculate stresses and strains. This chapter is separated into four 
sections. Each section is described below. 

The first section presents a description of soil specimens used in this study. Soil 
description includes sample depth and location, physical properties, and soil 
classification. 

The second section presents the RC/TS testing procedures. These procedures 
include the specimen preparation and RC and TS testing sequences. 

The third section presents the RC/TS analysis procedures. These procedures 
include correction for equipment generated damping, and application of the stress 
integration approach to develop modulus reduction and damping curves. 

And the fourth section presents a discussion of the factors affecting the modulus 
reduction and damping curves for Bonneville clay. These factors are plasticity index, 
confining pressure, and a cyclic degradation. The modulus reduction and damping of 


Bonneville clay are compared with generic modulus reduction and damping curves. 
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These comparisons show that Bonneville clay behaves more linear than generic curves 


predict. At the end of this section, conclusions drawn from the test results are presented. 


9.2 SOIL SPECIMENS 

Bonneville clay is very soft lacustrine silty-clay deposited in Ancient Lake 
Bonneville. It underlies the shores of both the Great Salt Lake and Utah Lake as well as 
the large portions of Salt Lake Valley and adjacent valleys. Seven samples of Bonneville 
clay were tested in this study. They were obtained from four different sites distributed 
across the Wasatch front. Three samples were collected at a site near the Provo River in 
Utah County. One sample was from a site at Nibley in Cache County. One sample was 
from a site in Logan in Cache County. And, two samples were collected from a site near 
the Salt Lake City Airport in Salt Lake County. Locations of each site are presented in 


Table 9.1. 


Table 9.1 Sample Sites Location 


Geneva Rd. Bridge over 
Cache 2835 S 1000 W 
Nibley, UT 
USU Drainage Farm Site 


BYU Research Site North 
Salt Lake West of Salt Lake City 

International Airport, 

Salt Lake City, UT 


2 
4 
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Undisturbed samples were recovered from the sites using thin-walled Shelby 


tubes. Radiographic (X-ray) images were made of all samples prior to opening sample 
tubes to assess the sample quality and to select portions of the sample for testing. 
Atterberg limits and hydrometer tests were performed on each sample to determine the 
soil classification and index properties of soil specimen. Other physical properties of the 
soil specimen were also determined prior to the RC/TS testing. The physical properties 
and unified soil classification are presented in Table 9.2. 

All of the soil specimens in this study are fine-grained soils with, the percent 
passing the No. 200 sieve all greater than 85%. The index properties of each soil 
specimen are plotted on the plasticity chart as shown in Figure 9.1. Specimens from the 
site near the Provo River, PRV37 and PRV67, are low to medium plasticity silty clay 
(CL), except for one sample, PRV93, from the depth of 93 ft which was clayey silt with 
very fine sands (ML). A specimen from shallow depth, SLC14, from a site near the Salt 
Lake City airport is a low to medium plasticity clay (CL) and a specimen from a deeper 
depth, SLC35, was low plasticity silty sand (ML). Soil specimens from Nibley, NBL24, 


and Logan, LGN14, are high plasticity clay (CH). 
9.3 RC/TS TESTING PROCEDURE 


9.3.1 Specimen Preparation 
The X-ray images of the sample tube were visually examined to select the best 
quality portions of the samples. Once a sample portion was chosen, each end of the 


selected portion was cut using a tube cutter. The side of the tube was then cut vertically 
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Table 9.2 Physical Properties of Soil Specimens 


% Clay Size Pegres 


Void of 
Ratio | Saturation 
%) 


Sample 
ID : 
PRV67 
PRV93 
NBL24 
LGNI4 
SLC14 
SLC35 


Unified Soil iqui Plasticity % Passing 
Classification imi No.200 Sieve 


EUDES < 0.002 Content | Weight 


TO 


E 89.94 =a 117.15 | 0.735 | 76.48 


007 


201 


Sample ID - USC Classification 
PRV37 - CL 
PRV67 - ML 


PRV93 - CL 
NBL24 - CH 
LGN14 - CH 
SLC14 - CL 
SLC35 - ML 


Plasticity Index 


MH or OH 


ML or OL 


0 20 40 60 80 100 
Liquid Limit 


Figure 9.1 Unified soil classification of Bonneville clay specimens tested in this study 


using a band saw. In most case, the tube sprung open after cutting, and the soil was easily 
removed. When the tube did not spring open, the opposite side of the tube was cut to 
remove soil with minimal disturbance. 

Soil specimens were carefully hand-trimmed to a diameter of approximately 1.4 
inches and a height of about 3.0 inches using a trimming device and wire saw. Water 
content was determined using the trimmings. The trimmed specimen was weighed and 
dimensions were measured before placing the specimen on the base pedestal of the 


RC/TS equipment. 


Next, two small circular pieces of filter paper were placed over the drainage holes 


on the base pedestal. The base pedestal was then secured to the base plate with four 
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screws. The specimen was then set on the base platen and the top cap was gently placed 


on top of specimen. Two layers of latex membrane were then placed over the soil 
specimen and each layer of membrane was sealed separately with O-rings on the top cap 
and the base pedestal. A silicon oil bath was then placed around the specimen to 
minimize the diffusion of air through the membrane. 

At this point the drive plate was attached to the top platen and the cell was 
positioned around the specimen and drive plate. The confining pressure was then applied 
and soil specimen was allowed to consolidate. When the specimen ake the end of 


primary consolidation, RC/TS testing was performed. 


9.3.2 RC/TS Testing Sequence 

A series of RC and TS tests were performed on soil specimens at each confining 
pressure. In this study, soil specimens were consolidated to the estimated mean in situ 
effective confining pressure (o'm) and two additional higher pressures at levels of 20n 
and 40'a. The RC and TS test are separated into two categories; low-amplitude tests and 
high-amplitude tests. The low-amplitude tests are tests at strains at which the soil 
specimen behaves linearly. The high amplitude tests are tests at strains where the soil 
acts nonlinearly. This nonlinearity can be identified from a decrease in resonant 
frequency in the RC test. Typically, the strain of a high-amplitude test is above 10°%. 

In this study, the TS testing was performed using loading frequency of 0.333 Hz 
for 5 cycles. The TS testing sequence is presented in Figure 9.2. For each confining 


pressure the TS tests were always performed first because the soil specimen was 
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--— Low-amplitude RC test 
—— High-amplitude TS test 


Time 


Shearing Strain 


Figure 9.2 Testing sequence for the TS test 


subjected to fewer cycles than in the RC test. However, low-amplitude RC tests were 
always performed before any high-amplitude TS tests to monitor any change in the small 
strain shear modulus (Gmax) of the soil specimen. If the Gmax changed significantly, the 
testing was paused until soil specimen returned to the original Gmax- 

After finishing the last high-amplitude TS test, a low amplitude RC test was 
Berane The RC test sequence is presented in Figure 9.3. Since the soil was subjected 
to high strain during high-amplitude TS test, a significant decrease in the value of Gmax 

_ was always observed after the last high-amplitude TS test. A series of low-amplitude RC 
test were then performed over period of time to monitor the increase in Gmax of the soil 


specimen. Once the value of Gmax was recovered to its original value (typically after at 
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-——— Low-amplitude RC test 
——— High-amplitude RC test 


Shearing Strain 


Figure 9.3 Testing sequence for the RC test 


least 4 hours), the series of high strain amplitude RC tests were performed. The RC 
testing sequence was similar, to the TS sequence, with low amplitude RC tests performed 
before any high amplitude RC test. If the value of Gmax did not return to the original 
value after a reasonably period of time, the change in the value of Gmax was then recorded 
and the next high amplitude test was performed. 

For both RC and TS tests the shear strains from less than10~ to more than 107% 
were applied to the soil, depending upon the soil specimen stiffness. When the high 
amplitude RC testing was completed the confining pressure was increased to the next 
level, the specimen was allowed to consolidated and the RC and TS testing sequence was 


repeated. 
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9.4 RC/TS ANALYSIS PROCEDURE 


9.4.1 Correction for Equipment Generated Damping 

The electromagnetic model presented in this study was used to predict the 
equipment generated damping, Deg, generated in the RC/TS drive system due to the back 
emf effect. As presented in Chapter 7, the back emf effect is modeled as an additional 
spring and an additional dashpot creating additional stiffness and viscosity in the system. 
The additional stiffness, called the equipment generated spring stiffness, Keq, is considered 
to be very small and can be neglected. The additional equipment generated viscous 
damping coefficient, Ceq, was used to calculate Deg. The ceq value can be predicted using 
Equation 5.31, and the model parameters; Le, Re, ki, kg. For the RC test, the value of Deg 


is calculated from: 
Da =. (9.1) 


As discussed in Chapter 7, the value of Ceg is reasonably constant over the a 
frequency range of 0-130 Hz, with a slight decrease at higher frequencies. The critical 
damping, c, is a function of resonant frequency and the polar moment of inertia of the 
drive plate, J, and top cap, J. Consequently, the variation of predicted Deg with resonant 
frequency mostly depends upon variation in critical damping, ce.. In this study, a 
variation of J, with frequency was observed due to the non-rigid body motion of the drive 
plate, as presented in Chapter 4. This J, variation was taken into account when 
calculating ce. Decreasing values of Deg with increasing resonant frequency was 


observed, as presented in Chapter 7. 
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In a series of RC test, the resonant frequency, fm at each stage of the test was used 


to determine the frequency dependent parameters, kg and Jo. Values of Ceq and Ce were 


then calculated using: 


and (9.2) 
c, = 4nZIf, (9.3) 


where f =—_—*—_., and 


" Wi-2D? © 


ZJ = total mass polar moment of inertia of the added mass (J,+J;). 

The model parameters, ki, Re and L, have constant values of 0.318 Ib-ft/rad, 45.6 
ohms, and 0.0256 H, respectively. Values of kg and J, vary with frequency. An example 
of the calculations for correcting D on specimen LGN14 is presented in Table 9.3. 

Increasing strain amplitude in RC tests results in decreasing resonant frequencies 
as can be seen in Table 9.3. The Deg generated in the RC/TS drive system is inversely 
proportional to the resonant frequency, so as the strain amplitude increases, the Deq also 
increases. Figure 9.4 shows a plot of the predicted Deq versus strain for the RC test of 
specimen LGN14. This soil specimen is very soft with low-strain resonant frequency of 
about 12 Hz, and a low-strain material damping of about 2%. As shown in Table 9.3 and 
Figure 9.4, the calculated Deg ranges from about 0.8% at low-strain amplitudes up to as 
high as 1.3% for high-strain amplitudes. The value of Deq is almost 40% of the actual — 
material damping at low-strains. As a result, correcting for the equipment generated 
damping, Deg, is critical when performing low-strain RC test on low-resonant frequency 


soils. 
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Table 9.3 Example of Calculations to Correct D for Equipment Generated Damping, 
Deg, for RC Test of Soil Specimen LGN14 from Logan, Utah . 


Resonant 


F Predicted | Calculated | Predicted | Corrected 
requency, 
fn 


(V/rad/s) 


(Hz 
1.286e-3 3.061 2.102e-3 | 0.3778 | 2.630e-3 | 3.178e- 0.827 2.234 
0.3778 | 2.630e-3 | 3.176e-1 


. 1 
| 2.102e-3 | 0.3777 | 2.631e-3 | 2.737e-1 | 0.961 | 6.151 | 
10.307 | 2.101e-3 | 0.3777 | 2.63le-3 | 2.479e-1 | 1.061 | 9.246 | 
15.855 
20.595 
Note: 1. J, = 0.0021013+4.673e-10freq”?™ 
2. kg = 0.3776+3.1227e-7freq? >” 


11.95 | 2.608¢-3 2.102e-3 | 0. 0.828 
6. 


1.4 


1.2 


1.0 


0.8 


Da %) 


0.6 


Sample ID: LGN14 
Confining Pressure: 6.2 psi 
Resonant Frequency Range: 8-12 Hz 


0.4 


0.2 


9.0 2 3 4 5 6789 2 3 4 56789 
0.001 0.01 0.1 


y (%) 
Figure 9.4 Predicted Deq for a series of RC test on soil specimen LGN14 
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Figure 9.5 shows the plot of measured D and corrected D versus strain amplitude 


for specimen LGN4 at confining pressure of 6.2 psi. The damping was corrected by 
subtracting the predicted Deq from the measured D. This procedure was used to account 


for the equipment generated damping for all of the RC tests throughout this study. 


9.4.2 Application of Stress Integration Approach 

The stress integration approach developed in this study was applied to all of the 
RC/TS tests on Bonneville clay. This approach was developed to account for the 
nonuniform stresses and strains occurring over the radius of the soil specimen. The basic 
concept of this approach is to find a t—y relationship described by an assumed soil model 
that results in a theoretical T-9 relationships that corresponds to the measured T-0 


relationship using curve fitting techniques. Details of this approach are presented in 


Chapter 8. 
Sample ID: LGN14 
Confining Pressure: 6.2 psi 
S © Corrected D 
a @ Measured D 
§ 
A 


2 3 4 5 6789 2 3 4 5 6789 
0.001 0.01 0.1 
y (%) 


Figure 9.5 Measured and corrected damping for a series of RC test on soil specimen 
LGN14 obtained from Logan, Utah at a depth of 14 feet 
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The stress integration approach was used to develop modulus reduction curves 


and damping curves using the following procedure. (1) Perform curve fitting on the 
measured T- O relation using either the closed form solution or numerical integration of 
each soil models following procedures presented in Chapter 8. (2) Select the best soil 
model that provides the best fit to the experimental data. (3) Develop the modulus 
reduction curves based on the selected soil model and the model parameters. (4) 
Determine values of Reg based on the hyperbolic model parameters over the range of 
rotations generating in testing. (5) Develop the damping curve by calculating the shear 


strain associated with each value of 9 using the appropriate Reg value. 


9.4.2.1 Curve Fitting for T-6 Relationships and Model Selection. The measured 
T-6 relationship from RC and TS test on seven soil samples at different confining 
pressures are presented in Figures 9.6-9.18. Curve fitting was performed on each 
measured T-O relationship using three soil models; the hyperbolic, the modified 
hyperbolic, and the Ramberg-Osgood. The closed form solution was used to predict the 
T-0 relationship for the hyperbolic model while the numerical integration was used to 
predict the T-9 relationship for the modified hyperbolic and the Ramberg-Osgood 
models. 

The small strain shear modulus, Gmax was evaluated from low-strain RC/TS tests. 
The Gmax value for each test was known, therefore it was constrained during the curve 
fitting. The best fit T-O relationships using the three soil models are also plotted in 


Figures 9.6-9.18. The model parameters and chi-square values determined from the 
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curve fitting are presented in Table 9.4. The shaded fields in Table 9.4 are the model 


parameter providing the best fit. 

In Figures 9.6-9.18, the measured and calculated values of T are plotted versus 
log 0 so the curve fitting can be evaluated over a wide range of strains. The values of 
Gmax Were constrained, therefore the curve fitting T-6 relationships always converge at 
low strains. All of the assumed soil models perform very well in the linear range 
verifying that the value Gmax used for curve fitting is close to the actual value of Gmax. At 
high strains, the stress integration approach performs effectively when the data points of 
the T-6 relationship are separated by small gaps. This can be observed by comparing the 
curve fitting results between specimen PRV93, in Figure 9.9, were the data points are 
separated by small gaps and specimen NBL24, Figures 9.11a and 9.11b, the data points 
are separated by large gaps. In Figures 9.11a and 9.11b, the modified hyperbolic model 
and Ramberg-Osgood performed equally well in fitting the last two points at high strains. 
However, in the region with no data points, the match is poor. The spacing between data 
points is a concern when applying the stress integration approach. 

Once the curve fitting was performed for each assumed soil model, the best fit 
model was selected to represent the t—y relationship for the soil. As discussed in Chapter 
8, chi-square is one indicator of how well the curve fitting matched the data. However, 
the chi-square value was always considered along with visual justification. The chi- 
square should be used alone only in the cases where all of assumed models match the 


measured T-@ relationship such that the differences cannot be distinguished visually. In 
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this study, the best fit soil model was selected for each test individually after careful 


consideration. 

The modified hyperbolic model fit was the mest commonly selected model for 
predicting the t~y relationships for Bonneville clay in this study. The Ramberg-Osgood 
model generally fit the T- relationship at high strains better than the other two models. 
Hne Ramberg-Osgood model fit the data poorly at strain levels close to the 
elastic threshold strain, y. The hyperbolic model fit the T-0 relationship well at low to 
medium strain levels but typically fit the data poorly at high strain levels. 


9.4.2.2 Development of Modulus Reduction Curve. The modulus reduction curve 


(G versus log y), and normalized modulus reduction curve (G/Gmax versus log y) were 
developed from the t~y relationship that provided the best fit to the measured T-6 data. 
The modulus reduction curve and normalized modulus reduction curve are generated as 
continuous functions. They can be generated for over any range of strains. However, it 
is only rational to develop the curve for strain levels in the range where the T-O 
relationship was measured. In this study, the modulus reduction curve was generated for 
each RC and TS test for the range of strain levels up to the strain level performed in each 
test which was corresponded to the maximum measured rotation, Omax measured for each 
test. The equivalent radius approach was used to determine maximum measured strain, 


Ymax, corresponded to that value of Omax. The Ymax, of each test is calculated from: 


Ymax =R (9.4) 
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Resonant Column Test 
Sample ID: PR V37 


A Measured T-0 


Torque (lb-ft) 


6 (rad) 
a) Confining pressure of 12.4 psi 
0.14 
a 0.12 © Measured T-9 
& 0.10 
2. 0.08 
A . 
e 0.06 
0.04 
0.02 
0.00 
10 


2 3 4 56789 2 3 4 56789 2 


10 ‘ais 10° 
b) Confining pressure of 24.8 psi 

0.35 

0.30 

0.25 A Measured T-0 


lb-ft) 


= 
ked 


0.20 -—— Hyperbolic Model 
0.15 ~——- Modified Hyperbolic Model 


0.10 
0.05 
0.00 


Torque 


——— Ramberg-Osgood Model 


9 (rad) 
c) Confining pressure of 49.6 psi 


Figure 9.6 Curve fitting T-@ relationships for the RC test on soil specimen PRV37 
obtained near the Provo River at a depth of 37 ft 
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10 
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> 
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Figure 9.7 Curve fitting for T-9 relationship for the RC tests on soil specimen PRV67 
obtained near the Provo River at a depth of 67 ft 


Reproduced with permission of the copyright owner. Further reproduction prohibited without permission. 


214 


0.25 
0.20 Torsional Shear Test 
des Sample ID: PRV67 
G 
2 0.15 @ Measured T-0 
3 
g 0.10 
© 
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0.05 
0.00 


6 (rad) y 
a) Confining pressure of 20.7 psi 


®& Measured T-6 
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10° 10° 10° 
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b) Confining pressure of 41.4 psi 


4 Measured T-0 


—— Hyperbolic Model 
—— Modified Hyperbolic Model 


Torque (lb-ft) 


~—— Ramberg-Osgood Model 


c) Confining pressure of 82.8 psi 


Figure 9.8 Curve fitting for T-9 relationship for the TS tests on soil specimen PRV67 
obtained near the Provo River at a depth of 67 ft 
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Figure 9.9 Curve fitting for T-0 relationship for the RC tests on soil specimen PRV93 
obtained near the Provo River at a depth of 93 ft 


Reproduced with permission of the copyright owner. Further reproduction prohibited without permission. 


216 


Torsional Shear Test 
Sample ID: PRV93 
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Figure 9.10 Curve fitting for T-O relationship for the TS tests on soil specimen PRV93 
obtained near the Provo River at a depth of 93 ft 
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Figure 9.11 Curve fitting for T-O relationship for the RC test on soil specimen NBL24 
obtained from Nibley at a depth of 24 ft 
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Torsional Shear Test 
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Figure 9.12 Curve fitting for T-0 relationship for the TS tests on soil specimen NBL24 
obtained from Nibley at a depth of 24 ft 
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Figure 9.13 Curve fitting for T-9 relationship for the RC tests on soil specimen LGN14 
obtained from Logan at a depth of 14 ft 
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Figure 9.14 Curve fitting for T-0 relationship for the TS tests on soil specimen LEONIS 
obtained from Logan at a depth of 14 ft 
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Figure 9.15 Curve fitting for T-0 relationship for the RC tests on soil specimen SLC14 
obtained near Salt Lake City Airport at a depth of 14 ft 
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Figure 9.16 Curve fitting for T-O relationship for the TS tests on soil specimen SLC14 
obtained near Salt Lake City Airport at a depth of 14 ft 
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Figure 9.17 Curve fitting for T-O relationship for the RC tests on soil specimen SLC14 
obtained near Salt Lake City Airport at a depth of 35 ft 
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Figure 9.18 Curve fitting for T-O relationship for the TS tests on soil specimen SLC14 
obtained near Salt Lake City Airport at a depth of 35 ft 
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The parameter y, from the hyperbolic model was used to obtain a normalized rotation 


(Omax/9,). Thus, the value of Req is obtained from the normalized Re, curve based on 
shear modulus developed in this study presented in Figure 8.13b. The modulus reduction 
curve at strains above Ymax is uncertain. 

The modulus reduction curve and normalized modulus reduction curve from RC 
and TS tests on each soil specimen‘are presented in Figures 9.19-9.25. The modulus 
reduction curves are plotted versus strain level up to the maximum measured strain. The 
curve was then extended to higher strains using the dashed lines. 

The results from RC and TS test were compared for all soil specimens except for 
specimen PRV37. On this sample, only the RC test was performed. It should be noted 
that the shear moduli obtained in the TS test, were evaluated from the first cycle of the 
hysteresis loop. For all of the specimens except for specimen NBL24, the shear modulus 
from TS test is higher than the shear modulus from the RC test. The two differences 
between the RC and TS tests are the frequency and the number of loading cycles. The TS 
tests were e at a frequency of 0.333 Hz for five cycles while the RC test was 
performed at resonant frequency with a high number of cycles. In general, the stiffness 
of a cohesive soil increases as the frequency increases. On the other hand, stiffness 
decreases as the number of cycles increases (Kim 1991). The lower values of shear 
modulus in RC test maybe explained by the effects of cyclic degradation outweighing the 


frequency effects. 
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Figure 9.19 Modulus reduction curves and normalized modulus reduction curves 
measured at various confining pressures for soil specimen PRV37 near the Provo River at 
a depth of 37 ft 
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Figure 9.20 Modulus reduction curves and normalized modulus reduction curves 
measured at various confining pressures for soil specimen PRV67 near the Provo River at 
a depth of 67 ft 
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Figure 9.21 Modulus reduction curves and normalized modulus reduction curves 
measured at various confining pressures for soil specimen PRV93 near the Provo River at 
a depth of 93 ft 
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Figure 9.22 Modulus reduction curves and normalized modulus reduction curves 
measured at various confining pressures for soil specimen NBL24 at Nibley at a depth of 
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Figure 9.23 Modulus reduction curves and normalized modulus reduction curves 
measured at various confining pressures for soil specimen LGN14 at Logan at a depth of 
14 ft 
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Figure 9.24 Modulus reduction curves and normalized modulus reduction curves 
measured at various confining pressures for soil specimen SLC14 near Salt Lake City 
Airport at a depth of 14 ft 
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Figure 9.25 Modulus reduction curves and normalized modulus reduction curves 
measured at various confining pressures for soil specimen SLC35 near Salt Lake City 
Airport at a depth of 35 ft 
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The effect of confining pressure can be investigated using the normalized 


modulus reduction curves. In Figures 9.19-9.25, the normalized modulus reduction 
curves show small effects of confining pressure with the curves falling in a very narrow 
band especially below the threshold strain yé. 

9.4.2.3 Development of Damping Curve. Material damping of the soil was 
determined using the half power bandwidth method and the free vibration decay method 
for the RC tests, and hysteretic damping for TS tests. The damping, D, was measured 
with an associated 6 for each method. To develop damping curves, the modified Reg 
approach presented in Chapter 8 was used. The Reg for damping can be obtained for any 
given 0 using the curve presented in Figure 8.13b, if the hyperbolic reference strain, yr is 
known. 

Procedures for developing damping curves are as follows. (1) Determine the 
value of y, from hyperbolic model using the stress integration approach. (2) Calculate 
normalized quantity 6/0, where 0, is calculated from Equation 8.26. (3) Determine the 
value of Req from Figure 8.13b. (4) Calculate y for a given 9 using the Req value. Figure 
9.26 shows an example of the Req values used to develop the damping curve for the TS 
test on sample PRV93. 

In this study, the half power bandwidth method was used in the RC test for 
determining damping of soil in its linear range (Dmin). The free vibration decay method 
was used to determine the damping at high strains in the RC tests. Free vibration decay 
tests were only performed on soil specimens obtained from Logan and near Salt Lake 


City Airport. Figures 9.27-9.29 shows the material damping curves for specimens from 
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Figure 9.26 Example for the evaluation of Req values for generating material damping 
curves for specimen PRV93 

Provo and Nibley. The high strain damping of specimens from Provo and Nibley were 
obtained only from the TS test. The damping for specimen PRV37 is not presented 
because the TS test was not performed on this specimen. 

Figures 9.30-9.32 shows a comparison of material damping from RC and TS tests 
for specimens obtained from Logan and near the Salt Lake City Airport. The Dmin 
obtained from RC and TS tests match well. However, the high strain damping measured 
in RC tests is higher than the damping measured in the TS tests. Typically, material 
damping for cohesive soil is independent of number of cycles but is significantly effect 


by the frequency (Kim, 1991). The RC and TS tests are performed at different 
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Figure 9.27 Damping curves measured on soil specimen PRV67 obtained near the Provo 


River at a depth of 67 ft 
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Figure 9.28 Damping curves measured on soil specimen PRV93 obtained near the Provo 
River at a depth of 93 ft 
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Figure 9.29 Damping curves measured on soil specimen NBL24 obtained from Nibley 
at a depth of 24 ft 
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Figure 9.30 Damping curves measured on soil specimen LGN14 obtained from Logan at 
a depth of 14 ft 
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Figure 9.31 Damping curves measured on soil specimen SLC14 obtained near Salt Lake 
City Airport at a depth of 14 ft 
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Figure 9.32 Damping curves measured on soil specimen SLC35 obtained near Salt Lake 
City Airport at a depth of 35 ft 


Reproduced with permission of the copyright owner. Further reproduction prohibited without permission. 


240 
frequencies. In the RC test, damping is measured at resonant frequency (above 10 Hz for 


this study) while the damping in the TS tests was measured at 0.333 Hz for this study. 
Typically, the damping curve is not presented in terms of a normalized quantity. 
The effect of confining pressure can be observed from the typical damping curves as 
shown in Figures 9.27-9.32. In general, the damping value decreases as the confining 
pressure increases. In this study, the confining pressure has a very small effect on the 
damping of Bonneville clay. In most cases, the damping values measured at different 
confining pressures at the same strain level are similar. An example of this is specimen 


PRV93 presented in Figure 9.28. 


9.4.3, Conclusions 

An improvement was achieved in obtaining the accurate material damping of 
Bonneville clay using the electromagnetic model. The model was used to correct for 
equipment generated damping, Deg, due to the back emf effect. The back emf effect is 
critical when performing low-strain RC test on low-resonant frequency material. 

The stress integration approach performs effectively when there are small gaps 
between the data points of the measured T-6 relationship. The modified hyperbolic 
model was most commonly the best model used for developing stress-strain relationships 
of Bonneville clay. The Ramberg-Osgood model provided best fit of the T-9 relationship 
at high strains and poorly fit of the T-6 relationship at low strains. The hyperbolic model 
fit the T- relationship well at low to medium strains by typically fit the data poorly at 


high strains. 
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9.5 DISCUSSION ON FACTORS AFFECTING 


MODULUS REDUCTION AND DAMPING 

CURVES 

This section contains a discussion of the effects of plasticity index and confining 
pressure on the normalized modulus reduction and damping curves. The normalized 
modulus reduction curves and damping curves are compared with generic curves used in 


practice for a given plasticity index and confining pressure. In addition, effect of cyclic 


degradation on Bonneville clay is also discussed. 


9.5.1 Plasticity Index Effect 

An important soil property affecting the nonlinear stress-strain behavior of clay is 
plasticity index (PI). The seven Bonnville clay specimens used in this study have PI’s 
ranging from 2 to 34 as shown in Table 9.2. The effect of PI can be evaluated using 
normalized modulus reduction curves for all soil specimens at their in situ mean effective 
confining pressure presented in Figure 9.33. 

The normalized modulus reduction curves presented in Figure 9.33 were 
determined from RC and TS tests. Soil specimens used in this study can be divided into 
four groups: samples with PI’s of 34, samples with PI’s of 17-20, samples with PI’s of 
10-11, and samples with PI’s of 2. In general, soil behaves more linearly (the elastic 
threshold strain, yé, increases) as the PI increases. This general trend of PI with soil 
linearity is not clearly observed in the RC test results. For the TS test results, the samples 
with PI of 34 behave more linear than the samples with PI’s of 10-17. An outlier was 


observed in sample with a PI of 11 that behaves as linearly as the samples with a PI of 34. 
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Figure 9.33 Variation of normalized modulus reduction curves with PI for RC and TS 
tests on Bonneville clay 
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Figure 9.34 Variation of damping with PI i t and TS tests on Bonneville clay 
The sample with a PI of 2 was expected to behave the most nonlinearly but, the results 
show that this sample behaves similarly to the samples with PI’s of 10-20. This might be 
because this sample has the interbeded layers of silt and clay. The clay layers may have 
more effect on damping than the silt layers. 

The variation of damping with PI for the Bonneville clay tested in this study is 
presented in Figure 9.34. In general, the damping decreases as PI increases. The 
measured damping in TS test on specimen SLC35 with PI of 2 has the highest damping 
value. In the TS tests, the damping in specimens LGN14 and NBL24 with PI’s of 34 
show the lowest values of damping. Therefore, the damping behavior of Bonneville clay 


is typical of cohesive soils. 
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Empirical curves used to estimate modulus reduction and damping curves based 


upon PI where compared with the results from RC and TS tests on Bonneville clay. The 
empirical curves used are curves by Sun et al. (1988), Vucetic and Dobry (1991), and 
Darendeli (2001). 

Sun et al. (1988) proposed empirical modulus reduction and damping curves for 
clays with different ranges of PI’s. Figure 9.35 and 9.36 show the comparison for the 
curves measured in this study and the curves proposed by Sun et al. (1988). In Figure 
9.35, three empirical curves are plotted for PI ranges of 5-10, 10-20, and 20-40. 
Bonneville clay specimens with PI’s of 2-20 plot close to the Sun et al (1988) curve for a 
PI range of 10-20. The modulus reduction curves for samples with PI of 34 are in 
reasonably agreement with the Sun et al curve for a PI range of 20-40. The empirical 
curve for a PI range of 5-10 underestimates the measured modulus for the soil specimen 
with a PI of 2. 

Sun et al. (1988) proposed the empirical damping curves with no correlation to PI 
as shown in Figure 9.36. The damping curve of Bonneville clay from TS tests plot close 
to the average curve at medium strains and close to the upper bound at high strains. 
While the small strain damping curve plots closer to the lower bound curve. However, 
these empirical damping curves span such a wide band that they are of little practical 


value. 
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b) Comparison with TS test results 


Figure 9.35 Comparison of normalized modulus reduction curves of Bonneville clay 
with empirical curves proposed by Sun et al. (1988) 
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Figure 9.36 Comparison of damping curves of Bonneville clay with empirical curves 
proposed by Sun et al. (1988) 

Comparisons of modulus reduction curves and damping curves to the empirical 
curves proposed by Vucetic and Dobry (1991) are presented in Figures 9.37-9.38. The 
Vucetic and Dobry (1991) curves for PI of 2, 10, 20, and 34 were generated by 
interpolation from their original proposed curves presented in Chapter 2. In Figure 9.37, 
the Vucetic and Dobry (1991) modulus reduction curve for PI of 34 shows good 
agreement with the curves of Bonneville clays with the same PI value. The Vucetic and 
Dobry (1991) curves for PI of 2 and 10 exhibit more nonlinearity than the actual behavior 


of Bonneville clay at the same values of PI. 
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Figure 9.37 Comparison of normalized modulus reduction curves of Bonneville clay 
with empirical curves proposed by Vucetic and Dobry (1991) 
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Figure 9.38 Comparison of damping curves for Bonneville clay with empirical curves 
proposed by Vucetic and Dobry (1991) 


Damping curves proposed by Vucetic and Dobry (1991) show higher values of 
damping at medium strain levels (10°-10°%) than the damping curves for Bonneville 
clay presented in Figure 9.38. The Vucetic and Dobry (1991) damping curves also show 
a more pronounced effect of PI than was measured on the Bonneville clay. The Vucetic 
and Dobry (1991) damping curves have been left undefined at low strains. At high 
strains, the Vucetic and Dobry (1991) damping curves underestimates the Bonneville clay 
damping curves. The measured damping for the soil sample with a PI of 2 matches well 


with the Vucetic and Dobry (1991) curve for a PI of 2 at strains above 0.1%. 
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Darendeli (2001) proposed a model to predict the modulus reduction curves and 


damping curves accounting for the effect of PI, confining pressure, OCR, frequency 
(freq), and number of cycles (N). The modulus reduction curves and damping curves 
were generated using the proposed model by Darendeli (2001) for PI values of 2, 10, 20, 
and 34 with an OCR = 1, on = 1 atm, freq = 0.333 Hz, and N = 5 cycles. The Darendeli 
(2001) curves were compared with Bonneville clay curves in Figures 9.39-9.40. 

The Darendeli (2001) modulus reduction curves show a smaller PI effect than the other 
empirical curves. In Figure 9.39, all of the Darendeli (2001) curves for different PI plot 
below the curve for Bonneville clay. This indicates that Bonneville clay behaves more 
linearly than predicted by the Darendeli (2001). It is interesting to note that the width of 
the band for effect of PI value of of 2 up to 34 predicted by Darendeli (2001) is close to 
the bandwidth of the measured curves. 

Damping curves proposed by Darendeli (2001) are compared to the measured 
damping in Figure 9.40. Again, the damping curves proposed by Darendeli (2001) shows 
narrow band of PI effect. Damping from RC test was higher than that predicted damping 
by Darendeli (2001). However, the damping from TS tests matches the Darendel (2001) 
damping curve at high strains. The Darendeli (2001) damping curve slightly 
overestimates the damping at medium strains. This is consistent with the more linear 
behavior of Bonneville clay. The Darrendeli (2001) damping curve for a PI of 2 


underpredicts the damping of the Bonneville clay at high strains. 
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Figure 9.39 Comparison of normalized modulus reduction curves of Bonneville clay 
with empirical curves proposed by Darendeli (2001) 
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Figure 9.40 Comparison of damping curves for Bonneville clay with empirical curves 
proposed by Darendeli (2001) 

Another model to predict the modulus reduction curves and damping proposed by 
Ishibashi and Zhang (1993). Ishibashi and Zhang (1993) developed a set of equations to 
generate the modulus reduction and damping curves for a given values of PI and on’. 
However, the equations proposed by Ishibashi and Zhang (1993) shows unrealistic 
modulus reduction and damping curves at high strain levels for the range of PI between 
10-40. A discussion about Ishibashi and Zhang (1993) curves is presented in Chapter 2. 
In this study, comparisons are not made between Bonneville clay Ishibashi and Zhang 


(1993). 
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9.5.2 Confining Pressure Effect 


In this section, the effects of confining pressure on the modulus reduction curves 

and damping curves for Bonneville clay is compared with the empirical curves developed 

_by Darendeli (2001). The empirical curves are compared with the results from TS tests 
on soil specimen PRV93. Three different confining pressures were applied to this 
specimens during testing. The sample PRV93 was confined at mean effective stresses of 
28, 56, and 112 psi. The modulus reduction and damping curves of the specimen PRV93 
aré presented in Figures 9.21 and 9.28, respectively. 

The normalized reduction and damping curves for sample PRV93 are compared 
with the Darendeli (2001) empirical curves presented in Figure 9.41. The normalized 
reduction curve becomes more linear as the confining pressure increases. Damping 
curves of specimen PRV93 shows a very small effect of the confining pressure. 

The Darendeli (2001) curves in Figure 9.41 were regenerated for the same 
confining pressure at which each specimen was subjected to in the RC and TS tests. In 
Figure 9.41, the measured normalized modulus reduction curve for sample PRV93 plots 
above the predicted curves. Less effect of confining pressure is seen in the normalized 
reduction curves of PRV93 than that predicted by the Darendeli (2001). For damping, 
the Darendeli (2001) curve plots above the measured curve at medium strains. With the 
normalized modulus reduction curve, the effect of confining pressure in damping for 


PRV93 is less than that predicted by Darendeli (2001). 
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Figure 9.41 Comparison of confining pressure effects from TS tests on specimen PRV93 


and curves proposed by Darendeli (2001) 
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9.5.3 Cyclic Degradation 


At high strain levels, the stiffness of soil is reduced due to cyclic degradation. 
This is observed as the slope of hysteresis loop begins to flatter with increasing number 
of cycles. This means that the secant shear modulus decreases as the numbers of cycle 
increases. 

The shear modulus and damping in this study was determined using the first 
loading cycle. However, the cyclic degradation of Bonneville clay was observed for 
most of the soil specimens at high strain levels. Figure 9.42 shows an example of cyclic 
degradation for sample PRV67 test at a confining pressure of 42 psi for strain level of 
0.06-0.22%. These hysteresis loops are plotted relative to T-0 measured directly from TS 
tests. These was little cyclic degradation in tests performed at strains of 0.06% and 
0.08% as shown in Figures 9.42a and 9.42b. However, increasing degradation occurred 
at strains of 0.12% and 0.22% as shown in Figures 9.42c and 9.42d. 

To quantify the effect of cyclic degradation, the slope of the T-9 relation for each 
cycle was determined. The stress integration approach was then applied to determine the 
modulus for each cycle. For this specimen PRV67, the modified hyperbolic model was 
used to describe the t—-y relationship. The normalized reduction curves for specimen 
PRV67 for each cycle are presented in Figure 9.43. Model parameters for the modified 
hyperbolic model for each cycle are in Figure 9.43. The value of y, decreases as the 
number of cycle increases which is expected because of the decreasing shear modulus. 
The change of parameter, a, does not show a consistent trend in this case however, the 


value of the first cycle was observed to be the lowest. 
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Figure 9.42 Hysteresis loops of five cycles of loading on soil specimen PRV67 obtained 
near the Provo River at a depth of 67 ft tested at a confining pressure of 42 psi 


From Figure 9.43, two threshold strains were determined. These are the elastic 


threshold strain, y and the cyclic threshold strain, y€. The y is defined earlier as the 


strain level where shear modulus begins to be affected by strain amplitude. The yê is 
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Figure 9.43 Variation of normalized modulus reduction curves with number of cycles 
for specimen PRV67 
defined as the strain level at which the shear modulus begins to be affected by the number 
of loading cycles. As shown in Figure 9.43, the y;° was selected at strains where the shear 
modulus was equal to 0.98 of Gmax and the y,° was defined at strain levels at which the 
difference between the shear modulus of first and fifth cycle is equal to 2% (Kim 1991). 
For this sample, the y€ and y were equal to 0.0006683% and 0.029853%. The different 
in shear modulus at first and fifth cycles was up to 8% at a strain of 0.22%. 

The effect of cyclic degradation on damping was also investigated. The damping 
was determined from hysteresis loops presented in Figure 9.42 for each cycle of TS 


loading. The variation in damping with strain level for each cycles is presented in Figure 


Reproduced with permission of the copyright owner. Further reproduction prohibited without permission. 


257 
9.44. The damping for each cycle is identical for all loading cycle at strains below the 


cyclic threshold, y,°. At strain above y;’, the first cycle of loadings exhibited lower 
damping than subsequent cycles. These are consistent with typical cohesive soils 
reported by Kim (1991). 

Because the number of cycles in this study was limited to five cycles, more study 
would be required to further investigate the effect of cyclic degradation. Kim (1991) 
investigated cyclic degradation using up to one hundred cycles and his results show a 


more pronounced reduction of shear modulus with large number of cycles. 
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Figure 9.44 Variation of damping for different numbers of loading cycle 
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9.5.4 Conclusions 


The overall behavior of Bonneville clay was found to behave more linearly than 
was predicted by the empirical curves for behavior of typical clay at the same values of © 
PI. Less confining pressure effect was also observed for Bonneville clay. The empirical 
curves proposed by Darendeli (2001) provide the best match wrth the behavior of 
Bonneville clay. Cyclic degradation of Bonneville clay requires more study. 

The more linear behavior of the Bonneville clay may be explained by geologic 
age effect. Bonneville clay is a relatively young soil, deposited 12,000-20,000 years ago 
during the Quaternary period. Zhang et al. (2004) compared modulus reduction curves 
from three groups of soils: Quaternary, Tertiary and older, and residua/saprolite. They 
found that Quaternary soils exhibit more linearity than soils of the other two groups. 
They also observed that most of the soils used to develop the Darendeli (2001) model are 
Tertiary soils. Accordingly, Zhang et al. (2004) showed a better agreement between the 
Darendeli (2001) model and Tertiary age soils than the Darendeli (2001) model and 


Quaternary soils. 
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CHAPTER 10 


SUMMARY, CONCLUSIONS AND RECOMMENDATIONS 


10.1 BACKGROUND 

A goal of this study was to achieve the accurate measurement of nonlinear stress- 
strain behavior for soil using the combined resonant column and torsional shear (RC/TS) 
testing device. Two weaknesses of the RC/TS tests are the back electromotive force 
(back emf) created in the coil and magnet drive system of the RC/TS apparatus and 
nonuniform stress-strain occurring over the radius of the soil specimen in the RC/TS test. 
This study focused on developing methods to deal with these problems more precisely. 
The methods developed in this study are the use of electromagnetic model to describe and 
predict the back emf effect, and the stress integration approach to account for the 
‘nonuniform stress-strain in the soil specimen. 

The electromagnetic model was developed by combining the electrical, magnetic, 
and mechanical properties of the RC/TS apparatus. The electrical part of model describes 
the current flowing in the drive coil circuit. The magnetic part of the model describes the 
development of torque and back emf in the coil-magnet system. The mechanical part of 
the model describes the RC/TS drive system and soil specimen as a single-degree-of- 
freedom oscillater. All independent measurement was made of each model parameters. 
Each model parameter was measured as a function of frequency, when possible, to 
identify and account for the frequency dependent parameters. The model was verified by 
comparing open and closed circuit measurements. And procedures were developed to 


apply the electromagnetic model to both the RC and TS tests. 
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Previous studies have developed an suggested the use of equivalent radius 


approach to account for the nonuniform stress-strain in soil in the RC/TS tests (Chen and 
Stokoe 1979). The equivalent radius approach proposed utilizes an equivalent radius 
ratio to calculate the stress and strain from the torque and rotation measured in RC/TS 
testing. This approach is adequate for the tests at low to medium strains however, it 
becomes less accurate for the tests at high strains. The stress integration approach was 
developed in this study to more precisely account for nonuniform stress-strain at any 
strain level. The stress integration approach is applied by integration using three different 
soil stress-strain models and to calculate a theoretical torque and rotation relationship. 
Curve fitting techniques are then applied to match the measured torque and rotation 
relationship with the theoretical torque and rotation relationship. The stress-strain model 
providing the best fit to the measured data is used to characterize the soil specimens. The 
stress integration approach was also used to account for nonuinform stress-strain in 


damping. 
10.2 SUMMARY AND CONCLUSIONS 


10.2.1 Polar Moment of Inertia of RC/TS Drive Plate 

The polar moment of inertia of the drive plate was calibrated in this study. 
Calibration specimens with different resonant frequency were used to investigate the 
moment of inertia of the drive plate at frequencies up to 200 Hz. An open circuit test was 
used to eliminate the back emf effect in the drive plate calibraiton. The drive plate 
calibration was performed on the USU drive plate and the UT drive plate No. 9 and 


results were compared. 
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A frequency effect on the polar moment of inertia of the drive plate was observed. 


Results showed that the polar moment of inertia of the USU drive plate increases 28% as 
the frequency increases from 14 to 188 Hz. This frequency effect in polar moment of 
inertia causes an error in shear wave velocity for up to 12% if not account for. The 
frequency effect on the polar moment of inertia of the UT drive plate was much smaller. 
The increasing polar moment of inertia of the drive plate with frequency is explained by 
non-rigid body motion of the drive plate. The non-rigid body motion of the drive plate 
was verified using a finite element model of the drive plate. 

The finite element analysis showed that stiffness of the connection of the drive 
plate and magnet plays an important role in frequency effect on polar moment of inertia. 
Different drive plates behave differently if they utilize different types of connections. 
Investigation of the polar moment of inertia of the drive plate with frequency should be 
performed to characterize any effect. To obtain high accuracy RC measurements, the 
shear wave velocity, the frequency effect on the polar moment of inertia of the drive plate 


must be taken into account. 


10.2.2 Electromagnetic Model 

10.2.2.1 Torque and Back Emf in the Drive System. The application of the 
principles of electromagnetism was used to describe the torque generated in RC/TS drive 
system and the back emf effect created due to the motion of magnets through the coils. 
In this study, the theoretical equations were developed to predict the electromagnetic 
properties of an ideal coil and magnet system. These are: 


k; = NI,Bd, (10.1) 
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k, =———_. 10.2 
are (10.2) 
where B = strength of magnetic field of the permanent magnet, 


N = total number of loops in the coil 

l = circumference length of coil loop, 

A. = cross section area of coil loop, 

y = width of coil loop in longitudinal section, and 

d= moment arm of the drive plate. 

10.2.2.2 Tranfer Function. The electromagnetic model was developed by 

combining the electrical, magnetic, and mechanical properties of to the RC/TS apparatus. 
The transfer function for voltage to rotation in the RC/TS tests is: 


k. 


1 


R, +joL, 


H, =2 = J> (10.3) 
V, lk, +k. )+ jalc, +C.4 =o J] 


o’kk,L, 


a RoT and (10.4) 


Cuq = TEn " (10.5) 
where H, = transfer function for voltage to rotation, 

6 = rotation, 

Va = measured voltage output from voltage source, 

ky = torsional spring stiffness, 


Ce = torisonal viscous damping coefficient, 


J = polar moment of inertia of the system, 
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Keg = equipment spring stiffness, 


Ceq = equipment viscous damping coefficient, 
R, = coil resistance, 

L. = coil inductance, 

ki = torque-current factor, 

kg = back emf-rotational velocity factor, and 
© = angular velocity (@ = 27f). 

This model describes where the back emf effect in terms of an equipment spring 
stiffness, keg, and an equipment viscous damping coefficient, Ceq. 

10,2.2.3 Model Parameters. There are four model parameters in the model 
presented in Equations 10.3-10.4. These are the coil resistance, Ro, the coil inductance, 
Le, the torque-current factor, kj, and the back emf-rotational velocity, kg. Independent 
measurements were made of the model parameters for the USU drive plate and the UT 
drive plate No.9. Where possible, the variation of each model parameter with frequency 
was evaluated. 

For the USU drive plate, a slight variation in R. and L. with frequency was 
measured, however, this variation is small and can be neglected. The value ofk; was 
measured in statically condition and compared with the back-calculated k; at higher 
frequency and it was found to be nearly constant at all frequencies. A crosstalk effect 
was observed during the measurement of kg. The cross talk interfered with the 
measurement of kp and resulted in unrealistic kg values. A measurement of crosstalk 
voltage was performed and the calculation of kg values was corrected to account for the 


crosstalk effect. An increasing of kg with frequency was measured. Additionally, there 
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was some scatter in the values measured with different calibration specimens. kg was 


characterized using an average value and upper and lower bounds. This uncertainty was 
found to have a minor effect on RC/TS results. 

For the UT drive plate No. 9, the variation of R. and Le with frequency showed 
trends similar to the USU drive plate. The measured value of ki was about twice the 
value measured for the USU drive plate. The value of kg was also about twice as high as 
the value for the USU drive plate with less variation with frequency than the USU drive 
plate. The higher values of kj and kg can be explained by stronger magnets used in the 
UT drive plate. 

Efficiency of torque and back emf was determined from the theoretical and 
measured k; and kg. The torque efficiency was found to be about 75-80% for both USU 
and UT drive plate. And, the back emf efficiency was about 32-50% for both drive plate. 
The lower back emf efficiency is due to the nonuniform magnetic field in the coil cross 
section. 

Based on the equations derived from the theoretical value of k; and kp, the 
capacity of the RC/TS drive system in generating torque can be improved by using the 
higher numbers of loops in the coil or stronger magnet. However, more back emf effects 
result from these modifications. 

10.2.2.4 Model Verification. In RC test, the increasing of resonant frequency, 
Afm, and the equipment generated damping, Deq, due to the back emf effect was 
quantified. The keq and Ceg were predicted using the measured model parameters and the 


predicted value of Afm and Deg were calculated. Open and closed circuit measurements 


Reproduced with permission of the copyright owner. Further reproduction prohibited without permission. 


265 
were performed on calibration specimens in the RC/TS apparatus to determine 


experimental values of Afm and Deg. 

The experimental values of Afm and Deq was used to calculate the measured Keq 
and Ceg. The measured and predicted keq increase with frequency in similar fashion. The 
model underprediction of ke, value at frequency between 60-130 Hz for the USU is of 
little concern because the effect of back emf on resonant frequency is very small and can 
be neglected. The slight decreasing value of the measured Ceq with frequency was 
observed for the USU and UT drive. The measured value of ceq scattered and reasonably 
constant over the frequency range of 0-130 Hz. The predicted Ceq value shows slightly 
underestimates the measured Ceg at higher frequency. This underprediction is explained 
by other sources of energy loss and the uncertainty in frequency effet on k; value. 

10.2.2.5 Model Application. fi RC test, the back emf effect on resonant 
Benay is very small and can be neglected. Effect of back emf on damping must be 
taken into account. The electromagnetic model with measured model parameters can be 
used to predict Ce and the equipment generated damping, Deg. The model predicts Deq as 
a function of resonant frequency and mass polar moment of inertia of the system. 
Therefore, the prediction of Deq can be performed for the RC testing on soil specimen 
that setup using any type of top platen. The Deg is then corrected with the measured 
damping. 

In TS test, the electromagnetic model verified that the effect of back emf on the 
torque measurement is very small and can be neglected. The discrepancy of equipment- 
generated damping in RC and TS test from previous work (Hwang 1997) was verified by 


the effect of phase shift between the voltage and current. 
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10.2.2.6 Model Limitation. The electromagnetic model prediction of Deg is 
limited for up to 200 Hz, however the Deg is very small above this frequency range. 
Therefore, the Deg at frequency above 200 Hz can be negligible. 

10.2.2.7_ Advantages of Electromagnetic Model. Three advantages of the 
electromagnetic model are summarized. First, the model successfully describes the 
behavior of back emf generated in the RC/TS equipment. Second, the model is applied 
for both the RC and TS tests. Third, the model can be applied for any typical RC/TS 


apparatus. 


10.2.3 Stress Integration Approach 

The stress integration approach is a new, more general approach to account for 
nonuniform stress-strain effect in RC/TS testing. The stress integration approach is 
applied by integration using three different soil stress-strain models to calculate a 
theoretical torque and rotation relationship. The stress-strain model providing the best fit 
of the theoretical torque and rotation relationship to the measured torque and rotation 
relationship is used to develop the stress-strain relationship of the soil specimen. The 
Three assumed model used for this approach are the hyperbolic model, the modified 
hyperbolic model, and Ramberg-Osgood model. The stress integration approach use 
either closed form integration or numerical integration to develop a theoretical torque and 
rotation relationship. The closed form integration was applied for the hyperbolic model 


and the closed form equation is: 


1 L L L ƏR 
T= Laaman e| 2R? -ardar |e 2 mart EP | ieri- a) . (10.6) 
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Numerical integration was applied for the modified hyperbolic model and 


Ramberg-Osgood model. 

The nonuniform stress-strain in damping is taken into account by using a 
modified equivalent radius approach. The conventional equivalent radius approach was 
modified to determine the strain corresponding to a measured damping value more 
precisely. Values of Req defined as ratio of the equivalent radius and the outside radius of 
soil specimen, were evaluated based on shear modulus and damping are determined 
separately. Results from this study are the plots of the Reg versus normalized rotation 
curve for shear modulus and damping shown in Figure 10.1. These Reg versus 
normalized rotation curves account for the soil nonlinearity. All curves converge to a 
value of about 0.8 at low strains. At high strains, the Reg curves based on damping are 
significantly lower than the Req curves based on shear modulus. Figure 10.1 suggests that 
using the same Req values for both shear modulus and damping is not appropriate. 

Advantages of using the stress integration approach are summarized as follows: 

1. The conventional equivalent radius approach uses single values of Reg for ranges 
of strain to calculate stress and strain. The actual Req decreases with an increasing 
of strains and the point where the Reg decreases is associated with the value of 
reference strain, yr. The stress integration approach effectively accounts for this 


variation of Req over the range of strains and the soil nonlinearity. 
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Figure 10.1 R. curves versus normalized rotation based on G and D 


2. The conventional equivalent radius approach was limited at strain level up to 10° 
1%. While the stress integration can be applied for higher strains. 

3. The conventional equivalent ae approach was developed based on the 
hyperbolic model. The stress integration approach is not limited to only one 
stress-strain soil model, therefore, the approach provide options to select the best 
model to match the experimental data. 

4. The conventional equivalent radius approach developed the value of Reg on shear 
modulus and this value is also applied for calculations of strain for damping. By 
applying the stress integration approach, the Reg value for damping was found to 
be different than the Req for shear modulus. The Req value for damping is less 
than the Reg ratio for shear modulus and the difference tends to increases as the 
strain increases. Therefore, the stress integration approach provides more 


accurate damping curves in RC/TS testing. 
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One limitation of the stress integration approach is that the method can only be 


applies with specific stress-strain models. In some special cases, no available stress- 
strain model provides a good best to the experimental data. 

Another limitation of the stress integration approach is that Req curve for damping 
was developed assuming that Masing behavior applies according to Masing rules. This 


assumption requires further study. 


10.2.4 Dynamic Properties of Bonneville Clay 

RC/TS testing of Bonneville clay was performed on seven soil specimens 
obtained from sites along the Wasatch front, Utah. High accuracy test results were 
obtained by applying the correction of equipment generated damping predicted by the 
electromagnetic model and the application of the stress integration approach. 

The effect of plasticity index and confining pressure effect on the nonlinear 
dynamic properties of Bonneville clay was investigated and compared with the empirical 
relationships developed by other researchers. Results showed that Bonneville clay 
exhibits more linearity than predicted using the empirical relationships. The more linear 
behavior of Bonneville clay might be explained by the aging effects. Effects of plasticity 
index and confining pressure on modulus reduction and damping curves observed for 


Bonneville clay was much smaller than that predicted by the empirical relationships. 


10.3 RECOMMENDATIONS 
Based on this work, the following issues require further investigations: 
1. An investigation of the effect of the non-rigid body motion of the drive plate on 


the polar moment of inertia in this study is only preliminary. Better modeling of 
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the connection between the drive plate and magnet requires further study. Other 


factors affecting the non-rigid body motion may need to be considered such as the 
accelerometer and counter weight. f 

2. An investigation into the effects of rotation on proximitor behavior is desirable. 
High accuracy measurement of small rotation angle is required to accomplish this 
task. 

3. An investigation for the variation of k; with frequency requires further study. A 
stiffer torque sensor might be used to measure torque at higher frequencies. 

4. The uncertainty in kg value needs to be investigated further. Explanations for the 
increasing of kg with frequency also requires further study. A more detailed study 
of magnetic effects may be required. 

5. Increasing capacity of torque generated in RC/TS equipment is recommended as a 
result of this study by the use of more number of loops in coil or stronger 
magnets. Alternatively, the connection of coil in parallel can also increase the 
torque capacity. When coils are in parallel, the voltage potential dropped across 
each coil is less than for the coils connected in series. This implies that for the 
same input voltage, higher torque can be generated in case of coils in parallel. 

6. Although the back emf effect at frequencies higher than 200 Hz is very small, the 
values of model parameters could be investigated at higher frequency. This woud 
require calibration specimens with resonant frequencies above 200 Hz. 

7. The problem with the available stress-strain soil models could not predict the soil 


behavior for the entire range of strain level requires further study. More versatile 
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soil model that can be fit the experimental data from very low strain level to very 


high strain level could be developed. 

8. Study of the non-Masing behavior in soil is required further study. Masing 
behavior was assumed in applying the stress integration approach to damping. 

9. The effect of plasticity index, confining pressure and cyclic degradation for 
Bonneville clay requires further study. Other factors affecting dynamic properties 
of Bonneville clay needs to be considered such as the water content, degree of 
saturation, cementation etc. More RC/TS testing on specimens with interbedded 


layers is recommended. 
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